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CONTROL OF WATER WAVES 


T. ALAZARD, P. BALDI, D. HAN-KWAN 


Abstract. We prove local exact controllability in arbitrary short time of the 
two-dimensional incompressible Euler equation with free surface, in the case with 
surface tension. This proves that one can generate arbitrary small amplitude 
periodic gravity-capillary water waves by blowing on a localized portion of the 
free surface of a liquid. 
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1. Introduction 

Water waves are disturbances of the free surface of a liquid. They are, in general, 
produced by the immersion of a solid body, the oscillation of a solid portion of 
the boundary or by impulsive pressures applied on the free surface. The question 
we address in this paper is the following: which waves can be generated from the 
rest position by a localized pressure distribution applied on the free surface. This 
question is strictly related to the generation of waves in a pneumatic wave maker 
(see [42, §21], [13]). Our main result asserts that, in arbitrarily small time, one can 
generate any small amplitude, two-dimensional, gravity-capillary water waves. This 
is a result from control theory. More precisely, this article is devoted to the study of 
the local exact controllability of the incompressible Euler equation with free surface. 

There are many known control results for linear or nonlinear equations (see the 
book of Coron [14]), including equations describing water waves in some asymp¬ 
totic regimes, like Benjamin-Ono ([33, 30]), KdV ([40, 31]) or nonlinear Schrodinger 
equation ([17]). In this paper, instead, we consider the full model, that is the in¬ 
compressible Euler equation with free surface. Two key properties of this equation 
are that it is quasi-linear (instead of semi-linear as Benjamin-Ono, KdV or NLS) 
and secondly it is not a partial differential equation but instead a pseudo-differential 

T.A. was partly supported by the grant “ANAE” ANR-13-BS01-0010-03. This research was 
carried out in the frame of Programme STAR, financially supported by UniNA and Compagnia di 
San Paolo; it was partially supported by the European Research Council under FP7, and PRIN 
2012 “Variational and perturbative aspects of nonlinear differential problems”. 
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equation, involving the Dirichlet-Neumann operator which is nonlocal and also de¬ 
pends nonlinearly on the unknown. As we explain later in this introduction, this 
requires to introduce new tools to prove the controllability. 

To our knowledge, this is the first control result for a quasi-linear wave equation 
relying on propagation of energy. In particular, using dispersive properties of gravity¬ 
capillary water waves (namely the infinite speed of propagation), we prove that, for 
any control domain, one can control the equation in arbitrarily small time intervals. 


1.1. Main result. We consider the dynamics of an incompressible fluid moving 
under the force of gravitation and surface tension. At time t, the fluid domain f l(t) 
has a rigid bottom and a free surface described by the equation y = rj(t,x), so that 

Q(t) = { (x, y) G M 2 ; -b < y < 77 (i, x) } , 

for some positive constant b (our result also holds in infinite depth, for b = oo). The 
Eulerian velocity field v is assumed to be irrotational. It follows that v = ’S7 X) y ( t> for 
some time-dependent potential (j> satisfying 

A x , y (j) = 0, d t (j) + i \V x , y 4>\ 2 + P + gy = 0, dy<f>\ y= - b = 0, (1.1) 

where g > 0 is the gravity acceleration, P is the pressure, V x , y = {d x ,d y ) and 
A Xj y = d x + dy. The water waves equations are given by two boundary conditions 
on the free surface: firstly 

dtJ) = \/l + (d x y) 2 d n cj)\y =v 

where d n is the outward normal derivative, so \Jl + ( d x r]) 2 d n <f> = d y 4> — (d x rj)d x 4>. 
Secondly, the balance of forces across the free surface reads 

P\y =r i — nH(r]) + P ex ^(t, x) 

where k is a positive constant, P e xt is an external source term and H(r]) is the 
curvature: 

H{v) := ^Oi+tL) 2 ) = PP0PW 

Following Zakharov [45] and Craig and Sulem [16], it is equivalent to work with 
the trace of <f> at the free boundary 


= 4>(t,x,r]{t,x)), 


and introduce the Dirichlet-Neumann operator G(r]) that relates ^ to the normal 
derivative d n (j) of the potential by 


(' G(rj)il>)(t,x ) = \/l + (d x rj) 2 d n (/)\ y=rt(t!x) . 

Hereafter the surface tension coefficient k is taken to be 1. Then ( 77 , i/j) solves 
(see [16]) the system 


' d t r) = 

l 

d t i> + gr) + ~{d x ^) 2 


1 (G(r))i/j + ( d x r])(d x ^)) 2 

2 1 + (d xV y 


H{rj) + Pext- 


( 1 . 2 ) 


This system is augmented with initial data 


r]\t=o = Vim Mt =0 = An- (1-3) 

We consider the case when 77 and ^ are 2 - 7 r-periodic in the space variable x and 
we set T := M/(27rZ). Recall that the mean value of y is conserved in time and 
can be taken to be 0 without loss of generality. We thus introduce the Sobolev 
spaces Hq(T) of functions with mean value 0. Our main result asserts that, given 
any control domain u> and any arbitrary control time T > 0, the equation (1.2) is 
controllable in time T for small enough data. 
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Theorem 1.1. Let T > 0 and consider a non-empty open subset co C T. There 
exist a large enough and a positive constant M$ small enough such that, for any two 

pairs of functions {rn n ,if in ), (vfinah'f’finai) in H^ +2 (T) x H a ( T) satisfying 

ll ? 7*n||^o-+^ + < Mo, \\lifinal ||^o- + l + || Tpfinal\\]-[v < Mq, 

there exists P ex t in C°([0, T\; H a (T)), supported in [0,T] x co, that is 

supp.P ext (f, •) Co;, Vt € [0, T\, 
such that the Cauchy problem (1.2)-(1.3) has a unique solution 

M) € C 0 ([0,T];F 0 ff+ ^(T) x H"( T)), 
and the solution (q,if) satisfies {g\ t=T ,if\ t=T ) = (r]finahipfinal)- 

Remark 1.2. i) This result holds for any T > 0 and not only for T large enough. 
Compared to the Cauchy problem, for the control problem it is more difficult to 
work on short time intervals than on large time intervals. 

ii) This result holds also in the infinite depth case (it suffices to replace tanh(61£|) 
by 1 in the proof). In finite depth, the non-cavitation assumption g(t, x) > —b holds 
automatically for small enough solutions. 

1.2. Strategy of the proof. We conclude this introduction by explaining the strat¬ 
egy of the proof and the difficulties one has to cope with. 

Remarks about the linearized equation. We use in an essential way the fact that 
the water waves equation is a dispersive equation. This is used to obtain a control 
result which holds on arbitrarily small time intervals. To explain this as well as to 
introduce the control problem, we begin with the analysis of the linearized equation 
around the null solution. Recall that G*(0) is the Fourier multiplier \D X \ tanh(6 |Z) X |). 
Removing quadratic and higher order terms in the equation, System (1.2) becomes 

f d t g = G(0)if, 

{ d t ip + grj - dig = Re¬ 
introduce the Fourier multiplier (of order 3/2) 

L:= {(g-dl)G( 0))U 

The operator G(0) -1 is well-defined on periodic functions with mean value zero. 
Then u = if — iLG{f))^ l g satisfies the dispersive equation 

dtu iLu — Text- 

To our knowledge, the only control result for this linear equation is due to Reid who 
proved in [41] a control result with a distributed control. He proved that one can steer 
any initial data to zero in finite time using a control of the form P ex t(t , x ) = g{x)U[t) 
(g is given and U is unknown). His proof is based on the characterization of Riesz 
basis and a variant of Ingham’s inequality (see the inequality (1.12) stated at the end 
of this introduction). In this paper we are interested in localized control, satisfying 
P ex t(t,x) = 1 u)P ex t.{t,x) where lo C T is a given open subset. However, using the 
same Ingham’s inequality (1-12) and the HUM method, one obtains a variant of 
Reid’s control result where the control is localized. 

Proposition 1.3. For any initial data Ui n in L 2 (T), any T > 0 and any non empty 
open domain co C T, there exists a source term f G C°([0, T]; F 2 (T)) such that the 
unique solution u to 

d t u + iLu = 1 w f ; u\ t=0 = u in , (1.4) 

satisfies u\t=r = 0 (here l w is the indicator function of uj). 
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Remark 1.4. i) (Null-controllability for reversible systems ). By considering the 
backward equation, the same results holds when one exchanges u\t=o and u\t=T- 
Now by using both results, one deduces that for any functions Uj n ,u/maZ in L 2 (T), 
there exists / such that the unique solution to (1.4) satisfies u\t=T = Ufinai■ This is 
why one can assume without loss of generality that the final state Ufi na i is 0. 

ii ) ( Real-valued control ). This result is not satisfactory for the water waves prob¬ 
lem since the control / given by Proposition 1.3 could be complex-valued. To obtain 
a real-valued control requires an extra argument. 

Step 1: Reduction to a dispersive equation. The proof of Theorem 1.1 relies on differ¬ 
ent tools and different previous results. Firstly, Theorem 1.1 is related to the study 
of the Cauchy problem. The literature on the subject goes back to the pioneering 
works of Nalimov [39], Yosihara [44] and Craig [15]. There are many results and we 
quote only some of them starting with the well-posedness of the Cauchy problem 
without smallness assumption, which was first proved by Wu [43] and Beyer-Giinther 
[10] for the case with surface tension. For some recent results about gravity-capillary 
waves, we refer to Iguchi [21], Germain-Masmoudi-Shatah [18], Mesognon [34], 
Ifrim-Tataru [20], Ionescu-Pusateri [22, 23] and Ming-Rousset-Tvzetkov [38]. 

Our study is based on the analysis of the Eulerian formulation of the water waves 
equations by means of microlocal analysis. In this direction it is influenced by 
Lannes [27] as well as [5, 2]. More precisely, we use a paradifferential approach 
in order to paralinearize the water waves equations and then to symmetrize the 
obtained equations. We refer the reader to the appendix for the definition of parad¬ 
ifferential operators T a . 

It is proved in [2] that one can reduce the water waves equations to a single 
dispersive wave equation that is similar to the linearized equation. Namely, it is 
proved in this reference that there are symbols p = p(t,x,tf) and q = q(t,x,£f) with 
p of order 0 in £ and q of order 1/2, such that u = T p ip + iT q r] satisfies an equation 
of the form 

P(u)u = P ext with P(u) := d t + T v ^d x + iL 2 (T^L* •), 

where L 2 = ((g — 3 2 )G(0)) 5 , Ty(u) and T c ( n ) are paraproducts. Here V,c depend 
on the unknown u with F(0) = 0 and c(0) = 1, and hence P(0) = dt + iL is 
the linearized operator around the null solution. We have oversimplified the result 
(neglecting remainder terms and simplifying the dependence of V, c on u ) and we 
refer to Proposition 2.5 for the full statement. 

We complement the analysis of [2] in two directions. Firstly, using elementary 
arguments (Neumann series and the implicit function theorem), we prove that one 
can invert the mapping (?y, ip) i->- u. Secondly, we prove that, up to modifying the 
sub-principal symbols of p and q, one can further require that 

/ lmu(t,x) dt = 0. (1.5) 

J T 

Step 2: Quasi-linear scheme. Since the water waves system (1.2) is quasi-linear, one 
cannot deduce the controllability of the nonlinear equation from the one of P(0). 
Instead of using a fixed point argument, we use a quasi-linear scheme and seek P ext 
as the limit of real-valued functions P n determined by means of approximate control 
problems. To guarantee that P ex t will be real-valued we seek P n as the real part of 
some function. To insure that supp P n C u we seek P n under the form 

Pn — Xu) Rg fn ■ 

Hereafter, we fix ui, a non-empty open subset of T, and a C°° cut-off function Xu, 
supported on uj, such that Xu(x) = 1 for all x in some open interval uq C oj. 
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The approximate control problems are defined by induction as follows: we choose 
fn+i by requiring that the unique solution u n+ \ of the Cauchy problem 

P ('U , n)'U j n-\-\ — Xu R6 / n .|_i, U n -^i\f—Q — Uin 
satisfies u(T) = Ufi na i. Our goal is to prove that 

• this scheme is well-defined (that is one has to prove a controllability result 
for P(u n ))\ 

• the sequences (/„) and (u n ) are bounded in C°([0,T]; H a (T))-, 

• the series XX/n+i ~ fn ) and Y,( u n +1 - u n) converge in C°([0, T]; FT 7- ! (T)). 

It follows that (f n ) and (u n ) are Cauchy sequences in C°([0, T]; H a ~2 (T)) (and in 
fact, by interpolation, in C°([0,T]] H a ' (T)) for any a’ < a). 

To use the quasi-linear scheme, we need to study a sequence of linear approximate 
control problems. The key point is to study the control problem for the linear 
operator P(u ) for some given function u. Our goal is to prove the following result. 

Proposition 1.5. Let T > 0. There exists so such that, if IMIc°([o t]-h s o) small 
enough, depending on T, then the following properties hold. 

i) (Controllability) For all a > s o and all 

Uin, Ufinai € H a ( T) := jic £ H a ( T); Im j w(x)dx = 0 

there exists f satisfying \\f\\co([o,T]-,H^) ^ K ( T )(\\ u in\\ H ^ + \\ u final\\ H <r) such that the 
unique solution u to 

P(u)u = XojRef ; u\ t=0 = u in , 

satisfies u(T) — 

ii ) (Stability) Consider another state v! with H^Hcouq t]-h s o) small enough and 
denote by f the control associated to vf. Then 

I \f - f llc°([0,T];/i' <T "l) — K ( T){\\ u in\\ H °- + \\ u final II H c- ) ||m “ M. |lc°([0,T];// s O) ‘ 

Remark 1.6. i) We oversimplihed the assumptions and refer the reader to Section 9 
for the full statement. 

ii) Notice that the smallness assumption on u involves only some iC s °-norm, while 
the result holds for all initial data in H a with a > so- This is possible because we 
consider a paradifferential equation. This plays a key role in the analysis to overcome 
losses of derivatives with respect to the coefficients. 

Step 3: Reduction to a regularized problem. We next reduce the analysis by proving 
that it is sufficient 

• to consider a classical equation instead of a paradifferential equation; 

• to prove a L 2 -result instead of a Sobolev-result. 

This is obtained by commuting P(u) with some well-chosen elliptic operator s of 
order s with 
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and depending on a small parameter h (the reason to introduce h is explained below). 
In particular s is chosen so that the operator 

P(u) ■= A h , B P(u)Ki 

satisfies 

P(u) = P(u) + R(u) (1.6) 

where R(u) is a remainder term of order 0. For instance, if s = 3m with m G N, set 

Ah,s = I + h s C~3 where £ := L 2 (T C L 2 •). 
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With this choice one has [A^ s ,£] = 0 so (1.6) holds with R(u) = [A/^g, T V(u)]^h?s- 
It follows from symbolic calculus that \\R(u)\\c(l 2 ) ~ ||V||wi,°o uniformly in h. 

Moreover, since V(u) and c(u) are continuous in time with values in H s ° (T) with 
sn large, one can replace paraproducts by usual products, up to remainder terms in 
C°([0,T];£(L 2 )). We have 

P(u) = d t + V(u)d x + iL? (c(u)L 2 ■) + R 2 (u) 

where 

R 2 {u) :=R(u) + (r v(M) -V(u))d x +iL^((T c ^ -c(u))L* •). 

The remainder R 2 (u) belongs to C'°([0,T]; C(L 2 )) uniformly in h. On the other 
hand, 

||[A/ ljS ,Xo;]A-i|| £{L2) = 0(h), (1.7) 

which is the reason to introduce the parameter h. The key point is that one can 
reduce the proof of Proposition 1.5 to the proof of the following result. 

Proposition 1.7. Let T > 0. There exists sq such that, if |Mlc°([o t] h s o) small 
enough, then the following properties hold. 

i) (Controllability) For all Vi n € L 2 ( T) there exists f with |]^llo°([o t]-l 2 ) — 

I\(T) \\vi n \\ L 2 such that the unique solution v to P(u)v = y^Re/, v\ t=0 = Vi n 
is such that v(T) is an imaginary constant: 

3b € R / Vx € T, v(T, x ) = ib. 

ii) (Regularity) Moreover \\f\\ CO{[0 T] . H ^ < K ( T ) \\ v in\\ H l ■ 

Hi) (Stability) Consider another state u' with ||m / |Ic 0 ([o t]-h s o) small enough and 
denote by f the control associated to u'. Then 

IU _ f llc 0 ([0,T];L 2 ) A A (T)\\ v in\\ H 3 ||« — U ||cO([ 0i T];H s O) ’ 

Let us explain how to deduce Proposition 1.5 from the latter proposition. Consider 
Uin ; ufinal in H a ( T) and seek / € C°([0, T]; H a (T)) such that 

p(u)u = Xoj Re/, u( 0 ) = u in => u(T) = Uf inal . 

As explained in Remark 1.4 it is sufficient to consider the case where u/maZ = 0. Now, 
to deduce this result from Proposition 1.7, the main difficulty is that the conjugation 
with A/,. s introduces a nonlocal term: indeed, A ^(Xuf) i s n °t compactly supported 
in general. This is a possible source of difficulty since we seek a localized control term. 
We overcome this problem by considering the control problem for P(u) associated 

3 

to some well-chosen initial data Vi n . Proposition 1.7 asserts that for all Vi n € H*( T) 
there is / € C 0 ([0, T]; i7f (T)) such that 

P(u)vi = Xw Re/, v 1 \t=o = v in => vi(T, x) = ib, b € M. 

Define = v 2 (0) where v 2 is the solution to 

P(u)v 2 = [A/j )S , xJ\ Aj/g Re /, v 2 \ t=T = 0. 

3 

Using (1.7) one can prove that the C(H a)-norm of /C is 0(h) and hence I + JC is 
invertible for h small. So, V{ n can be so chosen that Vi n + JCvi n = A h, s Ui n - Then, 
setting / := Aand u := A^g(ui + v 2 ), one checks that 

P(u)u = Xu Re /, u(0) = Ui n , u(T, x) = ib, b S M. 

It remains to prove that u(T) is not only an imaginary constant, but it is 0. This 
follows from the property (1.5). Indeed, P can be so defined that if P(u)u is a 
real-valued function, then (jj J T Irn u(t, x) dx = 0. Since Irn u(0, x) dx = 0 by 
assumption, one deduces that / T Irn u(T, x) dx = 0 and hence u(T) = 0. 
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Step 4 : Reduction to a constant coefficient equation. The controllability of P(u) 
will be deduced from the classical HUM method. A key step in the HUM method 
consists in proving that some bilinear mapping is coercive. To determine the appro¬ 
priate bilinear mapping, we follow an idea introduced in [1] and conjugate P(u) to 
a constant coefficient operator modulo a remainder term of order 0. 

To do so, we use a change of variables and a pseudo-differential change of un¬ 
knowns to find an operator M(u ) such that 

M(u)P(u)M (u)- 1 = d t + iL + U(u ), 

where \\P(u)\\ c(l 2 ) 1$ IMI.ff s o (and hence 7 Z(u) is a small perturbation of order 0). 

To find M(u) , we begin by considering three changes of variables of the form 

(1 + d x K,(t, x))5 h(t,x + n(t, x)), h(a(t),x), h(t,x — b(t)), 

to replace P(u) with 

Q(u) = 8 t + W8 X +iL + R 3 , (1.8) 

where W = W(t,x ) satisfies f T W(t,x)dx = 0, || W|| C '0([ 0 , T ] ; ipo-<*) < INIc°([o,T] ; .ff s o) 
where d > 0 is a universal constant, and R 3 is of order zero. This is not trivial 
since the equation is nonlocal and also because this exhibits a cancellation of a term 
of order 1/2. Indeed, in general the conjugation of L? (c(u)L 2 ■ ) and a change of 
variables generates also a term of order 3/2 — 1. This term disappears here since we 
consider transformations which preserve the L 2 (dx) scalar product. 

We next seek an operator A such that i [A, \D X \ 2 J +Wd x A is a zero order operator. 
This leads to consider a pseudo-differential operator A = Op(a) for some symbol 
a = a(x,£) in the Hormander class f) with p = namely a = exp(z|£|2/3(i,x)) 
for some function j3 depending on W (see Proposition 5.8 for a complete statement 
that also includes a zero order amplitude). Here we follow [1], To keep the paper 
self-contained (and since some modifications are needed), we recall the strategy of 
the proof in Section 5. 

Concerning the latter transformation, let us compare the equation P(u)u = 0 
with the Benjamin-Ono equation: 

dtw + wd x w + Tid x w = 0, (1.9) 

where P is the Hilbert transform. The control problem for this equation has been 
studied through elaborate techniques (see for instance the recent paper [30]) that 
are specific to this equation and cannot be applied to the water waves equations 1 . 
On the opposite, let us discuss one difference which appears when applying to (1.9) 
the strategy previously described. Given a function W = W(t, x) with zero mean in 
x , let us seek an operator B such that the leading order term in [B,Pd x ] + Wd x B 
vanishes. This requires (see [7]) to introduce a classical pseudo-differential operator 
B = Op(6) with b E S'/o- Then the key difference between the two cases could be 
explained as follows: For r large enough, 

• the mapping W i —> B is Lipschitz from H r into £(L 2 ); 

• the mapping W i —> A is only continuous from H r into C{L 2 ) (indeed, if 
\\W\\ H r = 0(5) then we merely have || A - I \\ c( ^ L2 . H -^ ) = 0(5)). 

This is another reason for which one cannot use a fixed point argument based on a 
contraction estimate to deduce the existence of the control. 


1 This can be seen at the level of the Cauchy problem: for the Euler equation with free surface, 

the well-posedness of the Cauchy problem in the energy space is entirely open. 
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Step 5: Observability. Then, we establish an observability inequality. That is, we 
prove in Proposition 7.1 that there exists e > 0 such that for any initial data vq 
whose mean value (vo) = 2 f T vo(x) dx satisfies 

|Re(u 0 )| > ^|(uo)|-e|KII L a, ( L1 °) 


the solution v of 


dtv + iLv = 0, u(0) = vq 


satisfies 

|Re(Au)(t,x )| 2 dxdt >K |uo (^)| 2 dx. (1.11) 

J T 

To prove this inequality with the real-part in the left-hand side allows to prove the 
existence of a real-valued control function; a similar property is proved for systems 
of wave equations by Burq and Lebeau in [11], 

The observability inequality is deduced using a variant of Ingham’s inequality 
(see Section 6 ). Recall that Ingham’s inequality is an inequality for the L 2 -norm of 
a sum of oscillatory functions which generalizes Parseval’s inequality (it applies to 
pseudo-periodic functions and not only to periodic functions; see for instance [26]). 
For example, one such result asserts that for any T > 0 there exist two positive 
constants C\ = C\{T) and C 2 = C 2 (T ) such that 



Ci ^ \ w n\ 2 



2 

dt < C 2 E I w n 

nEZ 


( 1 . 12 ) 


for all sequences (■ w n ) in £ 2 (C). The fact that this result holds for any T > 0 (and 
not only for T large enough) is a consequence of a general result due to Kahane on 
lacunary series (see [25]). 

Note that, since the original problem is quasi-linear, we are forced to prove an 
Ingham type inequality for sums of oscillatory functions whose phases differ from 
the phase of the linearized equation. For our purposes, we need to consider phases 
that do not depend linearly on t, of the form 

sign(n) [£(n)M + /3(t, x)\n\ 2 ], £(n) := ^(<7 + |n| 2 ) |n| tanh (6 |n|)^ 2 , 

where x plays the role of a parameter. Though it is a sub-principal term, to take 
into account the perturbation f3{t,x)\n\^ requires some care since — 1 is 

not small. In particular we need to prove upper bounds for expressions in which we 
allow some amplitude depending on time (and whose derivatives in time of order k 
can grow as |n| fc / 2 ). 


Step 6: HUM method. Inverting A, we deduce from (1.11) an observability result for 
the adjoint operator Q(u)* ( Q{u) is as given by (1.8)). Then the controllability will 
be deduced from the classical HUM method (we refer to Section 8 for a version that 
makes it possible to consider a real-valued control). The idea is that the observability 
property implies that some bilinear form is coercive and hence the existence of 
the control follows from the Riesz’s theorem and a duality argument. A possible 
difficulty is that the control P ex t is acting only on the equation for i/j. To explain 
this, consider the case where ( Tjfinahipfinal ) = (0,0). Since the HUM method is 
based on orthogonality arguments, the fact that the control is not acting on both 
equations means for our problem that the final state is orthogonal to a co-dimension 
1 space. The fact that this final state can be chosen to be 0 will be obtained by 
choosing this co-dimension 1 space in an appropriate way, introducing an auxiliary 
function M = M(x) which is chosen later on. 
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Consider any real function M = M{x ) with M — 1 small enough, and introduce 

L 2 M := € L 2 (T; C) ; Im J M(x)ip(x)dx = 0 

Notice that is an M-Hilbert space. Also, for any vq € L 2 M , the condition (1.10) 
holds. Then, using a variant of the HUM method in this space, one deduces that 
for all Vi n € L 2 (not necessarily in L 2 M ) there is / € C°([0, T]; L 2 ) such that, if 

Q(u)w = d t w + Wd x w + iLw + R 3 w = Xu Re /, w(0) = w in , 

then 

u;(T, x) = ibM[x ) 

for some constant b £ R. Now 

Q(u) = $(u) _1 P(u)$(u), 

where ^(u) is the composition of the transformations in (1.7). Since ^(u) and 
<h(u) _1 are local operators, one easily deduce a controllability result for P(u) from 
the one proved for Q(u). Now, choosing M = &(u(T, •))(!) where 1 is the constant 
function 1, we deduce from w(T,x ) = ibM(x ) that u(T,x ) is an imaginary constant, 
as asserted in statement i) of Proposition 1.5. Concerning M, notice that M / 1 
because of the factor (1 + d x n{t,x ))a multiplying h(t,x + n(t,x)) in (1.7). 

Step 1: Convergence of the scheme. Let us discuss the proof of the convergence of the 
sequence of approximate controls (f n ) to the desired control P ex t ■ This part requires 
to prove new stability estimates in order to prove that (/ n ) and (u n ) are Cauchy 
sequences. This is where we need statement ii ) in Proposition 1.5, to estimate the 
difference of two controls associated with different coefficients. To prove this stability 
estimate we shall introduce an auxiliary control problem which, loosely speaking, 
interpolates the two control problems. Since the original nonlinear problem is quasi- 
linear, there is a loss of derivative (this reflects the fact that the flow map is expected 
to be merely continuous and not Lipschitz on Sobolev spaces). We overcome this 
loss by proving and using a regularity property of the control, see statement ii) 
in Proposition 1.7. This regularity result is proved by adapting an argument used 
by Dehman-Lebeau [17] and Laurent [29]. We also need to study how the control 
depends on T or on the function M. 

1.3. Outline of the paper. In Section 2 we recall how to use paradifferential 
analysis to symmetrize the water waves equations. As mentioned above, the control 
problem for the water waves equations is studied by means of a nonlinear scheme. 
This requires to solve a linear control problem at each step. We introduce in Section 3 
this linear equation and state the main result we want to prove for it. In Section 4, we 
commute the equations with a well-chosen elliptic operator to obtain a regularized 
problem. Once this step is achieved, we further transform in Section 5 the equations 
by means of a change of variables and by conjugating the equation with some pseudo¬ 
differential operator. Ingham’s type inequalities are proved in Section 6 and then 
used in Section 7 to deduce an observability result which in turn is used in Section 8 
to obtain a controllability result. In Section 8 we also study the way in which the 
control depends on the coefficients, which requires to introduce several auxiliary 
control problems. Eventually, in Sections 9 and 10 we use the previous control 
results for linear equations to deduce our main result Theorem 1.1 by means of a 
quasi-linear scheme. 

To keep the paper self-contained, we add an appendix which contains two sections 
about paradifferential calculus and Sobolev energy estimates for classical or parad¬ 
ifferential evolution equations. The appendix also contains the analysis of various 
changes of variables which are used to conjugate the equations to a simpler form. 
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2. Symmetrization of the water waves equations 


Consider the system 
' d t rj = G(rj)ip, 

< i 

d t ip + 9V + 2 (d x ^) 2 


1 {G{rj)ip + {d x g)(d x ip))‘ 2 

2 1 + ( d x r]) 2 


H{rj) + P ex t(t,x). 


( 2 . 1 ) 


In this section, following [2, 5] we recall how to use paradifferential analysis to rewrite 
the above system as a wave type equation for some new unknown u. This analysis is 
performed in §2.2. In §2.1 and §2.3, we complement the analysis in [2, 5] by proving 
that all the coefficients can be expressed in terms of u only. 

We refer the reader to the appendix for the definitions and the main results of 
paradifferential calculus. 


2.1. Properties of the Dirichlet-Neumann operator. We begin by recalling 
that, if r] in W 1,00 (T) and ip in H 2 (T), then G(jq)ip is well-defined and belongs to 
H~ 2 (T). Moreover, if (r/, ip) belongs to H s ( T) x 77 S (T) for some s > 3/2, then G(rj)ip 
belongs to if s_1 (T) together with the estimate 

< C( IMI/p ) HV’llijs • (2-2) 

Then, it follows from usual nonlinear estimates in Sobolev spaces that the following 
result holds. 

Following [5, 2], the analysis is based on the so-called good unknown of Alinhac 
defined in the next lemma and denoted by ui as the notation for the control domain 
(both notations will not be used simultaneously). 


Lemma 2.1. Let s > 3/2 and ( 77 , -0) in H S (T) x H S (T), the functions 
n , n , G(rj)ip + (d x ri)(d x il>) 

B{ri)ip:= - 2 -, V{rj)ip ■= d x ip - {B(r])ip)d x ri 

1 + {o x i 7) 

u(v)*P :=ip- T B{r]) ^r) 


belong, respectively, to H s X (T), H s 1 (T), H s ( T) and satisfy 

\\B{ri)ip\\ HS -i + WV^ipWns-! + \\u(ij)ip\\ HS < c( IMI^) M\ H s • 


(2.3) 


(2.4) 


Proof. The estimate for B(rj)ip and V{rf)ip follow (2.2), by applying the usual non¬ 
linear estimates in Sobolev spaces, see (A.17) and (A.15). The Sobolev embedding 
then implies that B(r])ip belongs to L°°(T). As a paraproduct with an L°°-function 
acts on any Sobolev space (see (A. 9)), we deduce that 

WTbwvWhs ;$ \\ b (v)'I’\\l°° \\v\\h* <C(H\h^ IIV'II h > WvW H* ■ (2-5) 

This immediately implies the estimate for a j(rj)ip in (2.4). □ 


Consider a Banach space X and an operator A whose operator norm is strictly 
smaller than 1. Then it is well-known that I — A is invertible. Now write oj{rj)ip 
under the form (/ — A)ip with Aip = Tb(tj)^V- By applying the previous argument, 
it follows from (2.5) that we have the following result. 

Lemma 2.2. Let s > 3/2. There exists £0 > 0 such that the following property 
holds. If \\rj\\ H s < £ 0 , then there exists a linear operator 4 '(rj ) such that: 

i ) for any ip in H 2 (T), 

'F( rj)ui(ri)ip = ip ; 

ii) if ui in H s ( T) then T(ry)u; belongs to H s ( T) and 





Notation 2.3. Hereafter, we often simply write B,V,uj instead of B(r))ip, V(q)il>, 
oj(q)if. It follows from the above lemma that, if r/ is small enough in H s ( T), then B 
and V can be express in terms of 77 and u: 

B = B(rj)'S>{rj)u], V = V(rj)'$>{rj)u. 

We also record the following corollary of the analysis in [5, 2]. 

Proposition 2.4. Let s> sq with sq fixed large enough. There exists 6 E (0,1] such 
that 

G{q)iP = G{ 0)w - d x (T v ri) + F{g)fi (2.7) 

where F{rj)if satisfies 

\\F(rj)ilj\\ HS+ i < C (11^711 jjs ) \\v\\h° Ill’ll h* ■ ( 2 - 8 ) 

Proof. We prove that F(rj)if satisfies the following two estimates: 

\\ f { v )^\\ h ^ ^ c { ll^lljy-) Ill’ll h« 7 ( 2 - 9 ) 

\\ f (v)M\hs- 2 — ) \\v\\h> IIV ; IIh s • ( 2 -io) 

The estimate (2.8) then follows by interpolation in Sobolev spaces. 

Let us prove (2.9). In [5, 2] it is proved that, for any N, when s is large enough, 
G(rj)i/J = \D x \u - d x (T v ri) +F(rj)if where \\F(q)fi\\ Hs+N < C( ||r/|| HS ) \\fi\\ Hs . Now 
notice that (2.7) holds with F{rj)fi> = (\D X \ — G(0))oj + F(q)fi. Since G(0) = 

|L) X | tanh (6 |Z7a;|), the difference \D X \ — G( 0) is a smoothing operator. So using 

the estimate (2.4) for u, we find that \\F(r])ijj\\ Hs+ N is bounded by the right-hand 
side of (2.9). Taking N = 1 gives the desired result. 

We now prove (2.10). As for (2.5), using the paraproduct rule (A.9) and (2.4), 
one has 

\\ u -^\\H> + W d x{ T vv)\\ H .-i £ (II^IIloo + ||V'|| LO o) \\v\\ H s < C( \\r]\\ H s) \\rj\\ H , ||^|| H *, 

hence it is sufficient to prove that \\G(r))if — G(0)ip\\ Hs - 2 is bounded by the rhs 
of (2.10). This in turn will be deduced from an estimate of \\ip'{ t)\\ H s - 2 where 
tp(r) = G(Tr])if. Set B r = B(Tq)tl> and V T = V{Trf)ip. It follows from the com¬ 
putation of the shape derivative of the Dirichlet-Neumann operator (see [27]) that 
tp'(r) = —G(Tq)(B T q) — d x (V T r]). Now the estimate (2.4) implies that Wt'( T )\\h s ~ 2 < 
C{ I billies) 11 11 j/s llV’ll/fs- Integrating in r we complete the proof. □ 

2 . 2 . Symmetrization. As already mentioned, the linearized equations are 

( d t q = G(0)i/>, 

\ d t fi + gq - d 2 x q = Pext, 

where G(0) = \D x \tanh(b\D x \). Introducing the Fourier multiplier (of order 3/2) 

L,= ( {g-dl)G(fi)Y , 

with symbol 

KO '■= {(9+ |£| 2 )A (£)) 5 where A(0 := |C| tanh(61^|) ( 2 . 11 ) 

(so that L = £(D X )), and considering u = if — 7LG(0) _ 1 r/, one obtains the equation 

dtu -{- iLu — Pext- 

The following proposition contains a similar diagonalization of System (2.1). 

ii 


Proposition 2.5. Let a, <7o be such that a > cq with ao large enough. Consider a 


solution ( 77 , if) of (2.1) on the time interval [0,T] with 0 < T < + 00 , such that 




Introduce 


c := (1 + (d x rj) 2 ) 4 , 


( 2 . 12 ) 



where £,X are as in ( 2 . 11 ), % € C°° satisfies x{0 = 1 f or |£| > 2/3 and %(£) = 0 
for |£| < 1/2. Then 


u := T p uj — iT q rj 


satisfies 


d t u + T v d x u + iL 2 ( T C L 2 u) + R(r), -ip) = T p P ext , 


(2.13) 


for some remainder R{ 77 ,^) = Ri(j])if + R 2 (rj)r] with 



(2.14) 


for some 6 € (0,1] given by Proposition 2.4. 

Remark 2.6. Compared to a similar result proved in [2], there are two differences. 
We here obtain a super-linear remainder term (see (2.14)), and secondly we prove 
here that q can be so chosen that T q = 8 x Tq for some symbol Q] namely, 



(2.15) 


This will be used to obtain that f T q rj dx = 0. Since it is not a trivial task to obtain 


these additional properties, we shall recall the strategy of the proof from [ 2 ] and 
give a detailed analysis of the required modifications. 

Proof. The first step consists in paralinearizing the equation. We use in particular 
the paralinearization of the Dirichlet-Neumann operator (see (2.7)). Then, by using 
the paralinearization formula for products (replacing products ab by T a b + I/a + 
TZ(a,b)), it follows from direct computations (see [2]) that 



(2.16) 


where a denotes the Taylor coefficient, which is 


a — g + dtB + V d x B, 


and F 1 and F 2 are given by (see (A. 11 ) for the definition of TZ(a, b )) 
F 1 = F(r ,,)Vh 

F 2 = ( T v T dxV - T vdxV )B + ( T vdxB - T v Tg xB )r) 


+ ±k(b, B) - \n(v, V) + T v n(B, d x p) - n( b, vd x rj). 


On the other hand, the paralinearization estimate A. 13 applied with a = a — 1/2 
implies that 
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r := (1 + ( d x r /) 2 ) 2 








where / G L°°(0, T ; H 2a 2 ) is such that 


\f\\ H 2*-% < c (\\v\\ H ,+0 hW^+i, 


(2.17) 


for some non-decreasing function C. Hence, directly from (2.16), we obtain that 
( d t r) + T v d x r] - G( 0)w = f 1 , 

{ d t uj + T v d x u + grj — d x (T r d x r 7 ) = f 2 + P ex t, 

where 

/ 1 := ^ / 2 := F 2 + + T g _ a p. 

Then introduce ( := T q r\ and 9 := T p u. It is found that 

r dtC + T v d x C - t 9 g( 0 )w = f\ 

\ d t 0 + T v d x e + T p (grj - d x (T r d x rj)) = f 2 + T p P ext , 

where 

/ 2 : = 7p/ 2 + Tg tp UJ + [TyC^, T p ]cj. 

Assuming that (7 and p are as in the statement of the proposition, it easily follows 
from ( 2 . 8 ) and the paradifferential rules (A.3), (A.9) and (A. 6 ) (applied with p = 1 
to bound the operator norm of the commutators [Tyd x , T q ] and [Tyd x ,T p ]) that 

||(/\/ 2 )||„. < c(H„ +J ) Ikll" . { \m„- + IML. + j}• 


It remains to compute T q G(0)u; and T p {gr}—d x {T r d x j})). More precisely, it remains 
to establish that 
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\T q G(0)uj - l\t c l\t p u\\ H a < C{ \\ V \\ H(7+h ) | 
\T p (grj - d x {T r d x r ))) - LiT c L 2 T q r]\\ Ra < C( 


' H "+h IMI ^ CT ’ 


(2.18) 


l H a+ 5 


) 


1 H a+ 2 


'jP+S 


(We prove below these estimates with 9 = 1.) Then the estimates (2.14) follow from 
(2.4) which gives a bound for ||u;||# CT in terms of ||V ; ||_h^- 

To prove (2.18), it is convenient to introduce the following notation: Given two 
operators, the notation A ~ B means that, for any p G R there is a constant 
C{h\\ HtT+ i ) such that 

||(A - B)u\\ Hlt < C( \\v\\ H<7+ i ) Ihll^+i \M HI i ■ 

In words, A ~ B means that A equals B modulo a remainder which is of order 0 
and quadratic. 

For instance consider real numbers m, m! with m + m! = 2 and two operators 
A — T a (rn) i) and B — F^( m /) —i) where 




G T 


2 > 


(771—1) 


g r 


m— 1 
1 ’ 


b (m') g b (m'- 1 ) £ pn 


m'—1 


(see Definition A.2) with (see (A.l)) 

M m'( 6 (m')) + M™'- 1 ( 6 ( m '- 1 )) < C( 

M^(a (m) ) + Mf-^a^-^) < C 




)> 


'iP+2 ' 11 /n Jp+2 

By using (A.5) applied with p = 2 and (A. 6 ) applied with p = 1, we obtain that 




~ 0 , 


SO 


AB~T 


a <.™)b( rrl ') + jd£a( rn ')d x b( rn '')+a( rn '>b < - rn '- 1 ')+a( rn - 1 ')b ( rn '')' 
13 


(2.19) 


Using the previous notation, to prove (2.18) we have to prove that 

TqG{ 0) ~ L^T c L^x{D x )T p 

T p (gl - d x (T r d x -)) X (D x ) ~ L?T c L?x(D x )T q . 


( 2 . 20 ) 


Notice that x{^x)v = V and L 2 x(D x )u = L^u for any periodic function u. This is 
why we can introduce the cut-off function x hr the calculations. This cut-off function 
is used to handle symbols which are not smooth at £ = 0 . 

We remark that, by definition of paradifferential operators, we have 

T q G{ 0) = T q x(£), gl - d x (T r d x -) = Tg +r ^ 2 _^ dxr ^y 

Study of the first identity in (2.20). It follows from symbolic calculus (see (A.5)) 
that 

L^T c L^x{p x ) ~ T 7 with 7 = yc£ + — ^/Id x c. ( 2 - 21 ) 

Now we seek q under the form q = gl 1 / 2 ) _|_ ^(- 1 / 2 ) w j iere q( 1 / 2 ) i s Q f order 1/2 in £ 
(more precisely, q E T^ 2 ) and gl -1 / 2 ) is of order —1/2 (in Tj 1 ^ 2 ). Similarly, we seek 
p = p!°) + pl _1 ) with p E T® and pl _1 ) 7 T^ 1 . 

Also, it follows from (2.19) that L^T c L^x(D x )T p ~ TfiT p ~ T pi with 


pi = 7 + xc£p ( + -rX c { 9 ^£)d x p^ 


(the contribution of (dyx)d x p^ is in the remainder term). The first identity in 
( 2 . 20 ) will be satisfied if 


( 1 / 2 ) , = ( 0 )£(0 

Q ' XP A( 0 : 


(—1/2) _ X ^(0 
7 < A(0 


^(d x c)pl 0) + c5 x p (0) 


+ Xcp 


(-i)M 
A(0' 


Study of the second identity in (2.20). As above, it follows from symbolic calculus 
that L 2 T c L 2 x(D x )T q ~ T 7 T 9 ~ T p2 with (see (2.19)) 


-p 2 

1 


p 2 = 7 9 (1/2) + -xc(<%0<W 1/2) + xcf<? ( 1/2) - 


With gl 1 / 2 ) and gl 1 / 2 ) as given above, we compute that 

P2 = x{c^ (1/2) + f (cP^^c + c 2 ^^) + 
Moreover, by definition of £(£) one has 
% 


X c 2 p(-D^) 
* P A(f) r 


. £{t) 2 


m 


= 3^ + ri, 


m 


= + r 2 , ri, r 2 are of order 0 . 


Notice that the contribution of the term rficp^ d x c + c 2 <9 x pl°)) to T p2 and the one 
of r 2 C 2 pG 1 '> can be handled as remainder terms and hence 


L 2 TcL 2 x{Dx)T q 


with 

P 2 = Xi 
On the other hand, 


jc^gl 1 / 2 ! + -w£^cpl°)<9 x c + c 2 d x p(°^ + xc 2 pl ^£ 2 j. 
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By definition of c,£, g^ 1 / 2 ), recall that r = c 2 and £ 2 = (g + £ 2 )A(£) and hence 
p(g + rO - (d x r)(iO) = pc 2 0 + gp ~ p(d x r)(i0 

AO 2 


= pc 


m 


- gpc +gp-p(d x r)(i£). 


Since gB/ 2 ) := xcp(°)^|y, we deduce that 

p(3 + ^ 2 - (<V)(0)x = c^g (1/2) +x{p (_1) c 2 (g + £ 2 ) +gp{l - c 2 ) - zp(3 x r)£j. 

Since 1—c 2 and d x r depend at least linearly on g and since p and are symbols of 

order 0, it follows from the estimate (A.3) for the operator norm of a paradifferential 
operator that 
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\Tp(i-c 2 ) u \\ Hl i < 


l H CT+ 2 


) 




i « 


H a* > 


) IIT7H ..^i ||u| 


■ 


||^p( 1) (d :c r)£' u |l,H> — ^ ( ll r /ll^<r+j ) n-n\ H o+% 

Similarly, assuming that is a symbol of order —1 depending linearly on g (as 

this will be true, see (2.12)), we have 

II< C { hW H °+\) M H «+h IMI h* ■ 

Therefore, 

T p (gI-d x (T r d x -)) X (D x ) 

^ -^c£q , ( 1 / 2 )+xp( —1 )c 2 ^ 2 —(5®r)(i$) * 

Now since r = c 2 , with = c _ 3 ? we have 

-j) ( 0 ) (^r)(iO = +^(cp^d x c + c 2 d x p^, 

as can be verified by a direct calculation, so the second identity in (2.20) holds. 

It remains to compute q. We have 


q = x \cp 


M 


m . iw 
m i m 


^( d x c)p (0) + cd x p( 0) 


+ CP H)11 

P A(0/' 


Observe that 


L ^(<9 x c)p (0) + c^j) (0) 


= -c sd x c. 
o 


We now seek p( such that 


p m a i m dl 

for some constant a to be determined. We thus set 

„<-!>:= a hhM|> c - i BxC . 

I 1(0 

Then (replacing y 2 by y, to the price of adding a smoothing operator in the remain¬ 
der), we have 


, 2 m , 1 d^(0 

q -- =x ^W) + ~(w 


a + ^ ) c 3 9 x c 


Since 


2 V + T (0 

with r(0 is a smooth symbol of order 1/2, we have 

2 _ t 4 _ dJL 

3 x ®rni-r* xi0 m +r 
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where r' is of order —3/2. Then, choosing a such that a + 1/6 = 4/9, we find that 


where f is such that 


q = c s 


m 

m 


+ (d x c§) 


*eA(o 


+ r 


\\Tfu \\ Hll+ 3 < C( \\v\\ H „+i ) \\v\\ H *+l IMI/r#* ■ 

In particular, the contribution of r can be handled as a remainder term and the 
same results hold when q is replaced by the same expression without f, thereby 
obtaining (2.12). This completes the proof of (2.20) and hence the proof of the 
proposition. □ 


2.3. Invertibility of the change of unknowns. We thus obtained an equation 
of the form 

d t u + T v d x u + iL^ (T c L^u) + R(r), 7 / b) = T p P ext , 

where the coefficients V and c depend on the original unknowns ( 77 , 7 /). We conclude 
this section by proving that V and c can be expressed in terms of u only. We have 
already seen in Lemma 2.2 that these coefficients can be expressed in terms of 77 
and u). So it remains only to express (rj,ui) in terms of u. 

In this paragraph, the time is seen as a parameter and we skip it. 

Notation 2.7. Introduce the space H a ( T;C) of complex-valued functions u satis¬ 
fying 

/ lmu(x)dx = 0. 

J T 

Recall (see (A. 3)) that a paradifferential operator with symbol in T™ is bounded 
from any Sobolev space T) to Recall also that ui € H a ( T) whenever 

( 77 , 7/7) g Hq + 2 ( T) x H a (T). Since, as already mentioned, T q r] = T qx r] where x is 
as defined after (2.12) and since qx € Tq , we deduce that u belongs to H a ( T). 
Moreover, it follows from (2.15) that T p u — iT q r] belongs to H a ( T;C). 

We now introduce the mapping U : iLp +5 (T) x H a (T) —> H a ( T;C) so that 

1/(77, VO := T p u - iT q r]. 

The following result shows that this nonlinear mapping can be inverted. 

Lemma 2.8. Let ctq > 5/2. There exists eo > 0 and K such that the following 
properties holds. If\\ri\\ H a 0 < £0, then there exists 

Y: H ao ( T;C) -»• H° 0+ ^ ( T) x H ao (T), 
such that Y(u ) = ( 77 , 7 /) with u = U(r],ijj). Moreover, for any o > 5/2, 

^ ll’ u lli?°‘ 5 \H\\m < 2 ||u||^ ■ (2.22) 

Proof. Set u = 17 ( 77 , 7 /) := T p co — iT q ip Then T q r] = — Iiutt and T p co = Re u, where 
q and p depend on 77 . The only difficulty is to express 77 in terms of IniTX. Once this 
will be granted, to invert the equation T p uj = Re u, we use the fact that T p is a small 
bounded perturbation of the identity so that T p is invertible, indeed (recalling that 
M'pf a) is defined by (A.l)) 

\\T, - 7|| £(H .) < M°(p - 1) < C (||.;II H . +1 ) I|I)II H . + 1 . 

Now to solve the equation T q r] = — Im u, we use the Banach fixed point theo¬ 
rem. Denote by Q the Fourier multiplier with symbol Q(£) := x(£X(£)M(£) = 
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X&Va + e/VW)- The reason to introduce this symbol is that, with q as given 
by (2.12) one has 

M 0 1/2 (</(*,£) - Q{0) < C(M h „ + i)M h „ + i , (2.23) 

which is obtained by considering separately the principal and sub-principal terms in 
the definition of q. Then seek q in Hq +2 (T) such that &(q) = q with 

®(jl) '■= ~{a ~ dl)~^G(0)^({T q -Q)q + lmu). 

S+ 1 

It is easily verified that if $( 77 ) = q then T q q = — Imu and also that maps H 0 2 (T) 
into itself. To see that is a contraction, we use (2.23) to obtain 

l|T(r?i) - $( 772 ) 11^+1 < || {T qi -Q)(m -772)11#* + \\( T qi - T q 2 )m\\ H ° 

< M Q 1/2 (g 1 - Q) \\qi - q 2 \\ H<7+ 1 + M Q 1/2 (g 1 - q 2 ) IMI^+i 
<C{M)M\\q 1 -q 2 \\ Ha+h 

where M := ||r/i | CT .i+||r/ 2 [| CT+ i • If M is small enough, then <J> is a contraction. □ 

H "2 H 2 

3. The linear equation 

As mentioned in the introduction, we shall study the control problem for the 
water waves equations by means of a nonlinear scheme. This requires to solve a 
linear control problem at each step. We introduce in this section the linear equation 
we are going to study until section 10, emphasize one key property of this equation 
and state the main result we want to prove. 

We have seen in the previous section that one can express V = V(q)ip in terms of 
u only. To simplify notations, we write V = V(u), and similarly we write c = c(it). 
Also, one can write the remainder R(q, ip) under the form R{u)u where, for any u, 
the mapping u eA R(u)u is linear. 

We have proved that, for a large enough and a solution (q, ip) of (2.1) on the time 
interval [0, T\, satisfying 

( t ,i)eC°([0,r];^(T)xr(T)), 

the new unknown u satisfies u G C 0 ([0, T]; H a {T] C)) (where ff CT (T;C) is defined in 
Notation 2.7) and 

d t u + T v ^d x u + iL I ( T^L^u ) + R(u)u = T p ^P ext . 

We now fix u £ C'°([0, T];H a ( T; C)) and then set 

V = V(u), c = c(u), R = R(u ), p = p(u) (3.1) 

and consider the linear operator 

P = d t + Tyd x + id (: r c d -)+r. 

Except for the second condition in Assumption 3.1 below, we shall not use the way 
in which the coefficients depend on u and hence we shall state all the assumptions 
on V,c,p,R forgetting their dependence on u through (3.1). 

Assumption 3.1. i) Consider two real-valued, functions V, c in C'°([0, T\; H s ° (T)) 
for some sq large enough, with c bounded from below by 1/2. The symbol p is 

given by p := c“s + c~3d x c with x as in (2.12). R is always assumed 

3 

that the W^'°°-norm of c— 1 is small enough. 
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ii) If Pu is a real-valued function then 

d r 

— / x) dx = 0. 

dtjj V ' 


Fix an open domain u C T and denote by Xu a C°° cut-off function such that 
Xu{x) = 1 for x € oj. We want to study the following control problem: given an 
initial data V{ n find / such that the unique solution to 


Pv = T p Xoj Re /, v\ t=0 = Vi 


(3.2) 


satisfies v\ t=T = 0. The fact that the Cauchy problem for (3.2) admits a unique 
solution is proved in the appendix, see Proposition B.l. 

Our main goal until Section 10 will be to prove the following control result. 

Proposition 3.2. There exists sq large enough such that, for all T £ (0, 1] and all 
s > so, if Assumption 3.1 holds then there exist two positive constants 5 = 6(T,s ) 
and K = K(T, s) such that, if 


(1 < k < 3), 


(3.3) 


ll^llc°([0,T];H s O) + ll c 1 |Ig°([0,T];H s O) ^ d, 
ll^t ^llc 0 ([0,T];,fP) + 11^ C llc°([0,T];ff 1 ) — ^ 
ll-^llc°([0,T];Z:(i? s )) — d, 

then for any initial data Vi n £ H s { T;C) there exists f £ (7 o ([0, T]; H S (T)) such that: 
(1) the unique solution v to Pv = T p Xw Re /, v\t=o = v in satisfies v(T) = 0; 


(%) ll/llc°([0,T];if s ) — K \\ V ' 


in\\H s * 


Remark 3.3. Notice that the smallness assumption on V and c involves only some 
F s °-norm, while the result holds for initial data in H s with s > sq. We shall use this 
property with so = s — 2 in the analysis of the quasi-linear scheme. This is possible 
only because we consider a paradifferential equation. 

We conclude this section by proving that the second condition in Assumption 3.1 
holds when V,c,p,R are given by (3.1). 

Lemma 3.4. Consider u € C'°([0, T]; H S °(T] C)) with sq large enough and assume 
that V,c,p,R are given by (3.1). If Pu is a real-valued function, then 

d r 

— / Imit(t, x) dx = 0. 
dtjj V 

Proof. Set C = Im u. It follows from (2.17) that satisfies 

dtC + T v d x ( - T q G(0)u) = f\ 
f 1 = TqiPiv)^ - T d x vn) + T dtq ri + [: T v d x ,Tq]r) , 
where Fig)^ is given by (2.7). One can write this equation under the form 

dtC + T q {d x (T v r,)) - T q G(t))io = T q F{r,)^ + T dtqV . (3.4) 

Notice that T q G(t))uj and T q F(7j)ip are well-defined since G(0)w(0) = 0 = F(?y)V’(0) 
(this follows from the definition (2.7) and the fact that the mean values of Gfrfiif, 
G(0)uj and d x (Typ) are all 0). Using (2.15), one obtains that J T T q vdx = 0 = 
f T Tg tq v dx for any function v. So integrating (3.4) we obtain the desired result. □ 
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4. Reduction to a regularized equation 


In this section, we reduce the proof of Proposition 3.2 to that of a simpler result. 
We shall prove that: 

• it is enough to consider a classical equation instead of a paradifferential 
equation (this observation will be used below to simplify the computation of 
a change of variable); 

• it is enough to prove an L 2 -result instead of a result in higher order Sobolev 
spaces (this plays a crucial role). 

As explained in the introduction, the idea is to commute the equation with an 
elliptic semi-classical operator Ah lS of order s. To choose this elliptic operator, 
the key point is to prove that Ah, s can be so chosen that it satisfies the following 
commutator estimates: 


|| [A/i )S , P]\,sWc(l 2 ) ~ ^(l)’ 

|| [Afc,sj Xw]A h)S ||££ i2 ) = 0(h ), 

which is the reason to introduce the small parameter h. Some care is required to do 
so, and we introduce 

Ah,s = I + h s T*sL 2 f (4.1) 


Lemma 4.1. i) Assume that the L^° x -norm of c — 1 is small enough. Then A^ s is 
invertible from H s to L 2 and its inverse is denoted by A ^). 

ii) Moreover, for any s' < s, h s A^ 2 is uniformly bounded from L 2 to H s ■ 
is K > 0 such that for any h £ (0,1] and any u in L 2 (T), 


there 




< K llul 


H s ‘ 


L 2 ' 


(4.2) 


Proof. Set r = 2s/3. Statement i) is obtained writing A^ g as (I + B)(I + h s L r ) 
where B is a bounded operator from L 2 into itself. To do so, write 


A h)S = / + h s T c rL r = I + h s L r + h s T c r_ 1 L r , 

to obtain the desired result with B := h s T c r_iL r (I + h s L r )~ 1 . We now claim that B 
is a bounded operator on L 2 , with operator norm 0(||c — l|| i0 o). This follows easily 
from (A.9) (which implies that T c r_i is of order 0 with operator norm 0(\\c — l|| icx) )) 
and, on the other hand, from the fact that the operator h s L r (l + h s L r ) 1 is bounded 
from L 2 into itself uniformly in h (as can be verified using the Fourier transform). 

Now for ||c — 1|| LOO small enough, one has ||-B||£(£ 2 ) < 1/2 and one can invert 
(I + B) to obtain 

+ + (4.3) 

and statement ii) follows from the fact that h s (I + h s L r )~ 1 is uniformly bounded 
in C(L 2 ; H s ') for s' < s. □ 


The key property is that one has good estimates for the commutators of A h s and 
the various operators appearing in the equation. 

3 

Lemma 4.2. Assume that the W^’°°-norm of c— 1 is small enough. Then there is 
K > 0 such that for any h £ (0,1] and any u in L 2 (T), there holds 


[Afe )S ,A ft u < A ||R||jyi,oo IMIi 2 , 

(4.4) 

[^h,s,Xui]K}s u \\ L 2 < Kh\\xu\\ H ^ IMIl 2 ) 

(4.5) 

[A h!S ,L 2 (T c L 2 -))A~^ s u\\ L2 < Kh\\u \\ L 2 . 

(4.6) 
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Proof. Write 


[A h , s ,T v d x \A^ s = [T {ce)2s/3 ,T v d x ]h s A^ s , 


C(H S ,L 2 ) 


1^11 


C(L 2 ,H S ) W U Wl 2 ■ 


to obtain 

|| [h-h,s,T v d x ]A^ s u\\ L2 < K [T^2s/3 , T v d x 

It follows from (4.2) that ||/i s A)/ sII/^ls H s ) is uniformly bounded in h. On the other 
hand, the commutator estimate (A. 8 ) implies that 

II p n {ce) 2s / 3 i'^vd x ]\\ c ^ HSL 2^ < A || 1 ||jyi,oo, 

where K depends on ||c|| 3 ^ (which by assumption can be bounded by 2). 

To estimate the second commutator, we begin by establishing that 

|| [A/i )S , T Xu \< Kh\ 

To see this write 

= h[T (ct)a ,„,T x Jfc-'A^. 

Then we notice that, as above, 


H a 


\u\ 


L 2 ■ 


(4.7) 


[T ( c0*/3 > T Xu \ 


C{Hs-^,L 2 ) 
-1 A-l 


T K IIxojIIw 1 .' 


and we use that, thanks to (4.2), h s 1 A fl L g is uniformly bounded from L 2 to H s 1 . 

Now it remains to estimate [A^ s , (xu — T Xui )\. It follows from Proposition A .8 
(applied with (r,p, 7 ) = (s + 1 , 0 , s)) that 


^ -^(rf) 2s /3(Xu T Xui )A h ^ u 


L 2 


< h s 


( Xu> T XLL] )A h ' S u 


H s 


h s ||x, 


w||h s + 1 IIA/ i , !S 'w|| L 2 ^ h s Hx^ll^s+i \\u\\ L 2 , 


and similarly 


(x« - T xJ hST (c£) 2 ^\s U L 2 ^ Kh * IIXw|Ih«+i IML 2 • 


By combining these two estimates, we find that 


[A/i,s, (Xu -^xw)]A 


h,s l 


L 2 


< Kh s Hx^ll^s+i 11u11 j^2 ■ 


(4.8) 


By combining (4.7) and (4.8) we deduce (4.5). 

We now prove the last property (4.6). Write that L 2 (T C L 2 -) = T^ + Tp + i? where 
R is of order 0 and p = i~ 1 y/£(d^y/I)(d x c). Since A/j iS = I + h s T^ 2a / 3 , by definition, 
[A/j )S , Lz (T c L2 -)]A^g can be written as the sum (I) + (II) + (III) with 

(I) := [T {ce) 2s /3 , T a ] h s A-], (II) := [T {c£)2s/3 , T p ] h s A^, 

(///) := [T {c£)2s/3 ,R j^A-i. 

Since /i s A^|, belongs to C(L 2 ,H S ) uniformly in h, we need only estimate 

I [ T {cl) 2 °p > II c(H s ,L 2 y II [^(c^) 2s /3) Ap] 11 £(# 8 ^ 2 )- 

The second term is estimated by means of (A.5) applied with p = 1/2. To estimate 
the first term we notice that the Poisson bracket of the symbols vanishes: 

{(ct) 2s /\d} = l -((dt(cl) 2s ^)d x (cl) - (d x (d) 2s ^(cl)) = 0. 

Since ||c|| 3 :00 < 2 by assumption, it follows from (A.5) applied with p = 3/2 that 

II [^~(rf) 2s / 3 5 ||/;(H S ,L2) ^ 1- 

This completes the proof. □ 
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Next we conjugate P with s . Introduce 

P h := A Ks PA^ s . 

Then 

Ph = dt + T v d x + iL? (T C L 2 •) + Ri where 

R\ ■=Ah yS RA h s + [A/j )S , dt] A ftjS + [A^, T v d x ] A h<s + i [A h:S , L? (T c Li ■)] A h fj . 

3 

Lemma 4.3. Assume that the W T°° -norm of c— 1 zs small enough. There holds 
\\R h iu \\ L 2 < K( \\V\\ wllOB + \\dtc\\ Lao + h~ s \\R\\ cm ) ||«|| L 2 , (4.9) 

for some constant K independent of h. 

Remark 4.4. The constant h~ s is harmless since at the end of this section, h will 
be fixed depending only on T. 

Proof. We have 

11 Rh,sRRfi tS 11 £(£ 2 ) — 11 II II -^11 II 11 zz(z. 2 ; j/ s > — R^ 1 

since ||A/i iS || £(H8;L2) < 1 and ||A^|| £(i 2 ; ^ s) < h~ s . 

On the other hand, [A/j )S , (T C L ^■)] Af l s and [A/j )S , Tyd x ] are estimated by 
means of Lemma 4.2 and [A^ iS , dt\ Af l s is estimated by similar arguments. □ 

We further transform the equation by replacing Tyd x and Lz ( T C L 2 •) by Vd x and 
L 2 (cLs-) modulo remainder terms. Namely, write Ph as 

P h :=d t + Vd x + iL^(cL5 ■)+R% (4.10) 

where c stands for the multiplication operator by c and 

R. 2 'u = RiU + Tyd x u — Vd x u + i^L^TcL^u — L 2 ( cL^u )). 

3 

Lemma 4.5. Let sq > 2 and assume that the W 2 ’ 00 -norm of c— 1 is small enough. 
There holds 

\\R h 2 u\\ L2 < K{ ||I /||^ 0 + ||c - 1|U + \\d t c\\ m + h~ s p|| £( ^ } ) \\u\\ L2 , (4.11) 

for some constant K independent of h. 

Proof. We have already estimated R\, and the right-hand side of (4.9) is less than 
the one of (4.11) provided that so > 3/2. To estimate Tyd x u — Vd x ii, we apply 
Proposition A .8 with (r,/ 1 , 7 ) = (so,—1,0) (and so > 3/2) to obtain 

\\Tyd x u — Vd x u \\ L 2 < 11P”11 fls 0 \\d x u\\ H -i < IIVU^sq ||m ||^2 • 

The estimate for C — L 2 (cLz • ) = L 2 ((T C — cL)Lz ■ ) follows in the same way, 
assuming that sq > 2 . □ 

We are now ready to give the main reduction. Our goal in this section is to prove 
that one can deduce Proposition 3.2 from the following proposition. 

Proposition 4.6. Consider an operator of the form 

P :=d t + Vd x + iL?(cL^ •) +R 2 . (4.12) 

Let T E (0,1] and consider an open subset co C T. There exist an integer sq large 
enough and two positive constants 5 = S(T) and K = K[T ) such that, if 

ll^llc 0 ([0,T];.ff s O) + \\ c ~ l|lc' 0 ([0,r];/f s 0) ^ ^ 

ll^t ^llc°([0,T];iT) A ll^t C llc°([0,T];_ff 1 ) — ^ (1 < k < 3), (4.13) 

ll-^2|lc-0([0,T];/:(L 2 )) — 3’ 
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then for any initial data Vi n E L 2 (T) there exists f E C°([0, T]; L 2 (T)) such that: 

(1) the unique solution v to Pv = x^Re/, U| t=0 = V{ n is such that v{T ) is an 
imaginary constant: 

3b E M/ \/x E T, u(T, x ) = ib. 

(%) \\f\\c°([0,T];L 2 ) — K II || £2 • 

Remark 4.7. Notice that the final state v(T ) is not 0 but an imaginary constant. 

This result will be proved later. In the end of this section, we assume that 
Proposition 4.6 is true and prove Proposition 3.2. 

Proof of Proposition 3.2 given Proposition f.6. Proposition 4.6 holds for any P of 
the form (4.12). In particular, in view of (4.10), it holds for P replaced by Ph '■= 
Let us mention that h will be fixed at the end of the proof by asking 
that K'(T)h < 1/4 where K'(T) depends only on T. 

The idea is to apply control property for Ph associated with an unknown initial 
data to be determined. 

We shall prove that Proposition 4.6 implies that Proposition 3.2 holds with the 
conclusion 1 replaced by v(T) E iM. Then one deduces that v(T) = 0 by using 
condition ii ) in Assumption 3.1 and the fact that f T Vi n (x)dx = 0. 

Assume that 5 < h s 5 where 6 appears in the statement of Proposition 3.2 and <5 
is as given by Proposition 4.6. Then one has h~ s 5 < 5. Therefore, if the smallness 
condition (3.3) holds, then Lemma 4.3 implies that ||A ?2 Hc°([o t]-C(L 2 )) * s sma ^) an d 
hence the smallness assumption (4.13) holds. This explains why one may apply the 
conclusion of Proposition 4.6 under the assumption of Proposition 3.2. 

The assumption that Proposition 4.6 holds implies that for any y E P 2 (T) there 
is / E C'°([0, T]; L 2 (T)) satisfying 

I|/IU[0,T1;L*) £ «TT) IMb . (4-14) 

and such that the unique solution u\ to 

PhU\ = Xu Re /, ui |t =0 = y , 

is such that ui(T, x) = ib for some b E M and all x E T. 

Now introduce U 2 which is the unique solution of the Cauchy problem (with data 
at time T), 

P h u 2 = {A hjS T p Xu^-hl - Xu) Re/, u 2 (T) = 0. 

Again, the fact that the above Cauchy problem has a unique solution follows from 
Proposition B.l. One can then define the linear operator 1C by 

JCy = u 2 { 0). (4.15) 

The reason to introduce u 2 and the operator K. is that the function u defined by 
u := u\ + u 2 satisfies 

p h u = A/j )S (TpXu^hl Re 7 ) > U ( T ) = ib i u \t= o = y + Ky- 

Now, assume that / + TC is invertible with (/ + /C) -1 E C(L 2 ). Then y can be so 
chosen that y + ICy = A h, s Vi n - Using that A h )S b = b and hence A//J.6 = b for any 

constant 6, it follows that, with / := A/^/ and v := A 

Pv = T p Xu Re /, v(T) = ib , u(0) = v m , 

where P is the original operator, so that P^ = A^PA^. Moreover, it follows from 

the conclusion 2 of Proposition 4.6 that ||/|| c oq 0 t\l 2 ) — R 7I^Hl 2 R° m which we 
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deduce that ||/|| c o(r 0 t]-h s ) — AT(/i) 112/11 if»- The fact that the last constant depends on 
h is not a problem since h is fixed, depending on T. Now to see that Proposition 3.2 
holds, it remains to check that v(T) = 0. As already mentioned, the fact that 
v(T) = 0 follows from the fact that v(T) = ib as well as condition ii ) in Assumption 
3.1 and the fact that f T Vi n (x) dx = 0. 

Thus it remains only to prove that I + JC is a bijection from L 2 into itself. To see 
this, it is sufficient to prove that /C is a bounded operator whose operator norm in 
C(L 2 ) is strictly smaller than 1. In this direction, we first use the energy estimate 
(B.3) for the operator P : 

IMt )|| L 2 < e CT (jMT)|| i2 + ||Pu 2 || L 2 d/\ , 

for some constant C depending only on 

M so := sup { 11^(011^0 + \\c(t') - i11^ SQ + \\R 2 {t')\\ 2 }. 

t'e[o,T] v ' 

Since u 2 (T) = 0, this implies that, 

f T _ 

ll«2 Ilc°([0, T];L 2 ) < ^ J Q \\{ A h,s T pXu:^ s - Xu) f(t')\\ L2 dt' ■ 

To estimate the term (A/j^T/y^A^* — Xu)f we write it as 

[A/, )S , x»\K)J + a ks(t p - i)xuKlJ- 

It follows from (4.5) that 

||[A/ l ,s,Xw]A^ !S /|| i 2 < A/l HXwllifs+l ||/||_J,2- 

-—- 2s 

It remains to estimate Ah tS {T p —I)XuA^ s f. To do so, we write A/j iS = I+h s T c ( 2 s)/ 3 L^ 
to split this term as 

(T p - I) X uA^J + T c( 2 s) / 3 L% (T p - I)xUh s ^ s )f. 

For the first term we have (using (A.3) and (4.2) with s' = 0) 

lire, - i)x^j\\ L , < m§( P - 1 ) hxji^ ii7ii„. 

For the second term write (using (A. 9), (A. 16) and (4.2) with s' = s ) 

||T c(28)/3 Lf (Tp - I)xu(h s A^ s )f\\ L 2 < || (Tp - I)xUh s A^ s )f\\ Hs 

~ \\Tp-i\\c(H s ) IIXwlljjs II^Aft^/Hjys 

< Mq(p — 1) ||x«|| H » ll/L- 

It is found that 

|| A-h,s(T p — I)XujA h s f\\ L 2 1 $ (II c ~~ 1 || l°° A || 0 *c || L oo ) Hxwllijs ||/||x,2 ~ ^||/||z,2- 
This yields 

11^2 llc°([0,T];L 2 ) ~ (h + 5)e CT f ||/|| L 2 dt'. 

J 0 

In view of (4.14), we conclude that 

II U2 IIc°([o,t];L 2 ) — (T)(h + 5) ||y || L 2 , 

for some constant K'(T). Then chose h,6 such that I\'(T)h < 1/4, K'(T)5 < 1/4. 
We conclude that 

VtE[0,T], \\u 2 {t )\\ L 2 < ^ ||y || L 2 ■ 
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(4.16) 


By applying this inequality with t = 0, one obtains ||/Cy|| i2 < ^ ||y|| i2 which proves 
that I + K, is invertible in £(L 2 ). This completes the proof of Proposition 3.2. □ 

5. Further reductions 

Recall that until now we have reduced the study of the control problem in Sobolev 
spaces for P := dt + Tyd x + iL^T c L^ + R to the one of the control problem in L 2 
for P = dt + Vd x + iL i (cL2-) + R 2 . 

5.1. Change of variables. The goal of this subsection is to reduce the analysis to 
an equation where La (cLa •) is replaced with an operator with constant coefficient. 
To do so, we use three change of variables, which preserves the L 2 (cLx) scalar product. 
This allows us to conjugate P to an operator of the form 

d t + W3 X PiL + R 

where R is of order 0 and furthermore W = W(t, x) satisfies J T W(t , x) dx = 0. 

Proposition 5.1. There exist universal constants 5 q E (0,1), r > 2, C > 0 such 
that the following properties hold. Assume that c, V, R 2 satisfy 

\\ c - l|lc°([0,T];L°°) < $0, A/o < 1, (5.1) 

where 

A/() := ||C— 11|C°([0,T] ,H r ) + ll^llc°([0,T];.ff 1 ) + ll^ c llc , °([0,T];H' 1 ) + 11-^2 llc°([0,T];£(L 2 )) • 
Then there exist a constant T\ > 0 and a bounded, invertible linear map 
T: C°([0,T];L 2 (T)) -> C°([0, Ti]; L 2 (T)) 
with bounded inverse d> _1 such that 

Pu = 

where m = m(t ) is a function of time only, defined for t € [0, T], and 

P 3 = dt + W8 X + iL + R 3 . 

The function W = W(t,x) is defined for t € [0, Ti], it satisfies f T W(t,x) dx = 0, 
and 

ll^llc°([0,Ti];H 2 ) < C'dKc - l^)llc°([0,T];H 2 ) + l|3t c llc°([0,T];.ff 1 )) • ( 5 - 2 ) 

The operator R 3 maps C°([0, Ti]; L 2 (T)) into itself, with 

ll^3||c 0 ([0,Ti];£(L 2 )) < CAf 0 . (5.3) 

The constant T\ and the function m satisfy 

+ ll m— l|lc°([0,T]) < C\\c— l||c°([0,T];L°°)- 
The map is the composition of three local transformations where 

(’Pi h)(t,x) := (1 + d x j3i(t,x))%h(t,x + fafax)), ^ 

(^ 1/i )(A®) :=h(ij;- x (t),x), ((p~ 1 h)(t,x) := h(t,x - p{t)). 

Proof. This proposition is proved in Appendix §C. □ 

Remark 5.2. i) The proof is based on computations similar to the ones used in [1], 
However, the analysis in [1] used some special properties of the Hilbert transform 
which cannot be applied in the present setting. Instead, we shall rely on Egorov 
theorem. Moreover, we need in this paper to introduce a change of variables which 
preserves the skew-symmetric structure of the operator iLi ( cL 2 •)• This allows us to 
prove that some operator of order 1/2 vanishes, which plays an essential role below. 
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This in turn forces us to revisit the analysis of changes of variables, which explains 
why the proof is done in details in §C. 

ii ) In sharp contrast with other transformations that will be performed below, 
notice that a change of variable is a local transformation, thereby transforming a 
localized control in another localized control (this is used below to prove Lemma 9.2). 

In addition to Proposition 5.1, higher regularity and stability estimates are given 
in Proposition C.2. 


5.2. Conjugation. To study the control problem for the new equation 

d t + W8 X + iL + R 3 


we will use the HUM method. A key point is then to prove an observability inequality 
for solutions of the dual equation, which reads 

(~d t - d x (W-) — iL + R* 3 )w = 0. 

This equation can be written under the form Vw = 0 with 

Vw := dtw + Wd x w + iLw + Raw, 


where 

R a w := -R%w + (d x W)w. (5.5) 

The observability inequality will be proved later. As a preparation, in this section, 
we prove that V is conjugated to a simpler operator where dtw + W'd x w is replaced 
by dt . To do so, we use the analysis in [1], For the sake of completeness, we recall 
the strategy and the main steps of the proof. 

We often use below the following notation: given a function / with zero mean, 
dfi 1 / is the zero-mean primitive of /, defined by 

P-'f V 4''A /(,:) 

i¥= o lJ o 

We seek an operator A such that 

(pt + Wd x + iL + R,fi)A = A(dt + iL + Rp 
where R$ is a remainder term of order 0. By definition 

R 5 := A-'^duA] +RaA + + z [L, A] ). (5.6) 

Seeking A as a pseudo-differential operator, and trying to cancel the leading order 
terms (that is Wd x A + i [L , A]), it is natural to introduce A as follows. Let 

< f>(t , x , k) := kx + /3(t, x)|A:| 5 , 

for some function (3 to be determined. Consider also an amplitude q(t,x,k) to be 
determined. Then define the operator A(t) by setting 

Au(t, x) = u k q(t , x, , (5.7) 


for periodic functions u, where Wfc(t) are the Fourier coefficients of u, so that u(t, x) = 

Eke z Mt)e ikx . 

Below t is seen as a parameter and we omit it in most expressions. 


Assumption 5.3. Set 

N '■= ll^llc°([0,T];// s O) + ll c— l|lc°([0,T];H s O) + IM C°([0,T];// 1 ) + ll-^ 2 Hc°([0,T];£(L 2 )) > 

where so is some fixed large enough integer. In this section, we always assume that 
Af is small enough without recalling this assumption in all the statements. 
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Hereafter, so always refers to an index large enough whose value may vary from 
one statement to another. 


Lemma 5.4 (from [1]). There exists a universal constant 6 > 0 with the following 
properties. 

i ) Consider the case when the amplitude q is of order zero in k and is a pertur¬ 
bation of 1, 

q(x, k) = 1 + b(x, k ). 

Denote |6| s := sup fceZ ||6(-, &)lltf s (T)- V 

m\m + \b\ 3 <5, 

then A and A* are invertible from L 2 (T ) onto itself, with 

IM'Ml 2 + 11-A 1 u\\l 2 + ||H* u \\ L 2 + ||(v4*) 1 u \\ L 2 < C \\u \\ L 2 , 
where C > 0 is a universal constant. 

ii ) Consider the case when the amplitude q is small. That is, assume that 

WPWh 3 + Ms < b , 

then 

|| Au||£2 < C5 \\u\\ L 2 , 
where C > 0 is a universal constant. 

Proposition 5.5 (from [1]). Assume that ||/3||^i,oo <1/4 and ||/3||#2 < 1/2. Let 
r, m,so€M, m > 0, so > 1/2, Me N, M > 2(m + r + 1) + so- 

Then 

\D x \ r Au = J2 °P (^T (W) (q(x,k)e l \ k \^^u + R M u, 

where, for every s> sq, the remainder satisfies 

|| Rm \D x \ m u\\ H s < C'(s)|/C 2 (m+r+s 0 +l) \\ u \\h s + K-s+M+m+2 IM|_ffso j, (5.8) 

where K,f := |q - 1| M + \q\i\\/3\\ H ^+i and \q\^ := sup t sup fceZ || q{t, ■, k)\\ H r. 

We now deduce the following result (which is a variant of a result proved in [1] 
with a slightly different estimate for the remainder). 


Corollary 5.6. There exists a universal constant 5 > 0 with the following property. 
Assume that 


and let A be the operator A : 


q — 11 14 + ||/ 3 ||i/M < S , 

Op (q{xi . For any u in L 2 , there holds 


i[\D x \l,A\u 

= ^(d x /3)d x {Au) + Op j(d x /3) 2 q"j\f\h l ^^u +R a u (5.9) 

where Ra satisfies 

\\Rau\\ L 2 < C5\\u\\ L 2. 

Proof. Denote by p the symbol p = q(x, Set M = 8 and write 

2 

i |D x |i H = Op ( ^ dfp) + R 0 + Rm, (5.10) 

ct—0 
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where 


M—l 

Ro : =Op(E l^ d t\^ d >)- 

a=3 

3 

For any 3 < a < M — 1, the symbol |£|z d x p is a linear combination of terms of 

the form m(x, ^)b(x)e^ k ^ ^ where m is of order 0 (that is d^m(x,^) < |£| _i ) and b(x) 
is of the form (d x °q)(d x 1 (5) ■ ■ ■ (d x m /3). It follows from statement ii ) in Lemma 5.4 
that Ro is an operator or order 0, satisfying 

\\Rou\\ L 2 < C(5)S \\u\\ L 2 . 

We now estimate the operator norm of Rm- By applying (5.8) with s = So = 1, 
m = 1 and M = 8, then JC 2 ( m+r+S0+ i' ) < IC s+ M+m +2 = ^12 and the inequality 
simplifies to 

\/u € L 2 ( T), || Rm \D x \ u\\ h i < C(l)/Ci 2 ||ti||#i. 

Now we estimate the L 2 -norm of Rmv for v in L 2 . We can assume without loss of 
generality that v has zero mean (since RmC = 0 for any constant C) and then set 
u = lAz-l” 1 v. The previous inequality yields 

II-Rm^IIl 2 < \\Rmv\Ih 1 = II Rm \Dx\ u\\m < C(S)S ||u|| L 2. 

Therefore one has 

|| (Ro + Rm)u\\ 1,2 < C(S)5 ||w || L 2 ■ 

It remains to study the sum for 0 < |a| < 2 in the right-hand side of (5.10). One 
can split this sum into two symbols such that the contribution of the first symbol is 
the two terms in the right-hand side of (5.9) while the other symbol is of the form 

Q(x, £)e 1 ^ 2 ^ with Q of order 0. Therefore the contribution of the second symbol 
can be estimated by means of Lemma 5.4, so it can be added to Ro + Rm to obtain 
an operator Ra satisfying the estimate in the statement of the lemma. □ 


Notation 5.7. Set 


N ll^llc°([0,T];// s O- d ) + ll^3|lc°([0,T];£(L 2 )) i 

where so is the large enough integer which appears in the definition of J\f (see As¬ 
sumption 5.3) and d is an absolute number independent of so (as in the statement 
of Proposition 5.1). 


We now chose /3 under the form f3o(t) + f3\(t, x) for some function coefficient (3o(t) 
to be determined later and with (3 4 = ^d~ 1 W. Then (3 is such that 

\d x p = § 4/31 = W. 

Recall from (5.6) that 


R 5 := A" 1 ( [dt, A] + R 4 A + Wd x A + i [L, A] ). 


(5.11) 


Now we split the last term as i[L,A] = i \D X | 2 , A + i L — \D X | 2 , A . Then 

it follows from the previous corollary that the remainder R§ (as defined by (5.6)) 
satisfies 


Rs ■= A- 1 ( [dt, A] - Op - J (d x P) 2 q) 


9 i 


+ R^.A H - % 


‘2|C| 

L — \D X \% , A 


- Ra 


(5.12) 


27 









where Ra is as given by Corollary 5.6. Recall that R 4 is an operator of order 0. On 
the other hand 


[d t ,A] = Op (d t p) = Op {[d t q + i\i\ 2 {d t P)q)e l ^^j. 

So one can write R$ under the form R 5 = R^J 2) + R^ ] where i?~ 1 ' /2 ' ) (resp. R^) is 
of order 1/2 (resp. 0 ), 


R^ /2) ; = A~ l 

Rf ] := A- 1 




R4.A — Ra H - ^ 


L — | D x | 2 , A 


+ Op(($g)e^)). 


We claim that 

11-^5 II C°([ 0 ,T];C(L 2 )) ~' (5.13) 

Indeed, Ra has already been estimated and, directly from its definition (see (5.5)), 
the Sobolev embedding ||<9 x W|| L oo < ||IT ||#2 and (5.3), one has ||^ 4 ||co([o,r];£(L 2 )) ~ 
JV 7 . The last term is estimated by means of Lemma 5.4 and to estimate the com- 

3 3 

mutator [A,L — \D X \ 2 ] we notice that L — |D X | 2 is a smoothing operator. 

Now, in view of (5.2) and (5.3) one has N' < N and hence ||i4°' ) 11 c<o ([ 0 t]-£(l 2 )) ~ 

N. 

It remains to prove that /3 and q can be so chosen that R^ 2 ' 1 = 0. To do so, we 
first fix fio(t) such that 


27 rdtfio = - J {dtp 1 + < ^(d x /3 1 ) 2 ^(t,x)dx, (5.14) 

where recall that fii = — ^ c?~ 1 LW, so that 

/ { dt P + = °- 

Now dehne q as q = e 7 where 7 is such that 

7 = |*||| ^(dtP + ^xP) 2 )- (5.15) 

(Notice that the previous cancellation for the mean implies that 7 is periodic in x.) 
With this choice one has R^ 2 ' > = 0. 

By combining the previous results, we end up with the following proposition. 


Proposition 5.8. Assume that So is large enough. Consider the operator 

A := Op (q(t,x,0z m ' Xmh ) 

with 

P = Po(t) + \ d x l W 

where Po determined by (5.14), and q = e 7 where 7 is given by (5.15). Then 
(dt + W d x + iL + Rfij A = A (pt + iL + Rp 

with 

11-^5 llc' 0 ([0,T];£(I/ 2 )) 

where JV is as defined in Assumption 5.3. 
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6. Ingham type inequalities 


As already mentioned, the controllability of the linearized equation around the 
null solution is based on (a modification of) the following Ingham’s inequality: for 
every T > 0 there exist two positive constants C\ = C\ (T) and C 2 = ^(T) such 
that, for all (w n ) n£ z E ^ 2 (Z;C), 


Ci ^2 \w n \ 2 

tiEZ 



2 

dt < C 2 E I w n 

nEZ 


Hereafter, ( w n ) n& z always refers to an arbitrary complex-valued sequence in i 2 { Z). 

For our purposes, we need to consider more general phases that do not depend 
linearly on t. For some given real-valued function j3 E C 3 (M), set 


Hn(t) = sign(n) £{n)t + /3(t)\n 


£(n) = (g + n 2 ) 2 |n| 2 tanhz (b\n\), 


with fiQ = 0 and sign (re) = n/\n\ for re 7 ^ 0. We recall that l is the symbol of the 
linear operator L = (g — c^) 2 ( 7 ( 0 ) 2 obtained by linearizing the water waves system 
around the null solutions, see Section 2.2. We begin by proving a lower bound 
which holds for any T > 0 provided that the functions contain only large enough 
frequencies. 


Proposition 6.1 (High frequencies). Let T > 0. Let \dtfi\ < \ tanhz ( 6 ) and \d 2 /3\ < 
1 for all t E [0,T]. Then there exists Nq> 0 such that, for all N > No, 


T 

2 " 


2 K 

raSZ 

|n|>7V 


< 


y: w n e 

nGZ 
\n\>N 




dt. 


( 6 . 1 ) 


Remark 6.2. i ) For T small, one can take Nq = CT~ 2 ~ £ for some e > 0. See 
(6.8) for more details on this estimate, ii) For ||c^/3||£oo small enough and T large 
enough, the result holds with No = 0 . 


Proof. Splitting the sum into re = m and re 7 ^ rre, we write 


r T 2 pT 

/ w n ei,i " W dt>T y2 \w n \ 2 + V W n m^ / e )) ^ 

J 0 Jo 


nEZ 
\n\>N 

We have to estimate 


Integrating by parts, 
AT(re, m) = 
and therefore 


nSZ 

\n\>N 


n^m 
\m\>N,\n\>N 


r T 

K(n, m) := / dt _ 

Jo 


1 t=T ,T 


Ah'nit) ~ h' m {t))\t=0 


\K(n,m)\ < n(n,m) := 
Since n(n,m) = we have 


f 1 n" — n" 

_l_ / c hUnit)-limit)) rn rin ^ 

io ®04» - iC) 2 

T |u" _ n" I 


h'n- h'm 


+ 


L°°([0,T]) Jo \h'n h'm\ 2 


dt. 


Y J ^nWmK{n,m) < i ( \w n \ 2 + \w m \ 2 )n(n,m) < ^ \w n \ 2 «(re,m). 


n^m 

Hence 


E 

n£ Z 
|n|>JV 


n^m 

2 


n^m 


w n e 


iUn(t) 


dt> 22 ( T ~ k ( 

nEZ mGZ\{n} 

M>W |m|>IV 


re, 771) 7U 


29 
























We have to prove that N can be so chosen that 

T- K(n,m)>^. (6.2) 

mEZ\{n} 

|m|>7V 

To do so, we use the following lemma. 


Lemma 6.3. Assume that \dt/3(t)\ < ^ tanhz (6) for all t. Let e € (0, \). Then 
i ) There exists a positive constant K e such that, for all integers N > 0 and all 
n € Z with \n\ > N, 


E 


m£Z\{n} 

\m\>N 


1 


hn ~ IT 


< 


k f 


l°°([o,t]) (l + iV)2' 


ii ) For all integers n,m. in Z with n / m, and all t, 

1 


<2tanh-t(fr) \d?/3\. 

I hn hm I 


(6.3) 


(6.4) 


Proof. Let us prove statement i). Since n{—n,m) = n(n, —m), we can assume, 
without loss of generality, that n > 0. Let \dt/3\ < ^tanha( 6 ), and note that 
tanh( 6 ) < 1 < 1 + g. Then for all n > 0 

13 isll3 3 13 

tanh2(6)n2 < £(n) < (1 + 5)2712, - tanh2(6)n2 < < - (1 + <7)2712. 

For m < 0, m / n, one has 

1 1 3 3 1 1 3 

1/4 - 4ml = h'n + h'-m > If tanh " ( b ) ( n * + H*) > 4 tanlD ( 6 ) (1 + \m\ 2 ) 
and therefore 


E 

m<—N 


hn hm 


c 


< 


C' 


i°°([ OTD m l"L 1 + H 3/2 ^ 1 + iV 


for some constant C' > 0. We now consider the case m > 0 and split the sum into 
two pieces. For rri > An with A := (36(1 + g)/ tanh(fr )) 1 / 3 one has \J m > 2 fi' n , and 


| hn hm | — hm 



hm 


1 . 

— 2 hrr 



3 

2 


which again leads to a convergent series 


E 

m> 0 

m>An , m>N 


hr 


hm 


< 


L°°([0,T]) 


c 

VTTn' 


It remains to consider the sum over all m > 0 such that N < m < An. Denote 
cr(n) := \J{g + n 2 )n. Then 

a (m) — a(n) a(m) 2 — a(n) 2 m 2 + n 2 + nm + g 

m — n (m — n)(a(m) + a(n)) a(m) + a(n) 

Using the elementary inequality ab < l(o 2 + 6 2 ), one has 

(er(n) + cx(m)) = \/ n 2 + g n + 1 /m 2 + g\Jnm < m 2 + n 2 + nm + g 

for all m, n > 0. Therefore 

|<r(m) — (j{n)\ > max{n, m} 
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m — n 


(6.5) 



for all m, n > 0. Now suppose that m < n, with m > 0, n > 1. Then 

/j,' n — // m = (cr(n) — a(m)) tanhz (bn) + < 7 ( 777 .) (tanh 2 (bn) — tanha (bm)) 
+ (n 2 — mz)dt/3 

> (cr(n) — a(m)) tanhz (bn) — (nz — m^)\dtP\ 

m\ 


>n 2 (n — m) tanhs (b) ^1 — 
> — tanhz (b)y/n(n — m) 


\fn(yfn + y/m) • v /tanh( 6 ) 


if \dt/3\ < |-^/tanh( 6 ). We deduce that 

\n' n — p' m \ > C max{n, m} \n — m\ 
for all m,n> 0, with C = lj tanha (b). Now, for n > 1, we obtain 


E 


m> 0, m^n 
N<m<An 


1 


Tn - /4 


< 


1 


L°°([0,T]) 


E 


1 


Cy/n A ' / \n — m\ 

m> 0, mfn 
N<m<An 


< 


clog(cn) 


n 


( 6 . 6 ) 

(6.7) 


for some c > 0. For n > 1 and n > N, one has n > ^(1 + IV), and 

clog(cn)n _ 2 < C' e n _ 2 +£ < C £ 22 _£ (l +iV)~2 +e 

for e € (0, rj), for some C £ > 0. On the other hand, for n = 0 the first sum in (6.7) 
is zero because it has no terms. Thus the first sum in (6.7) is < C e ( 1 + N)~z +e for 
any n > 0. This completes the proof of statement i). Statement ii ) is proved using 
( 6 . 6 ). □ 

The previous lemma and the definition of n(n, m) imply that 


< \ ^ 2 A £ 

> ac (n,m) < - i- 

\m\>N 

Hence (6.2) is satisfied provided that 


(1 + T\\d 2 (3 \\l°°) 


4 Kr 


(l + T\\d?/3 \\l°°) < (1 + iV)^ _£ , 


and Proposition 6.1 is proved. 

From ( 6 . 6 ), applied for /3 = 0, we deduce that 


( 6 . 8 ) 

□ 

(6.9) 


| £(n) — l(m) | > C^niaxjn, m} \n — m\ 
for all rri. n > 0, with C = \ tanha (b). 

We now prove upper bounds. By contrast with the previous proposition, we shall 
see that these estimates hold for any function (not only for high frequencies). Also, 
a key point for later purpose is that one can add some amplitudes ( n depending on 
time (and whose derivatives in time of order k can grow with n as \n\ k / 2 ). 

Proposition 6.4. There exists C > 0 with the following property. Let T > 0. Let 
\dt/3\ < ^tanh 2 (&), and \d k fi\ < 1, k = 2,3 on [0,T]. Then, for all (w n ) € 1? 2 (Z;C), 


^ ^ IVnCnif ) f 




nEZ 


dt<CM(() 2 (l + T)^2\Wr 

nEZ 


where 


l\/r(/-\ II/- II I 11 Cn 11 L°° . || ^ Cn ||l° 

M(C) := sup HCnlU 00 + sup —r=== + sup || 

nEZ nEZ y 1 \Tl\ nEZ ' l^l 


( 6 . 10 ) 

( 6 . 11 ) 
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Proof. Splitting the sum into n = m and n / m, we write 

dt = yz([ \(n(t )\ 2 dt^\w n \ 2 + ^2 w n m^E(n,m) 
nEZ ® n^m 

with 

E(n, m) := [ T dt. 

Jo 

The first sum on the right-hand side is easily estimated. It remains to bound the 
sum for n/m, Integrating by parts twice, one has 

E(n, m) = [ fe lh dt = [e lh (-ifp + f'p 2 - fh''p 3 )\ ^ 

Jo 

+ [ e ih (f"p 2 - 3 fh"p 3 + Zfh" 2 p A - fh'"p 3 ) dt 

Jo 

with 

/ := CnCmj h '.= p n Pmi P := T~ ■ 

Pn ~ Pm 

Thus \E(n,m)\ < e(n,m), where 

e(n, m) := 2 ||/p|| L oo + 2\\f'p 2 \\ L oo + 2\\fh"p 3 \\ L °° ( 6 . 12 ) 

+ t(||/V||l« + s\\f'ti'p 3 \\L°° + 3||//i"V iu°° + \\fti"p 3 \\ L oo). 



We have to estimate the sum X)mez\{n} e ( n > m ), uniformly in n. First, we note that 

\\d?((nU)\\L~> = ||^/||l» < {(1 + M)* + (1 + |m|)l} fc M(C) 2 , fc = 0,1, 2. 

We have already seen in (6.4) that \h"p\ < 2\d 2 jJ\. Similarly, \h"'p\ < 2\df/3\. Also, 
applying (6.3) with N = 0, e = we deduce that X)mez\{n} \\P\\ L oo < C for some 
absolute constant C. Therefore the first, the third and the last two terms in (6.12) 
(i.e. those with /) are all bounded by CM(Cf) 2 { 1 + T). The remaining three terms 
of (6.12) are also bounded by CM(() 2 (1 + T) provided that 


E 

raEZ\{n} 


\n\ + \m\ 

{P'n - P'mY 


L°° 


< c 


(6.13) 


for all n E Z, for some C independent of n. The bound (6.13) is proved using the 
same splitting and estimates as in the proof of Lemma 6.3. □ 


By combining the two previous propositions with an induction argument (follow¬ 
ing [ 8 , 19, 37]), we now deduce the following result. 


Proposition 6.5 (Sharp Ingham type inequality). Let T > 0. Then there exist two 
positive constants C(T ) and 6(T) such that, if 

m\ x := sup \(dtP,$P,d$P)\<6(T), (6.14) 

te[o,T] 


then, for all (w n ) E ^ 2 (Z; C), 


C(T)^K| 2 < 

nEZ 



nEZ 


2 

dt. 


Proof. This proposition will be deduced from Proposition 6.1, the following claim 
and an immediate induction argument (with a finite number of steps). 
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Claim 6.6. Consider two subsets A. A of Z with A' = A U {iV} for some N € Z, 
and with \n\ > |JV| for all n in A. Assume that for every T > 0 there exist two 
positive constants 5(T) and K(T ) suc/i f/iat 


H0II* < m 


K(T) Y, K| 2 

nGA 


< 



w n e itin ® 


2 

dt. 


(6.15) 


Then for every T > 0 there exist two positive constants 5'{T) and K'(T) such that 


Ml x<6'(T) 


K'(T) Y K| 2 < 

nGA 1 


r T 

/ V w n e l ^ {t ) 

Jo nGA' 


dt. 


(6.16) 


Let us prove the claim. Introduce 

fit) := Y ^ n(t) , fit ) := Y Vne^, hit) := Y 

nGA tiGA' nGA' 

so that f = f + WNe l,1N , fi = e~ l ^ N f = fe~ tfJ,N +wn, and 


[ T \ hit)\ 2 dt= [ T \f'it)\ 2 dt 
Jo Jo 


JVn(t) 


dt. 


Y Wne ~ 

nGA' 

We prove that there exist two constants C i, Co (both depending on T) such that 
Ci V \w n \ 2 < [ \fit)\ 2 dt, C 2 \w N \ 2 < f \f'it)\ 2 dt. (6.17) 

nGA J ° 

Then (6.17) implies the second inequality of (6.16) with I\'{T) := ~ min{Ci, C 2 }. 
Let us begin with the first inequality of (6.17). Let t := 4 min{ 1, T}, and remark 
that 

ihit + rj) - hit)) dr] = Y Oft), (6.18) 

nGA 

(notice that the sum is over A and not A) with 

Onit) = 1 : ^ e KUri(t+7j)-ll n (t)-tl N (t+ri)+fl N (t)) _ ^ 

Assume that n,N are positive. We split 9 n = c n + ( n , where c n is a constant, 
independent of time (such that c n = 9 n for /3 = 0), and Cn is defined by difference, 
namely 


f 

Jo 


C77, • — 


== f (. 

==/ 


P i[£(n)—£(N)]T _ 1 

Mn)-emv _ 1 ) dr] = - T, 


i[iin) - £(N)\ 


Q n ._ I e i[£(n)-e{N)]ri ^ e i[P(t+r))-P(t)](y/n-VN) _ ^ 

Now we use the following elementary inequality: there exists an absolute constant 
Co > 0 such that, for all ileK, 

— l — id | 2 > co min{$ 2 , il 4 }. 

This inequality holds because |e l1? — 1 — id\ 2 = (1 — cos d) 2 + (i? — sin d) 2 is positive for 
all d h 0 and it has asymptotic expansion d 2 + o($ 2 ) for |i?| —>• 00 , and \d A + o(i? 4 ) 
for d —>• 0. We apply this inequality with d = [fn) — £iN)]r, and, using (6.9), we 
get 

|c n | 2 > cr 4 

for some c > 0 (note that min{r 2 ,r 4 } = t 4 because, by assumption, t < 1). 

It remains to estimate Cn and its derivatives. From the definition, 

ICnl < 2r, \d t Cn\ < 2\\d t f3\\L°°Ty/n, \d 2 Cn\ < 4(||a t /3|||oo + ||(9 2 ^|| L oo)rn. 
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However, we need a sharper bound on £ n which shows that ( n is small when /3 is 
small. Such a bound could be easily obtained by estimating \e l f —1| < |/|. However, 
this would make appear an extra factor y/n. Instead, we integrate by parts to obtain 


Cn(t) = - 


e i[e(n)—£(N)]T 


i[£(n) - £(N)j 

Jo 1 


J[P(t+T)-P(t)](y/n-VN) 


- 1 


r e i[e(n)-i(N)]ri 


- l) dr], 


i[£(n ) - £(N)] 

and it is easily checked, using (6.9) and the bound |/3(f + r) — /3(t)\ < T\\dt/3\\L°°, that 
ICnl < Crll^H L oo . By combining the previous estimates, we have M(() < Cr||/3||x 
where M(£) is given by (6.11), and C is independent on T,t. 

Set F(t) := YlnGA w n eitJ ‘ n ^&n(t) and split F = F\ + F 2 with 

F^t) := Y, uneven, F 2 (t) := w n e^C n (t). 

nGA nGA 

Since |c n | 2 > cr 4 , the assumption (6.15) implies that, if \\/3\\ x < 5(T — r), then 

ct 4 /\ (T -t)^ 2 \w n \ 2 < K(T -t)^2 \ w nC n \ 2 < / |-Fi (t)\ 2 dt. 

nGA nGA 0 

On the other hand, Proposition 6.4 applied with M(£) < Cr||/3||x implies that, if 
||/3||x < ^tanha( 6 ), then 

[ T T \F 2 (t)\ 2 dt < C 0 t 2 ||/?|||-(1 + T - r) V \w n \ 2 

Jo nGA 

where Co is independent of T, r. Therefore, if 

40)T 2 ||/1||^(1 + T - r) < ct a K{T - r), (6.19) 

then f^ T \F 2 \ 2 dt < \ \Fi\ 2 dt, whence J^ T \F\ 2 dt > \ J^ T \F 1 \ 2 dt. By 
(6.18), this implies that 

pT-r 

|2 


jCT A K{T-t) \w n \ 2 < [ |C(t)| 

! nGA J ° 


dt 


< 


\ T ~ T f T ih{t + d)-hit)) 

Jo Jo 


dr] 


dt. 


The condition (6.19) holds if 


lx< 


Ty/cK{T — T) 
2\/Co(l + T) ' 


( 6 . 20 ) 


and we set 6'(T) as the minumum among ^ tanha (&), 6{T), and the constant on the 
right in (6.20). Moreover, 

T 7” 7” 2 

/ [ (flit + 7?) - hit)) dri dt 

Jo Jo 

< [ t [ \ihit + V) ~ h(t))\ 2 drjdt 
Jo Jo 

pT—r rr pT—r pr 

/ / \hit + r])\ 2 dr]dt 

Jo Jo 


< 2 r 


< 2 Tt 


+ 2 r 


r \hit)\ 2 dt=2TT r \f\t )\ 2 dt, 

Jo Jo 


f\hit)Y 

Jo 


drjdt 


and we infer that the first inequality in (6.17) holds with C\ = |ct 4 A'(T — r)T 
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Now we prove the second inequality in (6.17). We have \wn\ 2 = | f'(t) — f(t)| 2 for 
any t, and so 

\w n \ 2 = ^J^ \f{t) - f(t)\ 2 dt < \f{t)\ 2 dt +\m\ 2 dt). 

It follows from Proposition 6.4 (applied with = 1) that 

[ T \f(t)\ 2 dt<(l + T)C V \w n \ 2 . 

J ° n&A 

Using the first inequality in (6.17), we deduce that 

£m\ 2 dt< {1+ c ^ )c £\f\t)\ 2 dt, 

where C is the constant of Proposition 6.4 and C\ has been found above. Conse¬ 
quently the second inequality in (6.17) holds with C 2 = \TC\[C\ + (1 + T)C] _1 . 
We set K'(T ) = \ min{Ci, C 2 } and obtain (6.16). This completes the proof of the 
claim in the case of n, N positive. The other cases are analogous. □ 


7. Observability 


We now use the previous inequalities for sums of oscillatory functions in order 
to prove an observability property. In particular, we prove that it is sufficient to 
control the real part of the solution to bound the initial data. 


Proposition 7.1 (Observability). LetT > 0. Consider an open subsetu C T and a 
constant 0 < c < 1. Then there exist positive constants K,e\ such that the following 
property holds. Consider a pseudo-differential Ag with symbol exp (i/3(t, x)|£| 2 ) for 
some function [5 satisfying 

sup sup \(d t /3(t,x),d^l3{t,x),df/3(t,x))\ < S(T), 

te[0,T] xG[0,2tt] 


where S(T) is the constant in Proposition 6.5. Then for every initial data vo € L 2 ( T) 
whose mean value (vo) = 9 ^ f T vo(x) dx satisfies 

|Re(u 0 )| > c|(u 0 )| -£i ||u 0 |lz ,2 , (7.1) 


the solution v of 


dtv + iLv = 0, u(0) = vq, 


satisfies 



|Re(v4ou)(i, x)\ 2 dxdt > K 



|u 0 (x )| 2 dx. 


(7.2) 

(7.3) 


Remark 7.2. The condition (7.1) cannot be eliminated. To see it, consider the 
simplest case (3 = 0, so Ao = I, and consider a constant solution v(t,x) = C of 
(7.2). Then (7.3) holds for some I\ if and only if the real part of C is non zero. This 
suggests to assume that 

|Re(u 0 )| > c|(u 0 )| . (7.4) 

In fact, it is sufficient to consider the weaker assumption (7.1). The advantage of 
assuming (7.1) instead of (7.4) is used below (see (7.14)). 


Proof. Write 


v(t,x) = J2 a n e inx e ie{n)t , a n 

nEZ 


1 /' 2 V m O; 0 (x)rix, 
^ Jo 
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where l{ri) = (g + n 2 ) 2 (|n| tanh( 6 |n|)) 2 is the symbol of L. Then set w = Aqv , given 
by 


W( 


— > a e inx e i W n ) t +^i. t < x ) l n l I ) i 


nEZ 


For n E Z, set 

A n = ^(n)f + /3(t, x)|n|^, (j, n = sign(n)A„, c n (x) = a n e m3: . 
Since /r n = sign(n)A n and H- n = — fi n , we write 

2 Re tu = 2 Re a 0 + ^ c n e iA71 + ^ c^e~ iXn + ^ c n e* A ™ + ^ ~-~ iXn 


Cn.e 


n> 0 


n>0 


n< 0 


n<0 


2 Re c 0 + c ^ n + c -^ n + ' 


n> 0 


n< 0 


n>0 


n<0 


to obtain 


2 Re re = ^ 'y n e lfln with 7 n = < 

nEZ 


c n + c_n for n > 0 , 
2 Re cq for n = 0, 


c n + c_ n for n < 0 . 
Consider an interval wq = [o, fe] C u. By Proposition 6.5, 


r 0 J U) 


R e(w(t, x ))\ 2 dxdt > 


^0 Jo 


> 


R e(w(t, x ))\ 2 dtdx 

^p! Eh-wi 2 *, 


(7.5) 


JuJ ° n£ Z 

where C(T) is the constant given in Proposition 6.5. For 0 we write 

1 2 w |2 . w |2 . ^ —- Plinx 


so that 


\ln(x)\ z = \a n \ Z + \a- n \ z + a n a- n e Zmx + a n a_ n e 2mx , 
/ |7n(x)| 2 dx > |w 0 |{K| 2 + |a—n| 2 } 

J OJQ 


I tin I tt-n 


( 

[ e 2inx dx 

+ 

[ e~ 2inx dx 

V 

JuJQ 


JUJQ 


Now 


[ e 2inx dx 

— 

f e~ 2inx dx 

— 

sin(n (6 — a)) 



J LJQ 


n 


Moreover there is a small universal constant > 0 such that, for all 6 E (0, Aq), 


sin(x) 


V|x| > 5, 

We can assume that 0 < b — a < 60 , so that 

V„€Z*, ( 6 -a)- Sin( " (t, - 0>) 

As a consequence, for all n 7 ^ 0, 


< 


sin(<5) 


n 


> (b — a) — sin (6 — a). 


/ |7n(^)| 2 dx > c'dttnl 2 + |a_n| 2 ), 

JuJQ 


where c' := (b — a) — sin (6 — a) >0. Then, recalling that 70 = 2Reao, it follows 
from (7.5) that 

j 


f T J< 


Ke(w(t,x ))\ 2 dxdt > C(T) (b — o)| Re ao | 2 + — |a 


neZ\{ 0} 
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Now, using (x + y ) 2 > \x 2 — y 2 and (7.1), one has |Re(uo )| 2 > \ c 2 |(uo )| 2 
namely 

c 2 

I Rea 0 | 2 > — |a 0 | 2 - 2?re? 2 _, |a n | 2 , 

nEZ 


and therefore 



Re(w(t, x ))\ 2 dxdt > K E |a n | 2 , 

nEZ 


IM 


2 

L 2 > 


with K = C(T) nrin {(6 — a)(|c 2 — 27r£ 2 ), |c' — (6 — a)27r£ 2 }. If £i is small enough, 
then K > 0, which completes the proof. □ 


Corollary 7.3. Let T > 0, let u C T be an open subset and let 0 < c < 1. T/ien there 
exist positive constants £o,£i ,r,I\ such that the following property holds. Assume 
that 

( W(t )) = 0 

for all t G [0, T] and 


sup 

te[o,T] 


E I d t W (t) 


l<k< 3 


+ sup \\W(t)\\ H r <£ 0 , 
H te[0,T] 


and consider the pseudo-differential operator A, given by Proposition 5.8, with sym- 
bol q(t,x,f)ex.p(i(3(t,x)\f\ 2 ). Then for every initial data vo G L 2 (T) whose mean 
value (vq) = ^ f T vo(x)dx satisfies 

|Re(u 0 )| > c|(u 0 )| — £\ ll^ollz ,2 , (7.6) 

the solution v of 

v(0) = v 0 , (7.7) 

/*27T 

> K / |uo(x)| 2 dx. (7-8) 

Jo 

(The constants £q,Ei,K depend on T,c, while r is a universal constant.) 


dtv + iLv = 0, 


satisfies 


f\ 

1 0 Jlj 


I \Re(Av)(t,x )\ 2 dxdi 
J U) 


Proof. Split A as Ao + A\ with 

A 0 := Op (exp (*/3(i, x)|£|^)), Ai := Op ({q(t,x,£) - l)exp (i/3(t,x)\^)). 

The contribution due to Aq is estimated by Proposition 7.1. Notice that, for £o 
small enough, the smallness assumption on j3 of Proposition 7.1 is satisfied because 

sup sup \(d t (3(t,x),d 2 (3(t,x),df(3(t,x))\ < sup V \\d?W(t)\\ Hl < £ 0 . 

tE[0,T] rrE[0,27r] tE[0,T] ^<^.<3 

On the other hand, it follows from the definition of q and /? and the estimate given 
by statement ii) in Lemma 5.4 that A\ is bounded from L 2 onto itself, with an 
operator norm of size 0(||IR||#r) = 0(e o). Then 

[ [ |Re(Aiu)(f,.x )| 2 dxdt < [ ||Riw(t )|| 2 2 dt < [ £o \\v{t)\\ 2 L2 dt. 

JO Juj Jo Jo 

Since \\v(t)\\L 2 = Ho)|| L 2 , by taking £o small enough, the desired estimate follows 
from the triangle inequality. □ 


We now want to deduce an observability result for equations of the form 

dtw + Wd x w + iLw + IZw = 0, 

where 1Z is an operator of order 0. In the appendix we prove that the Cauchy 
problem for this equation is well-posed (see Lemma B.3). 
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Corollary 7.4. Let T > 0, wCT be a non-empty open domain and let 0 < c < 1. 
Then there exist positive constants £ 2 , £ 3 , r, K such that the following property holds. 
Assume that 

(W(t )) = 0 

for all t £ [0, T] and 

sup + sup \\W{t)\\ Hr + sup \\n(t)\\ c{L2) < £ 2 . (7.9) 

te[o,T] 1 |^ 3 l te[o,r] te[o,T\ 

Then for every initial data wq £ L 2 (T) whose mean value (wo) = ^ f T wq (x) dx 
satisfies 

|Re(u> 0 )| > c|(u;o)| - £3 ||wo|lz, 2 , (7.10) 

the solution w of 

dtw + Wd x w + iLw + TZw = 0, w(0) = wo (7.11) 

satisfies 


r 2ir 


' 0 j LJ 


|Reu;| dxdt > K I |u?o(;r)| dx. 


(7.12) 


Remark 7.5. Corollary 7.4 also holds for data at time T, that is: If wq € L 2 ( T) 
satisfies (7.10), then the solution w of 

dtw + Wd x w + iLw + IZw = 0, w(T) = wo (7-13) 

also satisfies (7.12). Note that the datum in (7.13) is at time T instead of 0. To 
prove it, notice that the function w(t,x) := w(T — t, x ) satisfies 

—dtw + Wd x w + iLw + IZw = 0, 

where W(t.),TZ(t) stands for W(T — t),LZ(T — t). Since W and TZ satisfy the same 
assumptions as W,1Z, one can apply (7.12) with w replaced by w, noticing that 


|Reu?r dxdt = 


I Re w r dxdt. 


' 0 j LJ 


' 0 J LJ 


Proof of Corollary l.f. It follows from Proposition 5.8 that there is a change of 
unknown w = Av such that v satisfies an equation of the form 

dtv + iLv + D\v = 0, 

for some operator 91 of order 0, satisfying ||91(f)u|| i 2 < Ce 2 \\v \\ L 2 for all t £ [0, T], 
By a perturbative argument, we shall deduce observability for this equation from 
observability for the equation without 91. To do so, split v as v = iq + V 2 where vi 
and V 2 are given by the Cauchy problems 


dtv 1 + iLv 1 = 0 

wi( 0 ) = Vq 


I d t V 2 + iLv 2 + 91u 2 = — 9!ui 
1 ^ 2 ( 0 ) = 0 

and vq := u(0) = (>l _ 1 r/;)(0). We begin by estimating v\, claiming that its initial 
datum vq satisfies the hypothesis (7.6) of Corollary 7.3, which is 

|Re(u 0 )| > c|(u 0 )| -£i INIIl 2 (7.14) 

(where £\ is given in Corollary 7.3). To prove (7.14), we write v = w + (/ — A)v to 
obtain, at time t = 0 , 

|Re(u 0 )| = |Re(w 0 ) + Re((/ - A)u 0 )| > |Re(u;o)| - \{(I - ^4)^o>I • 

Thus, using the assumption (7.10), 

|Re(u 0 )| > c|(u;o)| - esll^ollz, 2 - |((7 - -4)u 0 )|. 
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Since w = v + (A — I)v, we have (wq) = (no) + ((A — I)v o), and 

|Re(n 0 )| > c|(n 0 )| - (c+l)\((A - I)v 0 )\ - e 3 |ko|| l 2 - 

By (7.9), |((A—7)no)| < C'eaH^oHz ^ 2 (see Lemma 7.6 below). Also, ||^oHi 2 < C'II^oIIl 2 
because A is bounded on L 2 (see Lemma 5.4). Thus 

|Re(n 0 )| > c|(n 0 )| - ((c+ l)Ce 2 + Ce 3 ) \\v 0 \\ L 2 , 

and the claim is satisfied if e 2 ,e 3 are small enough. As a consequence, from Corol¬ 
lary 7.3 we deduce that 


r 2 TC 


|Re(Aui)| dxdt > K / |no(x)| dx. 


(7.15) 


'0 J uj 


On the other hand, it follows from (B.ll) (applied with V = 0, c = 1 and R = 91) 
that 

11^2 llc°([0,T];L 2 ) — ^ 11^1 IliLlO,T];L 2 ) • 

Since \\‘d\(t)v \\ L 2 < C'e 2 ||n|| i 2 , by using (7.9), we find that the last quantity is 
bounded by C^TH^oll^. Since A is bounded on L 2 , we deduce that 

[ [ |Re(Au 2 )| 2 dxdt < [ ||Au 2 (t )||^2 dt <T sup ||Au 2 (t )|| 2 2 

Jo Juj Jo [o,T] (7.16) 

< CT |k 2 |lc , 0([0,T];L 2 ) — CT 3 £ 2 \\V 0 \\ 2 L2 . 

Using the elementary inequality (x + y) 2 > \x 2 — y 2 , for e 2 small enough we get 

K ^ 


Since Av = w and 

rT 


rl r Tf /*Z7T 

/ / |Re(An )| 2 dxdt > — |no(a :)| 2 

Jo Juj 4 Jo 

L 2 = ll^ollz , 2 < C'||no||L 2 , we obtain 


rT r r2n r2n 

/ / |Reu ;| 2 dxdt > — |uo(x )| 2 dx > K' / |tco(a ;)| 2 dx, 

Jo Ju) J Jo Jo 

which completes the proof. 

Now we prove a technical result used in the proof above. 


□ 


Lemma 7.6. Consider a pseudo-differential A with symbol q(x, £) exp(i/3(x)|£| 2 ). 
There exist universal positive constants 6, C such that, if ||/ 3||^3 + |g — 1 | 3 < 6 , then 
\{(A- I)u)\<C5\\u\\ L 2 for allue L 2 (T). 

Proof. Like in the proof of Corollary 7.3, we split A = Aq + A\, with 

A 0 := Op (exp (i/3(x)|£|2)), Ai := Op ((g(x,f) - l)exp (i/3(®)|£|3)). 

Directly from statement ii) in Lemma 5.4 we have ||^Ai ||/:(z. 2 ) < C5 , whence |(Aiu)| < 
CJUull^. To estimate (Aq — I), let u(x) = Ylnez u n£ mx , and calculate 


/ (A 0 - I)u dx = V u n c n , c n = I e l ( nx +\nf 3{x)) dx _ 
J T n*0 

Integrating by parts, one has 


/ 
J T 


Cn — 


/- 

J T l 


hi x ))\ 


dx = 


f 

J T 


-id xx / 3 (x) e i{ - nx+ \ n ^h(x)) 

s I — 

T |n| 2 (l + |n| ?n- l d x (3(x))‘ 2 


dx 


i{n + |n | 2 d x f5{x)) 
so that, for \d x (3\ < 1 / 2 , 

\c n \ < C\\P\\ H 2 \n\-l Vn€Z\{0}. 

Thus (JO |c n | 2 ) 2 < C\\( 3 \\ H 2 , and by Holder’s inequality the lemma follows. 


□ 
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8. Controllability 


Consider an operator of the form 

Q := d t + W3 X + iL + R, 

where W is a real-valued function and R is an operator of order 0. In this section 
we study the following control problem: given a time T > 0, a subset w C T and an 
initial data Wi n E L 2 (T), find a (possibly) complex-valued function / E C°([0, T]; L 2 ) 
such that the unique solution w E C°([0, T]; L 2 ) of 

Qw = Xu Re /, w(0) = w in (8.1) 

satisfies w{T ) = 0. We study this control problem by means of an adaptation of the 
classical HUM method. We need to adapt the standard argument since we want to 
prove the existence of a real-valued control, while the unknown is complex-valued. 
In particular, for this reason, one cannot obtain w(T) = 0. We prove instead that, 
for any real-valued function M such that the L°°-norm of M — 1 is small enough, 
one can find a control such that w(T,x ) = ibM{x ) for some constant b E M. We 
remark that, given / and Wi n , the existence of a unique solution w to (8.1) is proved 
in the appendix, see Lemma B.3. 

We prove not only a control result but also a contraction estimate , which is the 
main technical result of this section. This means that we estimate the difference of 
two controls / and f associated with different functions W, W' or remainders R, R'. 
This contraction estimate is the key estimate to prove later that the nonlinear scheme 
converges (using a Cauchy sequence argument). To prove this contraction estimate 
we introduce an auxiliary control problem which, loosely speaking, interpolates the 
two control problems. Since the original nonlinear problem is quasi-linear, a loss of 
derivative appears. This means that to estimate the C 0 ([0, T]; L 2 )-norm of f — f 
we need to have a bound for the C'°([0, T]; IL 1 )-norms of / and f. This is why we 
prove and use a regularity property of the control, namely the fact that the control 
is in C'°([0, Tj; whenever the initial data Wi n is in This regularity 

result is proved by an adaptation of an argument used by Dehman-Lebeau ([17]) 
and Laurent ([29]). Before stating the result, we recall the definition of the adjoint 
operator Q*, namely 

Q* = -Q, Q:=d t + wd x + iL + n , n-.= -r*+ (d x w). ( 8 . 2 ) 

Proposition 8.1. Consider an open domain w C T. There exist r and six increasing 
functions Fj\ M(j_ —>• (0 < j < 5), satisfying limr-^o (T) = 0, such that, for 

any T > 0, for any real-valued function M E iL 3 // 2 (T) with ||M — l||i,°o < Fq(T), 
the following results hold. 

i) Existence. Consider R E C°([0, T]; C(L 2 )) and a function W satisfying 

I W(t,x)dx = 0 

■It 

for any t E [0, T]. Assume that the norm 

IIWii)lkr:= £ 119? W \\c°([0,T]-,m) + |[R llc°([0,T];RU + ll^llc°([0,T];£(L 2 )) 

l<fc<3 

satisfies 

\\(W,R)\\ r>T <MT). (8.3) 

Then there exists an operator @m,t ■ L 2 —>• C°([0, T]; L 2 ) such that for any Wi n € L 2 , 
setting f := &M,T(yJin), the unique solution w E C°([0, T]; L 2 ) of 

Qiu = Xuj Re /, w( 0) = Wi, 
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(8-4) 


satisfies 


w(T , x) = ibM(x) 


(8.5) 


for some constant and 

ll/llc°([0,T];Z2) < ■ ( 8 - 6 ) 

ii) Uniqueness. For any Wi n £ L 2 (T) and any T > 0, @m,t (Wi n ) determined 

as the unique function f £ C 0 ([0, T]; L 2 (T)) satisfying the two following conditions: 

(1) There holds Q* f = 0 and Im f T M(x)f(T,x) dx = 0. 

(2) The solution w of (8.4) satisfies (8.5) for some constant b £ R. 

Hi) Regularity. Let p £ [0,3/2] and consider Wi n € 7R l (T). If 

\\(W,R)\\r,T + ll^llc°([0,X];£(/f>)) — F\{T), (8.7) 

then Qm,t{w in) is in C°([0, T]; H^ l (T)) and 

||@M,r(wi ri )|| c o([ 0)T ] ;//Al ) < • ( 8 - 8 ) 

iv) Stability. Consider two pairs (W, R) and ( W',R'), where (W, R) is defined 
for t £ [0, T\ and satisfies (8.7) with p = 3/2, and ( W' ,R') is defined for t £ [0, T'] 
and satisfies (8.7) (with p = 3/2 and T' instead of T). Denote by ®m,t and ®' M T , 
the operators associated to these two pairs. Consider the time-rescaling operator T 
defined by 

(Th)(t) := h(Xt), A:=^, (8.9) 

and let W := TW, R := TR, namely R(t) = R(\t). Then, given any Wi n £ L 2 (T), 

ll@A l,T'i w in) - 7'0M,T('Wm)||cO([O,T'];L 2 ) 

< ^ DfiT) 0 1 ~ ^ + 11^ _ ^ llc°([o,r , ];_H' 2 ) + \\R' ~ -R||c°([o,t'] ; £(l 2 )))- (8-10) 

v) Dependence in M. Consider M,M' in 77 3 / 2 (T) with ||M|| LOO + ||M'|| LOO < 
R 0 (T). If \\{W,R)\\ r}T < RiiT), then, for all w in € H 1 ^), 

||(©M,T — ©M',r)(^m)|I coq 0 ,T];L 2 ) — j~ r (J 1 ) W M ~ • ( 8 -H) 

In this section we often use the notation A < B to say that A < CB for some 
constant C depending only on T. The key result is the following lemma. 


Lemma 8.2. Introduce the space 


L 2 m := £ L 2 ( T; C) ; Im J M(x)ip(x) dx = 0 j. 

For any W{ n £ L 2 (T), there exists a unique f\ € I? M such that, 

V0i £ L 2 m , Ref f(t),(f)(t))dt =-Re(w in ,(f(0)), 

Jo 

where f and cj) are the unique functions in C^QO, T]; L 2 (T)) satisfying 

f Qf = 0 f Qfi = 0 

\f(T) = fi \<KT) = 0i, 

where Q is given by (8.2) (the existence of f and cf> follows from Lemma B.3.) We 
set 

®M,T(Win) ■= f ■ 


( 8 . 12 ) 


Moreover (8.6) holds. 
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Proof. The space L 2 M is an M-vector space. Introduce the M-bilinear symmetric map 
a(-, •) defined by 


a(fi,<f>i) := Re [ f Xu(x)Re(f(t,x))cf>(t,x)dxdt (8.13) 

Jo J T 

= / / Xu( x ) Re(/(i, x)) Re/>(t, x) dxdt. 

Jo J T 

This application is well defined and continuous. Indeed, it follows from the L 2 -energy 
estimate (see (B.ll)) that 

K/i,</h)| < [ [ |/| |0| dxdt 

Jo J T 

< T ll/llc°([0,T];L 2 ) ll^llc°([0,T];L 2 ) ^ C{T) ||/l || L 2 ||</l || L 2 • (8-14) 

Since Xu( x ) = 1 for x in an open subset uq C ui, one has 

a(/i,/i)> [ [ (Re /) 2 dxdt. 

J 0 J uj-\ 


If /i € L 2 m then Im f T M fi dx = 0 and we have 


Im 


/ /i( 

J T 


x) dx 


Im / (1 — M(x))f\{x) dx 


< || M - 1\\ LOO V2n 


from which (using | Rez| > \z\ — \ Ini 2 ;|) we deduce that 

|Re(/i)| > |(/i)| — \\M — l|| L oo v / 27t ||/i ||^2 ■ 


L 2 


For ||M — 1|| LOO small enough, one can apply the observability inequality proved in 
the previous section (see Corollary 7.4 and Remark 7.5) to conclude that 

C'i(T)||/ 1 ||i 2 <a(/ 1 ,/ 1 ). (8.15) 

On the other hand, (8.14) implies that a(/i,/i) < C(T) ||/i||^ 2 . Hence a(-,-) is 
a real scalar product on L 2 M which induces the norm iV(/i) = y/affi, /i), which 
is equivalent to the norm ||‘||x, 2 (tc) on ^M - Now, Lemma B.3 implies that the 
mapping (j>\ i —> (f)(0) is M-linear and bounded from L\j into L 2 and hence (f> i i— > 
A((/i) := —Ke(wi n , (f)(0)) is a bounded M-linear form on I? M . Therefore, the Riesz 
theorem implies that, for any M-linear form A on L 2 M , there is a unique f\ € L 2 M 
such that a(f\,(f>i) = A(cf>i) for all <f>i € L 2 M , together with 


II/iIIl 2 < 


l|A|| 

Ci(TY 


(8.16) 


Moreover ( 8 . 6 ) follows from (8.16) and the bound \\f\\c°([ 0 T]-L 2 ) ~ II/iIIl 2 already 
used. □ 


Proof of Proposition 8.1. Proof of statement i). We begin by proving that if M G 
i7 3 / 2 (T) then H 3 / 2 (T)nL 2 M is dense in ( L 2 M , ||-||i 2 ). To see this, let be the Fourier 
truncation operator defined by njvh(x) = S|j|<Ar hje l J x where h(x) = Y 2 ,j^z^ l j el2X - 
Given u € L 2 M , define un '■= jj^J-n(AIu). Since the operator n^r preserves the 
mean, one has that un € L 2 M . Moreover, since u € L 2 , one has Mu € L 2 ( T), 
n n (Mu) g C 00 ^), and hence M~ 1 U N (Mu) G H 3 / 2 ( T) since M~ l G H 3 / 2 ( T). 
Since (un) converges to u, this proves that PJ 3 / 2 ( T) n L 2 M is dense in L 2 M . 

Now let / be as given by the previous lemma. It is proved in the appendix 
that there is a unique solution w in C°([0, T]; L 2 (T)) of (8.4). Our goal is to 
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prove that w(T) satisfies (8.5). To do so we first check that (8.5) will be proved 
if Re(w(T), 0Q = 0 for all 0 i in l? M . Indeed, one has 


Re(w(T), </>i) = (Rew(T,x))Re0i(x)dx + / (Imw(T, x)) Irn 0\{x) dx = 0 


for all £ L 2 m . Therefore we obtain j (R ew(T,x))f(x) dx = 0 for any real-valued 
function / and f (Im M(x)~ 1 w(T, x))g(x) dx = 0 for any real-valued function g with 
f g(x) dx = 0. This implies that (8.5) holds. 

We now have to prove that Re(w(T), 0\) = 0 for any 0 i in L 2 M . By the density 
argument proved above, it is enough to assume that 0 i € l? M fl L7 3 // 2 (T). Given 
0i £ L 2 m n H 3 / 2 { T), let 0 € C°([0, T];H 3 / 2 (T)) be such that 


0.0 = 0, 0{T) = 0 1 . 


(8.17) 


Since Q = — Q*, multiplying the equation (8.4) by 0 and integrating by parts, we 
find that 

(w(T),0i) = (w(O),0(O)) + [ (xcj Re/ , 0) dt + [ (w , Q0) dt. (8.18) 

Jo Jo 

Notice that the integration by parts is justified since 0 £ G 1 ([0, T]; L 2 (T)). By 
definition of 0 the last term in the right-hand side vanishes and, by definition of /, 
the real part of the sum of the first and second terms vanishes. This proves that 
Re ( w(T ), 0 \) = 0, which concludes the proof of statement i). 

Proof of statement ii). Recall that Q = — Q* is given by (8.2). Consider 0± € L 2 M 
and denote by 0 the unique function in C'°([0,T];T 2 (T)) satisfying (8.17). As in 
(8.18), multiplying both sides of the equation Qw = XojRef by 0, integrating by 
parts one obtains (8.18). Since 0\ € L 2 M and w(T,x ) = ibM(x) for some constant 
b £ M, one has Re (w(T), 0\) = 0. Therefore, since Q0 = 0, 


Re 


[ (xujRef, 0) dt = -Re (w in ,0(O)). 

Jo 


Since Qf = 0 and f(T) € l? M by assumption, and since the function /i whose 
existence is given by Lemma 8.2 is unique, one deduces that f(T) = f\. Hence 
/ = &M,T(win) by uniqueness of the solution to the Cauchy problem (8.12). 

Proof of statement iii). We prove ( 8 . 8 ). In view of the energy estimate (B.3), it is 
sufficient to prove that ||/i||^ M is controlled by ||wi n ||^ M . We prove only an a priori 
estimate, assuming that f\ belongs to H^( T). To estimate ||/i||^ M , we adapt to our 
setting an argument used by Dehman-Lebeau (see [17, Lemma 4.2]) and Laurent 
(see [29, Lemma 3.1]). 

Consider the mapping 

S■ L\i —t L 2 (T), S : fi £ L 2 m h> / t—>■ w e -> w( 0) £ L 2 (T) 

where / and w are the unique functions in C'°([0, T]; L 2 (T)) successively determined 
by the backward Cauchy problems with data at time T 


Qf = 0 
f(T) = fi, 


Qw = XouRef 
w(T) = 0 . 


It follows from statements i ) and ii) that S is an isomorphism of L 2 M onto L 2 (T). 
As a result, with A M = (/ — d 2 )^/ 2 , one can write 


Hr = l|A A 1 /i || i 2 < ||5A M / 1 || 


L 2 ■ 
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(8.19) 


Now we have to commute S and A^. This amounts to compare (A^/, A^w) with 
( f',w ') defined by 

[Qf = 0 f Qw' = Xw Re f 

\f'{T) = A M /i, \w\T) = 0. 

We first estimate f — A^f and then deduce an estimate for w' — A^w. Write 

QW ~ A W) = [A**, n\f + [A M , W] 4 /, W - ^f)\t=T = 0, 

and use the energy estimate (B.13) to find that 

l|/ , - AM /|| c o ([0iT];i2) < ||[A^^]/ + [A^,W]a x /|| il([0;T];L2) . (8.20) 

Similarly, 

IK -A M Hlc°([o, t ]; l 2 ) < \W\\l^[0,T]-l 2 ) where ( 8 - 21 ) 

T := Re(/' - A"/) + [A'*, R]w + [A", W]d x w - [A", Xw ] Re /. 

By (8.20) and the obvious embedding C°([0, T]; L 2 ) C L 1 ([0, T]; L 2 ), we deduce that 

H-^11 £!([(), T];L 2 ) ~ || [A M , Xu] Re /|lc°([ 0 ,T];L 2 ) 

+ ||[A At ,7^]/|| c ,0([ 0]T ] ;i2 ) + || [A M , R]w\\c°([0,T];L 2 ) 

+ II[A M ,R^]<WIIcO([ 0 ,T];L 2 ) + \\[^iW}d x w\\ C 0(! [ Q^T]-,L‘ 2 ) ■ 

To estimate the commutators [A^,Xw] and [A M , W], we use the classical estimate 

s>^, 0 < /j, < s =>- \\[A^W]u\\ L 2 <K\\W\\ H s\\u\\ h »-i. 

On the other hand, to estimate the commutator [A^,TZ\ (or [A /x , R ]) we estimate 
separately A^IZ and 1Z A M . Recalling that 7Z = — R* + {d x W), we conclude that 

Ill'll£!([(),T];L 2 ) ~ l|/llc°([0,T];i?' 1 - 1 ) + a II (/> w )\\c°([0,T]-,HI*) 

where 

a := l|R'1lc , °([0,T];.ff 3 ) + l|R|lc°([0,T];£(H' i )n£(L 2 )) ’ 

Notice that a < ||(fU,i?)|| riT + ||-R|Ic , °([o,t];£(.H’ 0) where j|(B / ,R)|| rT is as defined 
above (8.3). 

Now, using the energy estimate (B.12), we have ||/||co([o,t];^) < 11/illifM- Using 
again (B.12) and the equation satisfied by w, we deduce that HHIcouq tyh^) ~ 

WhWHv- Thus - b y ( 8 - 21 )> we find 

Ih'-A^Hc°([0,T];£ 2 ) ~ ll/lll^-i + a ll/l||_ffM • 

In particular, at t = 0, we get ||u/(0) — A /x u)(0)|| L 2 < ||/i||b>-i +a||/i||^ M . Now, by 
definition, w'( 0) = SA ll f\ while A^ l w(0) = A^Wi n . Therefore, by triangle inequality, 



II SA»fi\\ L2 < || w in \\ Hlt + ||/i||^-i +a||/ 1 ||^ . 

(8.22) 

For fj, € [0,1], 

one has ||/i||^-i < ||/i||l 2 < W^in |l 2 , and therefore 



ll-SA^/ill^ < lluiinll^ + a ll/llljy#. • 

(8.23) 

Plugging this bound into (8.19) yields \\fi\\ H „ < ||wm||j*> + a\\fi\\m- 
small enough, we conclude that 

By taking a 


|/l lii> < \\ w in • 

(8.24) 


For fi G (1, |] we go back to (8.22) and note that ||/i||# M -i < because 

/a — 1 G [0,1]. Hence (8.23) holds, and we reach the same conclusion as above. This 
completes the proof of statement in). 
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Proof of statement iv). Given Wi n , let f\ and f[ in L 2 M be as given by Lemma 8.2, 
so that / := QM,T(win) and f := & MT ,(wi n ) are determined by the Cauchy prob¬ 
lems 

Qf = 0 on [0, T] ( Qff = 0 on [0, T'} 

f(T) = fi \f'(T') = f' l 

where 

q! := dt + iL + w'd x + n', n' ■.= ~{b!)* + (d x w'). 

Similarly, we denote Q' = — (Q!)* = dt + iL + W'd x + R'. By definition of /,// 
the unique solutions w € C°([0, T]; L 2 (T)) and w' £ C°([0, T']; L 2 (T)) of the two 
Cauchy problems 

Qw = ) 
w( 0) = 

satisfy w(T ) = ibM , w'(T') = ib'M for some b, b' £ M. The idea now is to introduce 
an auxiliary control problem. Let f" £ C°([0, T]; L 2 (T)) be the unique solution of 

Q/" = 0 on [0,T], f"(T) = f[, (8.26) 

so that f" solves the same equation as / and it has the same Cauchy data as f. 
Then introduce w" as the unique solution to 

Qw” = Xu,Re/" on [0,T], w"(T) = ib'M (8.27) 

and set w” n := w"(0). By uniqueness (see statement ii )) we deduce that f" is the 
control for the operator Q associated to w” n , that is 

f" = @M,T(w" n )- 

Then, by continuity (see statement i )) one has 

|| / / 11 C° ([0, T];L 2 ) ||®Af,T(^m ^in) \ | C°([0,T];L 2 ) ~ ||^m ^in||^ 2 - 

Let / := Tf and f" := Tf". Then 

11/ — f"\\c°([0,T']-,L 2 ) = 11/ — f"\\c°([0,T];L 2 ) (8.28) 

because 

\/h £ C°([0,T]-, L 2 ), ||T/i||co([ 0 ) t'];L 2 ) = INIc°([o,t];L 2 )- (8.29) 

It remains to estimate \\f— /"||c°([o,T'];L 2 ) and \\win — tcf n || L2 - We begin with f'—f”. 
Since f” solves (8.26), f" satisfies 

Qf" = 0 on [0, T '], f(T') = f[ 

where 

Q := d t + iXL + X Wd x + A K, 

and W := TW, 7 Z := T1Z (namely TZ(t) := 1Z(Xt)). By difference, one has 

Q(f - f") = f 0 , (/'-/")(t') = o, 

where 


'.u Re / on [0, T] 


U)ir, 


Q'w 1 = Xu Re f on[0,T'] 

w'(0) = Wi n 


(8.25) 


F 0 := (A - 1 ){iL + Wd x + K)f + (W - W')d x f + (K - 11') f. 

In order to apply the L 2 -energy bound (B.13), we estimate Fq. Using the regularity 
property of the control operator ©m,t (see statement in)) we have 


\\ L f'\\c°{[0,T']-,L 2 ) ^ ll/ , |l c o([ 0i T'];fl-|) ~ 


\Wi, 


Vln " 

II (IT - IU , )5 x / , || C 0([ 0jT /] ;i 2) < \\W - R // || C -0([ 0jr /].^l)||/ / || C -0([ 0)T /]. fl -l) 

< \\W - W'\\ c o^ 0T i]. H i' ) \\wi n \\ H i. 
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(8.30) 

(8.31) 


Similarly 


|| (^ - ^ / )/ , ||c°([0,T , ];£ 2 ) 

< 1172. - R , ||c 0 ([0,T'];£(L 2 )) ||^ 11 C° ([0,'Z’ / ] ;Z/ 2 ) 

< (ll-R - R'\\c°([ 0 ,T'}-,C(L 2 )) + \\W — ^Ilc 0 ([ 0 ,r'];_ff 2 )) ll^inll-L 2 

where R := TR (namely R(t) = R(Xt)). Using (B.13) we conclude that 

11/ ~ f"\\c°([o,T']-,L 2 ) 

/ ll^mll^-3 (|A — 1| + II W — W / ||c' 0 ([0,T , ];/f 2 ) + ll-R — R!\\c°{[0,T'}-,C(L 2 ))) ■ (8.32) 

It remains to estimate \\wi n — w” n \\ L 2 . Let w" := Tw". At t = 0 one has w'(0) — 
w"{ 0) = Wi n ~ w '( n ■ hence we study the difference w' — w". Since w" solves (8.27), 
w" satisfies 

Qw" = \xu Re /" on [0,T'|, w"(T') = ib'M, (8.33) 

where 

Q := d t + iXL + XWd x + XR. 

By difference, 

Q(w' - w”) = F, (w' - w")(T') = 0, 

where 

F := Xco Re(/' - A f") + (A - 1 )(iL + Wd x + R)w' + (W - W')d x w' + (R - R!)w'. 

In order to apply the L 2 -energy bound (B.13), we estimate F. First, f'—Xf" = A (/'— 
f") + (1 — A)/ 7 , and we have already estimated both f — f" (see (8.32)) and f. For 
the other terms in F we proceed as above, recalling that ||it/||c°([o,T'] ; l 2 ) 1$ llwmllz, 2 - 
Also, since w' solves the Cauchy problem (8.25), we deduce from (B.12) and the 
second inequality in (8.30) that H^HcoqoT']-// 3 / 2 ) ~ ll u ’m||//3/2. As a consequence, 
also ||-F||co([o,t'];L 2 ) is bounded by the term on the right-hand side of (8.32). Then, 
applying the energy inequality (B.13), we deduce that || w' — 'ih // ||c°([o,T , ];£ 2 ) satisfies 
the same bound. In particular, at time t = 0, this yields the desired bound for 

Win - Win = (' w ' - w"){R). 

Proof of statement v). We begin by introducing some notations which are used 
in the proofs of statements v ) and vi). As already mentioned, it follows from 
Lemma B.3 that there exists an operator Ex'- L 2 (T) —>• C°([0,T]; L 2 (T)) such that 
v = Et{v\) is the unique solution to the Cauchy problem Qv = 0 with data 
v{T ) = v\. Moreover 

II-Lt(^i)IIco([o,t];L 2 (T)) / II^iIIl 2- (8.34) 

Now recall that by definition 

ot(/i>i) := Re [ [ Xu(x)Re(E T (f 1 )) Et(4>i) dxdt. (8.35) 

Jo J T 

Also introduce the mapping A: L 2 ( T) —>• M defined by A(ui) = — Re(u;j n , Et(v i)( 0)) 
where Et(v i)(0 ) = Ex(vi)\t=o- It follows from Lemma 8.2 that there exist two 
functions /i € L 2 M and f[ € L 2 M , such that 

V(£i € L 2 m , a T {fi, 4>i) = A(0i), 

V</>i € L 2 m ,, ax{f[, 4>i) = A(t)>i). 

Then &M,T(wi n ) — ®M',T{wi n ) = Ft(/i — /()• In view of (8.34), to prove statement 
v) it is sufficient to estimate /i — /{. To do so, we need to compare elements in 
L 2 m and elements in L 2 M ,. Observe that, by definition of L 2 M , if tp € L 2 M , then 
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(M'/M)<p € L 2 m . Therefore ip\ := f\ — belongs to L 2 M and we can use (8.15) 
to deduce that 

2 


h - TfX 


L 2 


M' 


M' 


<a T [h- tt/;, fi ~ -tt/; =a T (h~ -rrf'nV 1 


M 


M 


M' 


M ' 


(8.36) 


Now write the last term as the sum (I) + (II) + (III), where 

(I) = a T (fi,<Pi)-a T > 

i 11 ) = a T (f[, i). 

(///) = ar(/i,<pi) - or (jjfi’Vi) ■ 

(Notice that both M/M' and M'/M appear.) Since (M/M')ipi belongs to we 
can write ar (/(, (M/M')ip{) = A((M/M')p\) to deduce that 


(/) = Afo) - A ( ^ | = A 


M' - M 


M' 


-<Pi 


so that \(I)\ < || M' - M\\ Loo \\w in \\ L2 ||^i||^ 2 - On the other hand, it follows from the 
easy estimates (8.14) and (8.34) that 


\(II)\ + \(III)\<\\M-M'\\ 
< \\M-M'\\ 


L°° 


h \\ L 2 mil L 2 


\Wt 


I l 2 II^iIIlz 


By combining (8.36) with the previous estimates we conclude that 


Ji M h 


< ||M - M'| 


L 2 


L°° ll u 'm|lL2 


Now write 


\h-f[\\ L2 <\\M-M'\\ LO 


m 


L 2 


+ 


f -—f 

M ,J 1 


L 2 


and 11 /( 11^2 < 11 Wi n 1 1^2 to complete the proof of statement v). 

9. Controllability for the paradifferential equation 


□ 


We now deduce from the results proved in the previous sections that the original 
equation introduced in Section 3 is controllable, together with Sobolev estimates for 
the control. 

Consider a paradifferential operator of the form 

P = dt + T v d x + ilk (T c d • ) + R (9.1) 

where R is an operator of order 0. Assume that P satisfies Assumption 3.1, so that 
as above V and c are real-valued, c — 1 is small enough and P satisfies the following 
structural property: 


d I* 

Pu real-valued =>■ — / Im u(t, x) dx = 0. 

dtJ T V ; 


(9.2) 


Introduce the norm 


II (c- 1,V, -R)||x 8 0 .s(T) II (c— l,<9iC, y)||cfO([0,T];ff s O) + ^2 11^ c llc'°([0,T n ];iY 1 ) 

k= 2,3,4 

+ ll^ y llcO([0,T];RQ + ll^llcO([0,T];£(R=)) + l|-R|l c , 0([aT];/;( ^ s +3 )) . (9.3) 

fc=l,2,3 

We recall that p is the symbol given by p := c“3 + Q xC ( se e (2.12)). 
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Proposition 9.1. Consider an open domain uj C T\ There exists sq large enough 
and for any s > s o there exist three increasing functions Tj : M)j_ —> M)j_ (1 < j < 3), 
with lim-T-H) Fj(T) = 0, such that, for any T € (0.1]. the following properties hold. 

*) If 

||(c-l,V,P)|| XSo , 8(T) < Ti (T), (9.4) 

then there exists a bounded operator 

0 8) r[(v; c, R )]: H s+ 1 (T) ->• C7°([0, T];H S+ 1 (T)) 
such that, for any Vi n € H s+ 2 (T) satisfying 


Im / Vi n (x) dx = 0 , 

J T 


setting f := © s> t[(V, c, R)](vi n ) one has 


'c°([o,T]-,H s+ i) - JF 2 (T) 
and the unique solution v to Pv = T p Xxc Re f, vu = q = satisfies 

v(T) = 0. 


Vjrt 1 3 


(9.5) 


n) Assume that the triple (c, V, P) satisfies (9.4) and 

ll^t 2c llc'°([0,T];H s O) + l|5t'^llc°([0,T];H s 0) + II^P||c , °([0,T];,C(.ff s )) < 1- (9-6) 

Let (c', V', R!) be another triple also satisfying the same (corresponding) bounds (9.4) 
and (9.6). Then 

||0s,t[(4^, C, R)](vin) - 0 SjT [(F / ,C / ,i? / )]Kn)||cO([O,T];^) (9-7) 

11 ^VTl 11 s _l ^ ^ 

< T-JfT) 2 l ll(c ~ c/, ^^ C “ C 0; V - Ollc'°([0,T];.ff s 0) + II R — R'\\c°([0,T]-,C(H s ))J- 

Proof. Let P be given by (9.1), with V,c,R satisfying (9.4). We begin by recalling 
how the various linear operators have been defined in the previous sections starting 
from P: 

P := Afo.sPA,;] = dt + T v d x + iL?(T c L*-) + R 1 = d t + Vd x + iL^(cL^-) + R 2 , 

P 3 : = 4>m _1 P4)~ 1 = d t + TV'S* + iL + P 3 , 

P := -(P 3 )* = d t + Wd x + iL + R 4 , 


where &,m,W are given in Proposition 5.1, 

Ri ■= A-h,sRA-h,s A [A/i,s> $t]A^ s + [Aft )S , Tyd x ]A h s + i[A/ liS ,L2(r c L2-)]A /ljS , 

R 2 u := i?iu + Tyd x u — Vd x u + i(L^T c L^u — L 2 (cL^u)), 

R 4 w := -R* 3 w + (4H)w, (9.8) 


and R 3 has a more involved expression, obtained in Appendix C. Moreover P = 
?n4 > -1 Lad). As a first step in the proof of Proposition 9.1, we study the control 
problem for P. 


Lemma 9.2. There exist so large enough and increasing functions Tj : Ml —>• M^j_ 
(j = l,2,3j, satisfying lim-r^o Pj{T) = 0, such that for any T > 0 f/ie following 
result holds. 
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i) If 


||(c — 1, dt.C, R)||c 0 ([ 0 ,T];.ff s O) + ^ 11 C 11 C° ([0,T] j-f/ 1 ) 

k= 2,3,4 

+ H^^HcOaO.T];// 1 ) + 11-^2 ||c°([0,T];£(Z, 2 )) < Fl(T), ( 9 - 9 ) 

k= 1,2,3 


then there exists an operator @t■ L 2 —>• C 0 ([0, T]; L 2 ) such that for any Ui n E L 2 , 
setting f := ©rC^m) one has 

ll/llc°([0,T];L 2 ) ^ ( 9 ' 10 ) 

and the unique solution u of 

Pu = Xu Re /, u(0) = n in 

satisfies u(T,x ) = ib for some b E R and all i£T, //, m addition, 

^ R2 ^C°([o,t]-,C(h%)) ( 9 - 1:L ) 

then 

11 11 3 

| I /| Io» (I 0,T I ;«4, £ ■ < 9 - 12 ) 

ii ) Assnme that (V, c, R 2 ) satisfies (9.9), (9.11) and 


l|dtR||c°([ 0 ,T];_ff 2 ) + H^llcOQO.T];/:^ 1 )) + ||^d?2 llc 0 ([0,T];£(L 2 )) < 1) (9.13) 

and consider another triple (V. c!, Rif) a ^ so satisfying the same (corresponding) 
bounds (9.9), (9.11) and (9.13). Then 


|| (@T - ®T)( u in)\\c°([0,T]-,L 2 ) 

11 'U'il'l 11 3 r 

~ PfT) l^ C _ c 'Wc°([o,T]-,H r + 1 ) + \\dtc — dtc'\\ c o([ 0iT \. H i) (9-14) 

+ \\V ~ R / ||c°([0,T];H 2 ) + 11^2 — R2llc°([0,T];£(L 2 ))}- 

Proof. Recall that the cut-off function Xu( x ) is supported on uj and Xu = 1 on the 
open interval uj± C uj. Consider another open interval u>2 and a cut-off function 
X‘i (x) such that 


(*) supp(x 2 ) C U 2 ; 

(ii) if supp(h) C uj 2 supp( < h _ 1 /i) C cui Vt E [0, T], V7i E L 2 ( T). 


(9.15) 


We want to apply Proposition 8.1 for Q = P$. The hypothesis (8.3) of Proposition 
8.1, namely the inequality || (W, d? 3 )|| r ,T < P±(T), follows from the assumption (9.9), 
by using (5.3) and (C.40) with a = 3/2. Hence, by the statement i) of Proposition 
8.1 (applied with T\ instead of T and X 2 instead of Xu), given Wi n € L 2 ( T), the 
unique solution w of the Cauchy problem 

P 3 W = X 2 ReC/ 2 ) Vi E [0,Ti], w( 0) = w in (9.16) 

satisfies w(Ti) = ibM for some real constant b if we choose f 2 = @M,Ti(wm)> where 
©m,Ti is the operator given by Proposition 8.1, and the function M will be fixed 
below in this proof. Also, by ( 8 . 6 ), /2 satisfies 

II/2 llc°([0,Ti];L 2 ) < ' 
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(9.17) 






Moreover, if (9.11) also holds, then, using (C.39) with a = 3/2, we deduce the bound 
(8.7) for W, with p = 3/2. Therefore, by the statement Hi) of Proposition 8.1, 


II /oil 3 < 

n C°([0,Ti];H?) - 


m, 


3 

1 H? 


MTi) 


(9.18) 


Now let Ui n £ L 2 (T) be given and define Wi n £ L 2 ( T) by Wi n := 5>|t=o Ui n . We apply 
the previous argument and obtain a function w satisfying (9.16) and w(T\) = ibM. 
Set u := 3> _ 1 w;. Since Pu = m^~ 1 P 3 ^u, it follows from (9.16) that 


Pu = m<5> l (x 2 Re(/ 2 )) Vf £ [0, T], u(0) = u in , (9.19) 

and u(T) = ~= T (ibM). Then we set / := m 4 > _ 1 (x 2 / 2 )) namely we define 


f = ® T (u in ) ■= m<$> 1 (X 2 ©M,Ti(^Wm)), (9.20) 


where &Ui n = &\ t= 0 Ui n . By the assumption (ii) in (9.15), / is supported in u i, and 
therefore / = Xuf- Then, since m<&~ 1 (x 2 Re(/ 2 )) = Re(m<h~ 1 (,\ 2 / 2 )) = Re(/), 


Pu = Xu Re /, u(0) = u in , 


and we have to choose M so that u(T) = ib. By definition of 4?, recall that w = &u 
means that 

w(t,x) = {l + d x p 1 ^j)~ l (t),x -p(t ) + j3i(% / -1 (f),z -p(t))J 

for t £ [0,Tj], x £ T. Since V’ _ 1 (Ri) = T, we see that u(T) = ib provided that 
w(Ti,x ) = ibM(x) with 

M(z) = {l + 9 x/ 9i(r,x-p(Ti))}3, (9.21) 

and p{T\ ) is given in (C.38). Now the estimates (9.10) and (9.12) follow from (9.17), 
(9.18) and Proposition 5.1. This completes the proof of statement i). 

ii) In what follows, we add the exponent ' to denote the objects associated to 
(V 7 ,d,R 2 ). Let / = &T(ui n ) be defined by (9.20), and let f = & T (ui n ) be the 
corresponding function obtained by taking (V. d , R 2 ) instead of (V, c, R 2 ). We have 
to estimate the difference / — f. If the constant P\ (T) in (9.9) is sufficiently small, 
then u> 2 , X 2 can be chosen so that (9.15) holds both for <!> and for < I >/ . Hence 

f - f = m$ _1 (x2©M,T 1 (^m)) - m , y~ 1 (x 2 ® , M',T[(^' U in))- 
We split this difference into the sum / — f = A\ + ... + Aq, where 
A- L := (m - m')®~ 1 (x2 < 3>M,T 1 (®u in )) 

A 2 := m / $ _1 [x2@M,Ti(^«m - &U in )\ 

A 3 := rn§ 1 [X2(©M,Ti — ©M'.TiX^^m)] 

A 4 := m'^f 1 - '&' 1 - 1 )l/>*<p*[X 2 @M',T 1 (®' u in)] 

A 5 := - ^/,T)[x2©tf',Ti($'»m)] 

Aq := m , $ ,_1 [x2{T©M',Ti( < L / «m) - ©V,r;( $/w m)}] 

and T is the time-rescaling operator defined above, namely ( Th)(t,x) := h(Xt,x), 
with A := T\ /T[. Let us estimate each A 4 . 

Estimate for A\. Apply (C.44). Estimate for A 2 . By construction (see Appendix 
C), ^ _ 1 (0) = 0, p( 0) = 0, and therefore <S>u in = <L| t =o u in = ( , Lr 1 )|t=o(' u m)- Hence 
the estimate for A 2 follows by (C.41) and ( 8 . 6 ). Estimate for A 3 . Apply (C.45). 
Estimate for A 4 . Apply (C.41) and ( 8 . 8 ) with p = 1. Estimate for Aq. Apply (C.42). 
To estimate dtf 2 , use that f 2 solves Pff 2 = 0 (statement ii) of Proposition 8.1), 
and similarly for f 2 . Estimate for Aq. The assumptions (9.9) and (9.11) imply that 
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W,Rz and W',R' 3 satisfy (8.7) with fi = 3/2, which is the hypothesis of statement 
iv) of Proposition 8.1. Then (8.10) holds, namely 

IIT0M',Ti(^«m) - ®'M',T{(^' u in)\\c°{[0,T[y,L 2 ) 

< \\&Ui n \\ H 3/2 (|1 — A| + \\W' — (T1T)|| C 0([ 0iT /].^ 2) + ||i ?3 — (Ti?3)||c°([0,T 1 , ];£(L 2 )))- 

Now ||d)'ttj n ||^3/2 < ||'Uin||#3/2, and the bounds for the last three differences are given 
in (C.44), (C.47) (with cr = 2) and (C.49). Note that assumptions (9.9), (9.11) and 
(9.13) imply (C.46), (C.48), which imply (C.47) and (C.49). □ 

Remark 9.3. The function W contains the terms dtc and V: see Appendix C (see 
also the bound (5.2)). For this reason we assume dfc and dfV to be bounded in 
(9.9) in order to get a bound for dfW, as required by Proposition 8.1. 

3 

Lemma 9.4. If the W 2 '°°-norms of c— 1 and c' — 1 are small enough, then 

ll-^h.s — ^h,s\\c{H^,L 2 ) + ll(Afe, s ) 1 — (A^ s ) 1 11£(7j2 ,h s ) ~ ll c “ ^IliT- 

Proof. By definition (4.1) of A/j s one has 

A/l,S — R' h)S = h b T c (2s)/3_ c /(2s)/3L^ 2b ^ A . 

So the bound for A/j iS — A/ g follows from the paradifferential rule (A.9) and the 
Sobolev embedding H 1 ( T) C L°°(T). To prove the other bound, we use the identity 
(4.3) to obtain that 

Afci - (A'm)- 1 = (/ + h'L% )- 1 [(/ + B)- 1 - (I + B 1 )- 1 ]- 

Recall that ||R||£(^ 2 ) < \ and < \ so the identity 

(I + B)- 1 - (I + B')- 1 = (I + B)~\B' — B)(I + B')~ l 

implies that 

HI + B)- 1 -(1 + B')- l \\c^) < m - B’\\ c{l2) , (9.22) 

and the bound follows from the definition of B , B' and the paradifferential rule (A.9) 
as above. □ 


End of the proof of Proposition 9.1. We recall that Qt is the control operator 
as given by Lemma 9.2, and the operator K, is introduced in (4.15), with ||(I + 
fc)~ l \\c(L 2 ) < 2 if (c, V,i?2) satisfy (9.9), (9.11) and ||c — l||c°([o,T]:H 2 ) is small 
enough. Set 

@s,t[(V, c, R)} := A- 1 s 0 t (J + (9.23) 

and let / := 0 S) t[(R, c, R)](vi n ). Then it follows from the previous construction (see 
Section 4, in particular the Proof of Proposition 3.2 given Proposition f.6) that the 
unique solution v to Pv = T p Xuif, v \t =0 = v in satisfies v(T ) = ib for some constant 
b € M. Since Im f T Vi n (x ) dx = 0 by assumption, by (9.2) we deduce that 


Im / v(T , x) dx = 0. 

J T 

Therefore 6 = 0 and v(T ) = 0. Thus it remains to prove (9.5). Following the 
same argument used in Section 4 to prove (4.16), one proves that WK-Wc^h 3 / 2 ) ^ 1/2, 
whence ||(I + K-)~ l \\c(H 3 / 2 ) — 2- By combining this estimate with (9.12), we have 


C°([0,T];iP 




< 


||© t (/ + /C)- 1 A m u 


m|l C°([0,T];i?i) 


< 


l|Aft, ; 


3 

1 H? 


< 


'H‘ 


•s+ 
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which is (9.5). Finally, we observe that 


ll^2||c°([0,T];£(Z,2)) + H i?2 ll C 0 ([0 ,T]-C(H%)) 

< || C - 1 ,d t C, V 11 C 0 ([o,T] ;// s 0) + H-R|lc°([0,T];£(.f/ s )) + H- R ll C 0([ 0j T];£(H s +f))' ^ 9 ' 24 ) 

This bound for R2 follows easily from the arguments used in the proof of Lemma 4.5 
and Lemma 9.4. Hence, if ( c,V,R ) satisfy (9.4), then (c, V, R2) satisfy (9.9), (9.11). 
This completes the proof of statement i). 

ii ) Given y £ 77 3 / 2 (T), we have to estimate the difference 0 S t[(U, c, R)\(vi n ) — 
0 Si t[O C,d,R!)]{v in ), which is, by definition, 

A hl®T{I + iq^K^Vin - (A; iS )- 1 0 , t (I + K')~ l K' h s Vi n . 

We write it as the sum B\ + ... + B4, with 

Bi := {A^ s - (A^)- 1 }©^/ + Kiy'Ah'sVin, 

B 2 := (A^J-^Qt - &' T )(I + 1Q- l k h , s v in , 

B3 ■= (A^)- 1 ©^/ + /C)- 1 — (/ + JCT^KsVin, 

B 4 := (A^)- 1 ©^/ + K')~\ A hiS - A' h s )vi n . 

If ( c,V,R ) satisfy (9.4), then (c, V, -R 2 ) satisfy (9.9) and (9.11), and ||/C||£(£ 2 ) < 
see (4.16). Then, using Lemma 9.4 and (9.10), we bound the G°([0, T]; 77 s )-norm 
of B\ and B4 by ||c — d\\ H i ||uj n ||tfs. To estimate B2, we want to use (9.14), which 
holds provided that (c, V, R2) and (d,V',R 2 ) satisfy (9.13). One proves that, if 
||c — l||c°([0.T] ; H 3 ) is small enough, then 

11-^2 - ^2llc°([0,T];£(L 2 )) < ll( c ~ c', d t (c — c'), V — V"')|| C -0([ 0 i t] ; H s 0) 

+ II R - R'\\c°([o,T]-,c(H s )) (9.25) 

||^-R2||c°([ 0,T];£(L 2 )) < ||(c— l,dtC,d 2 C,V,dtV)\\ c O([ 0:T ]. HSO ) 

+ || (i?, ^i?)|| c -o([ 0i T];£(i/ S )). (9.26) 

These bounds follow from the arguments used in the proof of Lemma 4.5 and 
Lemma 9.4. Hence assumptions (9.3) and (9.6) imply (9.13), which implies (9.14). 
We have ||(I + tC)- 1 A htS v in \\ H 3/2 < 2||A^^-u^||^ 3 / 2 < Hvmll^+ 3 / 2 . Using (9.25) to 
estimate the last term in (9.14), we deduce that 

11^2 l|c°([0,T];J? s ) < H't’inll^+S {||(c — c', d t (c — d), V — U')|| C - 0 ([ 0)T ]./ ?S0 ) 

+ ||R - R \\c°([o,t]-,c(hs))}- 

It remains to estimate B 3 . The difference 1C — 1C satisfies 

II (A - fc')y\\ L 2 < ||y|l H 3 {||(c- c', V - U')||cO([0,T];i?l) + 11-^2 - ^2 ||c 0 ([0,T];£(L 2 ))} - 

To prove this bound, recall that /C is defined by solving an evolution equation, 
and then, as above, use the energy estimates proved in the appendix to bound the 
difference of two solutions satisfying evolutions equations. Since ||(I + K,)~ l \\c{L 2 ) < 
2, ||(I + 1C) 1 \\c(h 3 / 2 ) — 2) an< i 

(I + 1C )- 1 - (/ + 1C )- 1 = (I + JC)-\IC - K){I + K')-\ 

we deduce that £>3 satisfies the same bound as i? 2 . The proof of Proposition 9.1 is 
complete. □ 
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10. Iterative scheme 


In this section we conclude the proof of Theorem 1.1. As explained in Remark 1.4, 
it is sufficient to prove this result with (jgfinal,'tp final) = (0,0). Also, as explained in 
the introduction, we seek P ex t as the real part of the limit of solutions to approximate 
control problems with variable coefficients. 

Consider the unknown u = T p u — iT q r\ as introduced by Proposition 2.5. As 
proved in §2.3 (see also Section 3), this new unknown u solves an equation of the 
form 

d t u + T V ( u )d x ii + iL? (T c(n )L 2 u)u + R(u)u = T p ^P ext , (10.1) 

where, with a little abuse of notation, we write V(u), c(u), ... as shorthand notations 
for V(r])ip (see (2.3)), c = (1 + (d x r]) 2 )~ 4 , ... where ( 77 , V>) is expressed in terms of u 
by means of Lemma 2.8. 

Fix T > 0. We claim that there is e > 0 such that, for all initial data whose 
iL s (T)-norm (with s large enough) is smaller than e, and all source term P ex t whose 
L 1 ([0, T]; LP(T))-norm is smaller than e, the Cauchy problem for (10.1) has a unique 
solution in C'°([0, Tj; H S (T)). The existence of a solution follows from the analysis 
given below. The uniqueness is obtained by estimating the difference of two solutions 
(as in [ 2 ]) and we omit its proof. 

Recall that T;C) denotes the space of //^-functions whose imaginary part 
have zero mean (see Notation 2.7). 

Proposition 10.1. Let T > 0. For all Ui n € H a { T;C) for some a large enough 
such that ||uin||#<T is small enough, there exists a real-valued function 

P ext eC°([0,T}-,H*(T)), 

such that the unique solution u € C°([0, T]; H a (T)) to (10.1) with initial data Ui n 
satisfies u(T ) = 0. 

Before proving this proposition, let us explain how to deduce Theorem 1.1 from it. 
Recall that it is sufficient to consider the case where ( Ofinal-,^ final ) = (0,0). Once 
P ex t is defined by means of Proposition 10.1 applied with ui n = T Pin Ui n —iT qin rii n , we 
solve the water waves system ( 2 . 1 ) for ( rj, i/j) with data (■ rji n ,ipi n ) with this pressure 
seen as a source term. Then u = T p u — iT q rj solves (10.1), so u(T) = 0 which in 
turn implies that (r/,^)(T) = 0 in view of Lemma 2.8. 

Proof of Proposition 10.1. Set s = cr — 3/2. Given Ui n € H s+ 5(T;C) and T > 
0 , introduce the following scheme: define (uo,/o) := ( 0 , 0 ), and then, for n > 0 , 
(u n+ 1 , f n - i-i) are defined by induction in this way: f n +\ is determined by asking that 
the unique solution u n+ \ to the Cauchy problem 

clfUri+l T Ty( Un )d x Un-{-l T iL 2 T c ( Un ^L 2 U n _|_i T R(u n )u n - 1_1 — /n +1 

Rn+l|t =0 — u in: ( 10 . 2 ) 

satisfies u n +i(T) = 0. 

Our goal is to prove that this scheme converges. Then we define P ex t as the limit 
of (R e f n ) when n goes to +oo. Using the operator Q s ,t defined by Proposition 9.1, 
the scheme corresponds to define (u n ) and (/ n ) as follows: 

fn+l ■= @s,T[X n \(uin) where X n := (V(u n ), c(u n ), R(u n )) (10.3) 

and u n+ 1 is dehned as the unique solution to the Cauchy problem (10.2); by definition 
of f n +i we then have u n+ i(T) = 0. Our goal is to prove that, for any T > 0, if Ui n is 
small enough, then this scheme is well-defined and ( u n , f n ) converges to a solution 
(ti, /) of the desired nonlinear control problem. This will be a consequence of the 
following result. 
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Lemma 10.2. Consider T > 0. There exists sq large enough and for any s > sq + 6 
there exist £o > 0 and positive constants K\,, K 7 such that, for any e € (0, £ 0 ]; */ 

< £ 


\Uj, 


'H s 


■i(T) 


then, for any n > 0 , there holds 


l Un llc'°([0,T];/r + §) — K l£ ’ 




IC°([0,T];H s O) 


< K 2 E for 1 < k < 4, 


Moreover, for any n > 0, 

||^n+l — u n\\c°([0,T];H s ) < K 3 e2~ 


%('U“n +1 Ur 


In 3 ]\^E2 


and for any n > 1, 


l C 0 ([0,T];i? s+ 7) 


< K s e, 


dffn 


C°([0,T],H s o) 


< Kqe for 1 < k < 3, 


||/n+l fn\\c°([0,T];H s ) — ^7 £ 2 


2 «■)—n 


(10.4) 

(10.5) 

( 10 . 6 ) 

(10.7) 

( 10 . 8 ) 

(10.9) 

( 10 . 10 ) 


Proof. For this proof we denote by C various constants depending only on T, s, so 
or u>. Also we denote by T various increasing functions T : M + —>• M+ depending on 
parameters that are considered fixed. 

Step 1: proof of (10.4), (10.5), (10.8) and (10.9). 

We prove these estimates by induction. They hold for n = 0 since (uq, fo) = (0, 0). 
We now assume that they hold at rank n and prove that they hold at rank n + 1. 

We begin by checking that the fact that the properties (10.4)-(10.5) hold at rank 
n implies that one can apply Proposition 9.1 to prove that the scheme is well-defined. 
This means that we have to prove that the smallness assumption (9.4) is satisfied. To 


do so, we first recall that (see (2.4)) HV^ 


\H s o 


< H 


'n II J^so+i 


) 


J n\\H s o+ 1 ' 


the estimate ( 2 . 22 ) (applied with s replaced by So + 1 ) implies that ||P(m« 
||u n ||jy S 0 +i. Similarly, the estimates (2.14) and (2.22) yield 


Then 

\h s o 


||i7K)|| £( ^ + 3 ) <^( 




) 


'H s 


< T 




\U r , 


' H s+ § ’ 


and, directly from the definition c = (1 + ( d x r /) 2 ) 3 / 4 , one has 

||c(tfn) — 111 H s o — ll^nllif^+i ) II ^ II H s o+! ~ ll^nll^so+i • 

Gathering these estimates and recalling that so + 1 < s, we conclude that 

11 V (u n ) 11 + ||c(u n ) — 1 || H s 0 + ||-R(^n)H ^ /tJ s+ 4 \ ~ 11 11 zjs+ 4 • ( 10 . 11 ) 


] C(H S+7 1) ~ 11 n ''H s+a 2 ' 

Consequently, the property (10.4) at rank n implies that the part of the smallness 
condition (9.4) concerning V, c, R is satisfied. Concerning the estimates of the time 
derivatives djfV and d^c, we use the equations (2.1) and the rule (see [27]) 

d t G{rj)ij) = G{rj){d t il)- (B(r))ip)d t r)} - d x ((V(r})ij>)rj) 

(where B(r])ijj and V(rj)ip are given by (2.3)) to express time derivatives d^V and 
d^c in terms of spatial derivatives and in terms of the operators B(r]),V(r]) (see 
Appendix A.3 in [4] or [28, 34]). Then, as above, the desired estimates then follow 
from (2.4) and the usual nonlinear estimates in Sobolev spaces. 

We now prove (10.4) and (10.8) at rank n + 1. By (B.2) we obtain that 


ll u n+l|l CO ([0 )T ] ;i? s+§) — 


u iri\\ H s+% T ||7p„X(ii Re fn+l || 




), ( 10 . 12 ) 
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where the constant C depends on s, T (by (10.11) and (10.4) at rank n, the constant 
M in Proposition B.l is bounded by 1 if K\ is large enough and eo is small enough). 
Now observe that, since T Pn acts on any Sobolev space with operator norm bounded 
by M$(p n ) < FiWitnWns) < 1), one has 

ll^PnXu; R e fn+1 Hf70Qo,T];i? s+ t) — ^ C°([0,T];H s+ ^) ' 

Moreover, by (9.5), ||/ n +i || c0 < K 0 \\u in \\ HS+ 3 for some K 0 depending 

only on T. We conclude that, choosing K\ large enough and eo small enough, (10.4) 
holds at rank n + 1. Also (10.8) at rank n + 1 follows by the same argument. 

It remains to prove (10.5) and (10.9). Directly from the equation (10.2), expressing 
dtu n+ 1 in terms of u n ,u n+ i and f n + i and using the operator norm estimate (A.9) 
for paradifferential operators, one deduces (10.5) for k = 1 from the bounds (10.4) 
and (10.8). We next prove (10.9) for k = 1. To do so, the key point is to make 
explicit the equation satisfied by f n +\ ■ We recall from (10.3), (9.23) and (9.20) that 

fn +1 := (AS lfl )- 1 K(^)- 1 ( X2 /n+l)), fn+1 ■= ©M n ,T P (*"(/ + A^)” 1 AJ^) , 

where Af s , <L n , rri n , M n , /C n , T'{ L are given by replacing (V,c) with (V(u n ), c(u n )) in 
the definition of A^ s , 4>, m, M, 1C, T\. By definition of ®m,t (Lemma 8.2) one has 


dt.fn+i T ('U>n)dxfn +1 T iLfn- 1-1 T R^f^n)fn+1 — 0, 

fn+l\t=T™ fn+1’ 


(10.13) 


where R^(u n ) is given by (9.8) and the initial data /^ +1 is given by Lemma 8.2. It 
follows from ( 8 . 8 ) that 


H/rH- 1 !! 




3 <KU n (I + JC n )~ 1 Al s u in \ 


3 

H? 


< K II Ur 


'H s+ % 


Using the equation (10.13) we thus estimate the (7 0 ([0, T]\ L 2 )-norm of dtf n +i from 
which we estimate dtf n +1 in C°([0, T]; H s ). This gives (10.9) for k = 1 since s > so- 
Now we obtain (10.5) for k = 2,3,4 as well as (10.9) for k = 2,3 by differentiating 
in time the equations satisfied by u n+ \ and / n +i- 


Step 2 : proof of ( 10 . 6 ), (10.7), ( 10 . 10 ). 

The estimate (10.10) will be deduced from (10.6) and (10.7). To prove (10.6) and 
(10.7) we proceed by induction. We assume that they hold at rank n — 1 and prove 
that they hold at rank n. 

The key point is to estimate 5 n := u n+ \ — u n . Write 

dt$n + T v(un) d x 6 n + iL 2 T c ( Uri )L 2 S n + R[u n )5 n = G n (10.14) 

with 

G n •= (Tv(u n _i) — Tv(un))dx'U'n + iE 2 (T c ^ Unl ^ — T c ( Un j)Liu n (10.15) 

+ (R{u n — i) R(u n ))u n + Tp( Un ')Xoj(fn +1 fn) + (^p(u„) ^p(ti„_i))Xu/n' 

As in the previous step, it follows from Proposition B.l (noticing that <5 n +i(0) = 0) 
that ||5 n ||c 0 ([o,T];iP) < Co\\Gn\\ C o([o, T ]-, H s) for some C 0 depending on s,T. 

Estimate for G n . We claim that 

ll < ^"-llc°([0,T];P s ) — £ K(T) ||<5 n —l|| c o ([ 0 ,T]- iH s)+ £ K(T) ||5t5„_i||^ 0 ^ o ^^ s _3^ . (10.16) 
Let us prove this claim. At each t £ [0,T], using (A.9) one has 

||(^V(u„_l) ~ -^V (u„) )^ X Un 11 Jfs ~ ||U('li n _l) — V(u n ) ll^oo ||5j;tl n ||jj s . 
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It follows from (10.4) that \\d x u n \\ HS < K\£. To estimate V{u n -\) — V{u n ) we use 
the following consequence of Lemma 5.3 in [3]: Assume s > 3/2 and consider ( 771 , 772 ) 
such that \\rn\\ Ha + ||» 72 II< 1- Then 

II - G( m )f 2 \\ HS _s < K \\ m - m \\ RS _r ll/rll^ + K ||/r - f 2 \ |^_i . 
Then, directly from the definition of U0)0 one deduces that 

ll^(’7i)'0i - Vim^Wni < K 11 771 - 7721^2 |0i|05/ 2 + K H-01 - -02102 • 

Since H 1 ( T) C L°°(T), we then conclude that 

10 071— l) V 077 ) II £00 0 1071 Vn—1 |0S T 11 -071 "071 — 1 |0s 0 |071 ^ 77—1 10s ■ 

The estimate of the H s norm of L^(T c ^ Un _ 1 ^ — T c / Un ^)L^u n is similar. To estimate 
(R(u n -i) — R(u n ))u n recall that R(u)u is as given by Proposition 2.5. This operator 
is defined by means of the remainder E0)0 in (2.7) and also in terms of explicit ex¬ 
pressions involving symbolic calculus or the paralinearization of products. The only 
delicate point is to estimate i ? (?y n )0 n — F To do so one uses Lemma 6.8 

in [ 2 ]. 

It remains to estimate the last two terms in the right-hand side of (10.15). Directly 
from (A.3) we find that 

IK T p(n„) - T p{u n -i))Xufn\\ HS < M$(p(u n ) - p0„_i)) ||Xw./n||jf» ■ 

Now ||xa;/7i|0s < ||x008 11010s < £ by (10.8), and M$(p(u n )-p(u n -i)) is bounded 
by K | \u n - iiji — 1 1| H s • Eventually, to estimate the iL s -norm of T p ( Uri )Xoj(/«+i - In) 
we use again (A.3) to bound this expression in terms of \\f n +i — / n |0 s . We use (9.7) 
to obtain 

II/tH-1 — fn\\c°([0,T];H s ) 

I077i||^-s+| {ll( c 7i — C n —li dt(c n — C n —i),V n — V n — 1 ) |0°([O,T];7/ S o) 

+ | 07 i - -Rn-l||c°([0,T];£(H s ))} 

10771 II { 1077 — u n—l I 0 s o 4" 11 <0 — U n — l) | 0 s o } > (10.17) 

and then we use (10.6) and (10.7) at rank n — 1. 

Estimate for u n+ \ — u n . For eqK(T)Cq < 1/2, it follows from (10.6) at rank n — 1 
and (10.16) that the desired result ( 10 . 6 ) at rank n holds. 

Estimate for f n+ 1 — f n . The estimate (10.10) follows from (10.17) and the assump¬ 
tions (10.6)-(10.7) at rank n — 1. 

Estimate for dt{u n+ \ — u n ). By (10.14), 

d t 5 n = -T V ( Un) d x 5 n - iL^T c{;an) L^6 n - R(u n )S n + G n . (10.18) 
As above, one has 

\\ T V(u„)\\c{H s ,H s - 1 ) + \\ L2 T c(u„)L 2 \\ c( ^ Hb hS —3 + 10071)10(77+77®) < C'|0ti| 0 s - 

Therefore one can use (10.6) and (10.4) to estimate the first three terms in the right- 
hand side of (10.18). The last term G n is estimated by means of (10.16) and the 
induction assumptions. Consequently, we get ||di<5 n ||^ 0 ^ o — Ce 2 2~ n , and 

for e < £q, with eo small enough, we deduce (10.7). □ 

We can now conclude the proof of Proposition 10.1. 

Recall that s = a — 3/2 by notation. By (10.6) and (10.10), we deduce that 
07 i) 77 gN and (f n )n£N are Cauchy sequences in C°([0, T]; H 8 ) and therefore converge 
to some limits u and / in C°([0, T]; H s ). Using the uniform bounds (10.4) and (10.8) 
and the interpolation inequality in Sobolev spaces, we infer that 0 n ) n£ N and (f n ) n eN 
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converge in C'°([0. T]: H s ' + 2 ) for all s' < s. Furthermore, we get that u and / belong 
to C°([0,T]-, H s + %) fl L°°([0, T , ];i7 s+ i) for all s' < s. Passing to the limit in (10.2), 
we conclude that u and / satisfy (10.1) and u(T ) = 0. Eventually, using Lemma B.l 
(seeing (10.1) as a linear equation of the type (B.l) with unknown u and coefficients 
in L°°([0,T];FP)), we deduce u € C°([0, T]; H s+ l). 

It remains to prove that / € C'°([0, T]; H s+ z ). We know that u n —> u in 
C°([0,T]-, H s ) C C 0 ([0, T]\H So+6 ). As a consequence, V(u n ) —X H(w), c(u n ) —> 
c(u), dtc(u n ) —»■ dtc(u), p(u n ) —>■ p(u) in C°([0, T]; FZ" 80 ), and R(u n ) —»■ iZ(w) in 
(7 0 ([0, T]; £(H S )). Now consider /oo := 0 s ,t[V(u),c(w), iZ(w)](«j re ), and recall the 
dehnition (10.3). By (9.7), ||/„ — /oo||c°([ 0 ,T];i 7 s ) -1 0 as n —>• 00 . On the other hand, 
/ = lim/ n in C ,o ([0, T]; FF S ), and therefore / = /oo. By statement (i) of Proposition 
9.1, /oo € C°([0,T]; FP + 2 ), with estimate (9.5). 

This concludes the proof of Proposition 10.1 and hence the proof of Theorem 1.1. 

□ 


Appendix A. Paradifferential operators 


Notation A.l. For p € N, we denote by W P, °°(T ) the Sobolev spaces of L°° func¬ 
tions whose derivatives of order p are in L°°. For p e] 0 , + 00 [\N, we denote by 
W p ’°° (T) the space functions in VF^’°°(T) whose derivatives of order [p] are uni¬ 
formly Holder continuous with exponent p — [p\. 


Definition A.2. Given real numbers p > 0 and m € R, F” 1 denotes the space of 
functions a(x ,£) onTxR which are C°° with respect to £, and such that, for all 
a € N and all £, the function x e-x d^a(x,f,) belongs to W P, °°(T ) and 

Definition A.3. For m € R, p € [0,1] and a € r™(R d ) ; we set 


M?{a) 


sup sup 
M<6+p ?6K 


(l + |ei) H “ m ^a(-,0 


WP’°°( T) ’ 


(A.l) 


Now consider a C°° function x homogeneous of degree 0 and satisfying, for 0 < 
£1 < £2 small enough, 


X{0, rj) = 1 if \6\ < £1 17/|, x(0, V) = 0 if \°\ > £2 M ■ 


Given a symbol a, we define the paradifferential operator T a by 

T a u(f,) = (2vr) _1 x(£ - V, V)a{£ ~ V, v)u{v), (A.2) 

jj ez 

where a(0,£) = J e~ lx ' e a(x, f) dx is the Fourier transform of a with respect to the 
first variable. 

The main features of symbolic calculus for paradifferential operators are given by 
the following theorem. 


Definition A. 4. Let m in R. An operator T is said of order m if, for any /i£R, 
it is bounded from H p { T) to H p ~ m (T). 

Theorem A.5. Let m € R. 

(i) If a € T™, then T a is of order m. Moreover, for any /i£R there exists K > 0 
such that 

II T 

|| - L a 
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(A.3) 






(ii) Let ( m,m') E M 2 and p E (0, +oo). If a E T ™,6 € T™* i/ien T a T b — T a $ is of 
order m + ml — p where 


» b =E 7 


ano = 


M<P 




Furthermore, for any p E M t/iere exists K > 0 such that 

\\T a T b - T a (j fe || £( . HMiJ/Al _ m _ m / +p ^ < KM™(a)M™ ( b ). 

In particular, if p E (0,1], a E T™, 6 E T ™ 2 f/ien 

||r a T fe — T ab || £( .^ M < KM™(a)Mp 1 (b). 


(A.4) 


(A.5) 


(A. 6 ) 


(m) Lei m E 1, /) > 0 and a E r" l (M d ). Denote by (T a )* the adjoint operator of T a 
and by a the complex-conjugated of a. Then (T a )* — T a * is of order m — p where 


a = 


£ ra*** 


M<P 

Moreover, for all p there exists a constant K such that 


II(Ta)* - T a - II H ^H,- m+P < KM™(a). (A.7) 

Remark A.6. These properties are well-known when Sobolev spaces of periodic 
functions are replaced by Sobolev spaces on the real line. To prove these results for 
periodic functions, one can use the results proved in [3] about the general case of 
uniformly local Sobolev spaces Namely, in [3], the above results are proved 

to hold when H S (T) is replaced by In particular it is proved that 

\\{T a T b - T ab )u\\-.»- m - m ' +p < KM™(a)M™\b ) \\u\\ H , . 

n ul r a ut 

Since ||tt||^- s < ||n|| ffS ( T ), it follows that 

II (T a T b - T ab )u\\ K - m - m . +P(m < KM™(a)Mf(b) ||n||^ (T) . 

Now, if u is a periodic function and a and b are periodic in x, so is ( T a T b — T ab )u 
and we deduce that 


I( T a T b T ab ')u\\ji^— rn — rn i+p(Y') ^ || {T a T b T ab )u\\+p, 


By combining the previous estimates we obtain (A. 6 ). The other estimates are 
proved in a similar way. 


It follows from (A. 6 ) applied with p = 1 that, if a E b E Tf 1 ' then 

IIM] KM™(a)M™' (b). (A. 8 ) 

If a = a(x) is a function of x only, then T a is called a paraproduct. We often 
use that the following consequence of (A.3): if a E L°°(T) then T a is an operator of 
order 0 , together with the estimate 


Vcr E M, \\T a u\\ H a < \\a\\ LOO \\u\\ H a ■ (A.9) 

If a = a(x) and b = b(x) then (A.4) simplifies to afb = ab and hence (A.5) implies 
that, for any p > 0 , 


T a T b | Ip I'm* ff/j-m-m'+p) ^ 


\C(H^,H^- rn - rn 
provided that a and b are in to W p,co (T) 


< K ||a| 


\\ U \\WP’ C ° > 


(A.10) 
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Theorem A.7. i) Given two functions a,b defined on M we define the remainder 

lZ(a,u) = au — T a u — T u a. (A.11) 

Let a £ M + and (3 £ R be such that a + fi > 0. Then 

\\K(a,u)\\ Ha+/3 _i < I< \\a\\ Ha \\u\\ H p . (A.12) 

ii) Let a > 1/2. For all C°° function F with F( 0) = 0, if a £ H a ( T) then 

||F(a) - T F , {a) a\\ H2a _i < C (||a||^ s ) ||o|| Ha . (A.13) 

Proposition A. 8 . Let r, /i £ R be such that r + fi > 0. // 7 € M satisfies 

7 < r and 7 < r + p — 

then there exists a constant K such that, for all a £ Ff r ( T) and all u £ T), 

11ait — T a u\\ Hl < K ||a||jj r ||u||# M . (A.14) 

We also recall two well-known nonlinear properties. Firstly, If u\, U 2 £ H s ( T) n 
L°°(T) and s > 0 then 

Iliti^H^s < K ||iti|| LOO + K ||«2|| i0 o ||iti|| HS , (A.15) 

and hence, for s > 1/2, 

||'ui'u 2 ||# s < K HitiH^s \\u2\\ H s ■ (A.16) 

Similarly, for s > 0 and F £ C°°(C N ) such that F(0) = 0, there exists a non¬ 
decreasing function C : M+ —> M + such that 

\\F(U)\\ HS <C{\\U\\ Loa )\\U\\ HS , (A.17) 

for any U £ (H S (T) n L°°(T)) JV . 

Appendix B. Energy estimates and well-posedness of some linear 

EQUATIONS 

Recall the linearized equation dpu + iLu = 0, where L := fig — df)\D x \') A 
We gather in this section Sobolev energy estimates for linear equations of the 
form 

dtip + Vd x ep + iLa ( cL^ep ) + Rep = F, 

where V = V(t,x) is a real-valued coefficient, c = c(t,x ) is a real-valued coefficient 
bounded from below by 1/2, F = F(t,x) is a given complex-valued source term 
and R is a time dependent operator of order 0 which means that Rep is defined by 
( Rep)(t ) = R(t)ep(t) and R belongs to (7 0 (M + ; £(77^)) (for some fi) where C(Hfi de¬ 
notes the set of bounded operator on H^{ T). Below we consider various equations of 
this form where, for instance, R is either a multiplication operator by some function 
or the commutator between Vd x and a Fourier multiplier. 

We also consider paradifferential equations of the form 

d t ep + Tyd x ep + iCep + Rep = F, 

where C = La (T c Lz ■ ) and V, c, R are as above. 

Proposition B. 1. Let T > 0 and p £ [0, + 00 ). Consider R £ C°([0,T]; C(F[^)) 
and real-valued coefficients V, c satisfying 

V £ C'°([0,T];IF 1 ’ oo (T)), c £ C°([0, T]; W*'°°(T)), 

with the Lff-norm of c— 1 small enough. 
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For any ipi n € H fi ( T) and any F € L 1 ([0, T]; H^iT)), there exists a unique <p € 
C°([0,T]-H^(T)) such that 

d t ip + T v d x ip +R(p + i£(p = F, P\t=o = Pin- (B.l) 

Moreover, for any t > 0, 

— eCt (ll^mll.ffM + II-^IIl 1 ^,*];/*>)) > (B-2) 

for some constant C = C(p,M ) depending only on p and 

M = Jot] { ^ dxV ^ L °° + + II^WII c{m) }■ 


Remark B.2. We often use energy estimates for backward Cauchy problems, that 
is for Cauchy problems on time intervals [0, T] with a data prescribed at time T. 
Then the energy estimates read 

11^)11/^ ^ eC (ll¥ ? (^~')ll+ II^ ? IIl 1 ([0,T];H'')) ' (B-3) 

Proof. As already seen in (2.21), C = L^T C L^ = T 7 + R' where R! is of order 0 and 

7 = cl+ \(d^y/J)Vld x c. 

Up to replacing in (B.l) the remainder R by R + iR', we prove the existence of the 
solution as limits of approximate problems of the form 


d t p> + Tyd x J £ p> + iT~fJ e ip + Rip = F, <£>| t=0 = JePin , (B.4) 

where J £ are smoothing operators. Then (B.4) is an ODE in Banach spaces and 
admits a global in time solution denoted by ip £ . 

Set (f, x, £) = c(t,x)d.{tf), which is the principal symbol of 7 . As in [ 2 ], 

consider the paradifferential operator A^ with symbol 1 + (c(t, x)£(£)) 2/i//3 and, given 
e € [0,1], define J £ as the paradifferential operator with symbol j £ = j e (t,x,£,) given 
by 

3 c = = exp ( - S 7 < 3 ' 2 >) - \(B x ap exp ( - <r 7 W 2 >). 

Recall that the Poisson bracket of two symbols is {a, b} = (d x a)(8^b) — (d^a)(d x b). 
Then 

{i 0) ,7 (3/2) } = 0, {i»),(rf)W 3 } = 0,{ 7 W 2 >,(rf)« 3 }=0, (B.5) 

and 

Imji -1) = ~{d x d^)jf\ 

Of course, for any e > 0, j £ € C°([0,Tj;T^^M" 1 )) for all m < 0, so that T Je u € 
C°([0, T]; H°°(T)) for any u € C°([0, T]; H~°°(T)). Also j £ is uniformly bounded in 
C°([0, Tj; r^ 9 (M d )) for all e € [0,1]. Hence, using (A.5) with p = 3/2 or (A. 6 ) with 
p = 1 , we have the following estimates (uniformly in e): 


II \Jsi ^ c ||w|| , 


II (Je 


*u — J £ u\\ 3 < C I 


U 


Ha- i 


A M ,Li(T c Li) 


u 


L 2 


— C IMI_H> ) 


II [A^, Je] ’ u 'll^§ — C |M|/f> , 


|| [A^, T v d x j e ] U \\ L 2 < C II V II yyl, 00 |M|#> 1 || [Je: < C || U || ^- 1,00 , 

(B.6) 

for some constant C depending only on ||c|| ^ 3 / 2,00 and uniform in e € [0,1]. 

Recall that, by notation, ip £ is the unique solution to (B.4) and introduce <p e := 
A tl p> £ . Using the fact that A^ is invertible (for c — 1 small enough) and the preceding 
estimates, we deduce that 


dtp’s T r J'v'd x J £ dp £ + A^RA^ ip £ + iTryJ £ f> £ — F £ , 
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Ps\t =0 — A^J £ (pi n , (B.7) 








where 


ll-^e IIL 1 ([0,T];Z, 2 ) ^ C(M){ II < / , eIIl 1 ([0 ,T]-Hr) + II^11 V-([0,T\-,Hn) }' ( B ' 8 ) 

Write i mh = 2 Re (dt^pe, <p £ ), where (•,•} denotes the scalar product in L 2 (T), 
and hence 

^\\Ve\\ 2 L 2 =-{{P+ P*)Ve,<Pe) with 
P = Tyd x J e + A^RA ^ 1 + iT 7 J £ . 

To estimate the operator norm of P + P*, there are two ingredients. Firstly, we 
replace J* by J e + ( J* — J e ) and commute J e with Tyd x and T 7 . This produces 
remainder terms that are estimated by means of (B. 6 ). The proof is then reduced 
to the case without J £ and it suffices to estimate the operator norm of P + P* where 
P = Tyd x + A^RAj 1 + iTj. Since A^PA ” 1 is bounded from L°°([0, T]\ L 2 { T)) into 
itself with an operator norm estimated by M, it remains only to estimate Tyd x + 
iTry + ( Tyd x + iP,)*, which can be done directly by means of the paradifferential 
rule (A.7). We conclude that 

J t Hollis < C(M) ll^llia + \(2F e ,<p e )\. (B.9) 

We thus obtain a uniform estimate for the P°°([0, P]; L 2 )-norm of tp e (from Gron- 
wall’s inequality and (B. 8 )) which gives a uniform estimate for the L°°([0, T]; iP 4 )- 
norm of <p £ . From this uniform estimate and classical arguments (see [35]), one 
deduce the existence of a solution in L°°([0,T]; P M (T)). The uniqueness is obtained 
by considering the equation satisfied by the difference of two solutions and perform¬ 
ing an L 2 -energy inequality (using similar arguments to those used above). The 
continuity in time of the solution is proved as in [2, §6.4]. □ 

Lemma B.3. Consider real-valued coefficients V, c satisfying 

V € C°([0, T]; W 1,00 (T)), c € C°([0, T]; W?'°°(T)), 

with the Lff-norm of c— 1 small enough. Consider also R € C 0 ([0, T]; C(L 2 )). 

i) For any ipi n € L 2 (T) and any F € L 1 ([0, T]; L 2 (T)), there exists a unique 
ip € C°([0, T]; L 2 (T)) such that 

d t y + Vd x ip + Rp + iL^ (cL^ip) = F, <p \ t=0 = ip in . (B.10) 

Moreover, for any t > 0, 

M*)llz,2 < eX P M(t') dt'^j (\\<Pin\\ L 2 + ||P|Il1 ([0,t];L2)) > (B- 11 ) 

with M{t') = \\d x V(t')\\ LOO + \\R{t')\\ c{L 2 y 

ii) Let Pi G [0,3/2], Assume that V € C°([0, T]; H 2 (T)), c € C°([0, T]; P 3 (T)) 
and R G C°([0, T]; C{H»)). If ip in € 7P 4 (T) and F € P 1 ([0, T]; iP l (T)), then, for 
any t > 0 , 

\\p(t)\\ H » < exp (J M(t') dt'^j (\\win\\ H y- + ||P UpdO,*];^) ) , (B.12) 

with M(t') = \\V(t ')\\ H 2 + ||c(f )||^3 + \\R(t')\\ c{my 

Remark B.4. Consider a backward Cauchy problem, that is a Cauchy problem 
with a data prescribed at time T. Then (B.ll) implies that 

M*)IIl 2 < exp (£ M{t')dA (\\<p(T )\\ L2 + ||P|| L 1 ([ 0 ,T];L 2 )) , (B.13) 

with M{t') = \\d x V{t')\\ Lao + ||P(t')ll£(L 2 )- 
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Proof, i ) The existence of the solution can be deduced from the previous proposition, 
writing 

dtp + Vd x p + Rip + iLz [cL^ip) 

under the form 

d t p + T v d x p + R'p + iL? (T c L 2 p^j 

where 

R'p = Rp + ( Vd x p — Tyd x p) + 1 L 2 ((c — T c )L 2 p ). (B.14) 

Indeed, R' belongs to C°([0,T]-, C(L 2 )) in view of (A.12) and (A.14). 

In order to see that the energy estimate does not depend on the norm of c, start 
from i Mb = 2 Re (dtp, p) ■ Since R e(iL 2 (cL 2 p) } p'j = 0, we obtain that 

^ MlL = 2R e(-Vd x p- Rp- F,p) . 

Hence, integrating by parts, 

j t \\p\\l 2 = Re (((d x V) - 2 R)p - 2 F, p) , (B.15) 

and the result easily follows from Gronwall’s inequality. 

ii) This follows from (B.2) and the fact that the remainder R! in (B.14) belongs 
to C°([0,T]-, C(H f1 )) in view of (A.12) and (A.14). □ 

Appendix C. Changes of variables 


Recall the operator 

P:=dt + Vd x + iL*(cL*-)+R 2 , (C.l) 

where L := (g — d xx )zG(0) 2 , the operator R 2 is of order zero, and c(t,x), V(t,x) 
are real-valued functions. Consider a time-depending change of the space variable 
(namely a diffeomorphism of T) and its inverse, 

x = y + Pi(t,y) O y = x + /3i(t,x), 

x,y E T, t E M, with ||||, 1111z,°° < 1/2- Introduce a self-adjoint variant of 
the pull-back operators, defined by 

(^1 h)(t,y) := (\ + d y Px(t,y))^h(t,y + P\(t,y)), (C.2) 

(’Fj/ 1 h)(t,a;) := (l + d x Pi(t,x))^h(t,x +Pi(t,x)), (C.3) 

and note that are self-adjoint with respect to the standard L 2 ( T) scalar 

product in space, for any t. We want to compute when Qo is a Fourier 

multiplier (the analysis below applies more generally assuming only that Q 0 is a 
pseudo-differential operator). 

C.l. Change of variable as a flow map. Introduce a parameter r E [0,1] and 
consider a diffeomorphism of T (depending on (r, t)) and its inverse, 

x = y + P(r,t,y) y = x +P(r,t,x), 

x,y € T, r E [0,1], f E I, where /3 and /3 are such that ||9y/S||x,oo, ||9 x /3||ioo < 1/2 
and 

/3| r= o = 0, (3\ t =o = 0, /3| t =i=/3i, /3| t =i = /3i- 

We denote 

(T(r)/i)(f,y) := (1 + d y j3(r, t, y))^h(t, y + /3(r, t, y))\ (C.4) 

(T(r)' 1 h)(t,.r) := (1 + d x /3(T,t,x))2h(t,x + /3(r,t,x)). (C.5) 
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Then T] = 'L(l). The reason to introduce the parameter r is that 'F(r) satisfies an 
equation of the form 

d r T(r) = F(r)T(r), T(0) = /, (C.6) 

namely <9 r (T(r)/i) = F(r)('I'(r)/i), T(0)/r = h for all h, where 


F(j) = b 0 {T,t,y)d y + -{d y b 0 )(r,t,y), b 0 (T,t,y) 


d T P(r,t,y ) 

1 + d y P(T,t,y) 


(C.7) 


Assume that Q o is a Fourier multiplier with symbol qo{Q of order m < 3/2. We 
seek a pseudo-differential operator Q(t) of order m such that the difference 


R(r) := Q(r)tf (r) - T(r)Q 0 (C.8) 

is an operator of order 0. Commuting Q(r) with the equation = i ? (r)'k(r) 

one obtains 


a T (g(r)T(r)) = Q(r)F(r)T(r) + (5 r g(r))T(r) 

= F{t)Q{t)^{t) + ([Q{t),F(t)]^(t) + (d T Q{r ))'k(r)). 

On the other hand d T {^>{r)Qo) = F(t)^{t)Qq. By combining both equations we 
obtain that R satisfies 


d T R{r) = F(t)R(t) + 77i(r)^(r), 77i(r) := [Q(r), F(t)] + <9 t Q(t). (C.9) 

The analysis is then in two steps. The main step consists in proving that Q(t) can 
be so chosen that Q(t = 0) = Qo (then -R(O) = 0) and 1Z\ (r) is of order 0. Then, 
by using an L 2 -energy estimate for the hyperbolic equation d T u = b^dyU + /, one 
deduces an estimate for the operator norm of R{t ) uniform in r (and hence the 
desired estimate for r = 1). Here we describe in details only the main step, as the 
L 2 -energy estimate is a standard argument. 


C.2. Expansion of the symbol. Let p(r, t, x, £) be the symbol of Q(r). To obtain 
1Z\ of order zero amounts to seek p such that d T p — has order zero (where 

&[f,q] is the symbol of [F, Q]), and p\ T =o = Qo- The asymptotic expansion of 

is 

oo 1 

a [F,Q] ~ E {(3rm a p) - (c.io) 

a=l 

where f(r,t,x , £) := ibo(r,t,x )£ + ^{d x bo){ T ,t,x) is the symbol of F(t) (we rename 
x the space variable). Since m < 3/2 < 2 by assumption, it is enough to determine 
the principal and the sub-principal symbols of p. Thus we write p = po + pi, where 
Po has order m and p\ has order m — 1. The equations for po,Pi are 

d T Po = b 0 d x p 0 - £(<9x& 0 )%>0, Po|r=o = qo, (C.ll) 

d T Pi = b 0 d x pi - £(<9 x 6 0 )%h + z, pi| T=0 = 0, (C.12) 

where 

^ := ^(dxxb 0 )(d^p 0 + £%p 0 )- (C.13) 

If po,pi satisfy (C.ll),(C.12), then it follows from standard symbolic calculus for 
pseudo-differential operators (similar to (A. 5)) that TZ\(t), defined in (C.9), is an 
operator of order 0 satisfying 

\\Ki(t)\\ c{l2} + \\Fi(t)\\ c{hI) < (M r m (po) + M™~\pi)) ||6o(r)||^,cc (C.14) 

with r large enough (here the semi-norms are as defined by (A.l); one has to 
consider r large enough because we are here considering pseudo-differential operators 
instead of paradifferential ones). 


63 




Equation (C.ll) can be solved by the characteristics method: if x(t),£(t) solve 
A x ( r ) = -b 0 (r,t,x(T)), J^£(t) = £(T)(d x b 0 )(T,t,x(T)), (C.15) 

then 

Po(r,t,x(r),^(r)) =p o (O,t,x(0),£(0)) Vr. (C.16) 

Now, by (C.7), the first equation in (C.15) is 

0 = {1 + (d x fj){T,t,x{r))} x'{ t) + (. d T /3)(T,t,x(r )) = ^{x(r) + P{r,t,x(r))}, 
whence 

x(t) + P(r,t,x(T )) = x(0) + (3(0,t,x(0)) = x(0). (C.17) 

Applying the inverse diffeomorphism, we get x(r) = x(0) + /3(r, t, x(0)). This is 
the solution x{r) of the first equation in (C.15) with initial datum x(0). Also, one 
verifies that 

£( r ) = £(0)(l + (9*/9)(r,t,x(r))) (C.18) 

satisfies the second equation in (C.15), because x(t) satisfies the first equation in 
(C.15), bo is given by (C.7), and 


d x b 0 (T,t,x) = 


d rx f3{T, t, x) d r (3(T, t, x)d xx p(r, t, x) 


1 + d x P(r,t,x) [1 + d x /3(r, t, x)] 2 

Hence we deduce a formula for the backward flow of (C.15): fixed any n E [0,1], and 
given any (aq, £i), the solution (x(t),£(t)) of (C.15) with initial datum (x(0), £(0)) = 
(xo,^o) satisfies (x(ri), £(ti)) = (aq,£i) if the initial datum is 

C 


X 0 = Xl+P(Tl,t,Xl), £,0 = —^ n , 

l + d x p(n,t,xi) 

As a consequence, using (C.16) and the initial datum in (C.ll), we get 

Po( r i,f,aq,£i) =po(0,t,xo,fo) = qo{t,x 0 ,£ 0 ) 

= qo(t,xi + P{ri,t,x{),- -—-- 

' 1 + d x /3(T 1 ,t,xi) 

We have a formula for the solution po(r,t, x, of (C.ll): 

Po{r, t, x, £) = q 0 (t, x + P(t, t, x), -- —S -- 

^ 1 + d x /8(T,t,x) 

Now we study equation (C.12). By the definition of (x(r),£(r)), 


(C.19) 


Pi(r,t,x(r),^(r)) = / z(s,t,x(s),£(s))ds, 
Jo 


(C.20) 


(C.21) 


where 2 is given in (C.13). We examine 2 in detail. By (C.20), for k = 1,2, 


dj;Po(T,t,x,€) = 


t,x + (3(r,t,x), ■ 


£ 


1 


l + d x /3(r,t,x)J [1 + d x $(r,t,x)] k 
for all r,t,x,£. Hence along the curves (x(s),£(s)), by (C.17),(C.18), one has 

fak w , , , „ (dcqo)(t,x o,£o) 

{9fPo)(s,t,x(s),^(s)) = --—=--, 

? [1 + d x /3(s,t,x(s))] k 

where (xo,£o) := (a;(0), £(0)), and therefore, using (C.18) again, 

%lo(C x 0 , g 0 ) + £o%gp(t, x 0 , £q) 


(%>o + £%Po)(s,i,x(s),^(s)) = 


l + d x /3(s,t,x(s)) 
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Now we note that 

(d xx bo)(s,t,x(s),£(s)) _ d f (d xx j3){s,t,x(s)) \ 

1 + (d x j3)(s,t,x(s)) ds \ [l + (d x j3)(s,t,x(s))] 2 J ’ 

as it can be verified by a straightforward calculation, using also (C.15) and the 
definition (C.7) of 6 q- Hence, recalling the definition (C.13) of z, 


z(s,t,x{s),£{s)) 
and, by (C.21), 


1 

2 


{d^q 0 (t, x 0 ,Co) + €od&qo(t, xo, £o)} 


d_ 

ds 


(d xx /3)(s,t,x(s)) \ 

[l + (d*0)(s,i,a;(s))] 2 J ’ 


Pl(T,t,x(T)^(T)) 


1 

2 


{d^qo(t, x 0 , Co) + Cod^qoit, x 0 , £ 0 )} 


(d xx p)(T,t, x(t )) 

[1 + (d x p)(T,t,x(r))] 2 


because (3 | r= o = 0. We use the backward flow as above: given ti,xi,£i, the 
solution (x(t),£(t)) of (C.15) with initial datum (,t(0),£(0)) = (xo,Co) satisfies 
(x(ri), £(ti)) = (xi,^i) if the initial datum is (C.19). Therefore, replacing (xo>£o) 
by (C.19) in the last equality, we get a formula for p\, which, writing t,x,C instead 
of is 


Pi(r,t,x,C) = ^{(d£q 0 )(t,x + P(r,t,x), 

c 


+ 


1 + d x /3(r,t,x) 


(d^qo)[t,x + /3(r,t,x), 


1 + d x /3(T,t,x) 

C 


(C.22) 


d xx P(r,t,x) 


1 + d x /3(T,t,x)' > (1 + d x f3(r,t,x)) 2 


C.3. Conjugation of L. We fix qo(C) to be the symbol of L (see (2.11)) with a 
cut-off around £ = 0, namely 

?o(0 := (g + C 2 )^(C)^x(C) = (g + C 2 )^- 1^1^ tanh^(fo |CI)x(0, 

where x(C) is the cut-off function of Proposition 2.5. Note that Op(go) = L on the 
periodic functions, as their symbols coincide at any ^£Z, and therefore no remain¬ 
der is produced replacing L by Op(go)- In the previous section we have constructed 
Po,Pi, and we have defined p := po + pi, Q(t) := Op (p). Then TZ\(r) defined in 
(C.9) is an operator of order zero and it satisfies estimate (C.14). Now observe, 
in view of (C.9), that for any function u o € L 2 (T), R{r)u o solves an hyperbolic 
evolution equation. Using the energy estimate (B.ll), we deduce that the difference 
R{t) := Q(r)T(r) — 4 ’(r)L (see (C.8)) is also of order zero, and it satisfies the same 
estimate (C.14) as TZi (r). As a consequence, the conjugate of L is 

T(r)TT(r) _1 = Q(r) + 77 2 (r), 77 2 (t) := (r)T(r)- 1 (C.23) 

and 77 2 (t) satisfies the same estimate (C.14) as R(t). By formula (C.20), po = 
go(£(l + dxfi)^ 1 )- We expand 

po = (1 + d x j3)~%q 0 + r, (C.24) 

where the remainder r satisfies || Op(r)|| £ ^ M h ^+i/ 2 ^ < \\d x j3\\ H ^.+ P for all p > 0, for 
some absolute constant p large enough, because 

g + fh 2 = h 2 {g + C 2 )(l + 2 ) ), h := ( 1 + d xP)~ 1 , 

and then use Taylor expansion for the square root of the last factor. The second 
component p\ is given by formula (C.22). By Taylor expansion, 

ko© - n«rk| < a + ifiri K© - !i«C| < (i + leiri 

so that we calculate 

Pi = il{d X xf3){l + d x (3)-^\C\~Hx(C)+r, 
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(C.25) 











where the remainder r satisfies || Op(r)|| jC ^ (J j^+ 3 / 2 ) < \\d x p\\ H ^+ P for all p > 0 , for 

some p large enough. Assume that ||3 r /3||i2> < ||/0 ||#m (this bound holds for the 
choice of P we make below). By (C.23), (C.24), (C.25), we have 

^(r)L^(r)- 1 = (1 + dJ)~?L + §(d xx 0)(l + dJ)~l\D x \~^d x + 72 0 , u (C.26) 

where 72o,i is defined by difference and it satisfies ||72o,i \\c(w i ,Hi i ) 1$ \\dxP\\Hv+p for 
all p > 0, for some p large enough. With similar calculations, one proves that for 
any r € M 

^(r)|£) x r®(r )- 1 = (1 + dJ)~ r \D x \ r + 72 0 , 2 (C.27) 

where 72-0,2 is defined by difference and it satisfies ||72o,2|| J c(.H>,.ff* i - r + 1 ) < ||d x /3|| h»+p- 


C.4. Conjugation of P. We conjugate the operator in (C.l) by VEb := ’F(l) = 
4/(r)| r= i. From symbolic calculus it follows that 

11 3 1 

LzcL? = cL — — (d x c) d x \D x | 2 + 72o,3, (C.28) 

where 72 o ; 3 is defined by difference and it satisfies ||7^.o,3 II + i < ||<9 x c||^ M + P for 

all p > 0, for some p large enough. We recall that c — 1 is small, and therefore d x c 
is small. By definition (see (C.4),(C.5)), and recalling that P\ T =\ = Pi, P\r=l = Pi, 
we directly calculate 

’FiSt’Fj" 1 = d t + aid x + n, 4'1^4'j' 1 = a 2 d x + r 2 , 

where 

ai(t,x) := (d t Pi)(t,x + Pi(t,x)), a 2 (t, x) := (1 + d x pi(t, x)) -1 , (C.29) 


and 

ri(t,x) := ~(dt x Pi)(t,x + Pi(t,x)) (1 + d x Pi(t,x)), 

r 2 (t,x ) := -(1 +d x p 1 (t,x)) {d xx Pi)(t,x + f3i(t,x)). 

The conjugate of any multiplication operator h 1 —> ah is the multiplication operator 
h i->- ( Ba)h , 

Tiaikj" 1 = (Ba), ( Ba)(t , x) := a(t, x + p± (t, x )). 

Thus 

'FiPTj" 1 = d t + a 3 d x + ia^L + ia 5 d x \D x \~% + R 3 

where 

a 3 := ai + (BV)a 2 , a 4 := (Bc)( 1 + d x P )~^, 

as := -^| - | (Bc)(l + d x (3)~i(d xx p) + (£?(d x c))(l + <9x/?)~5 j ; 

3 11 

#3 := n + (BV)r 2 + z(77c)72 0 ,i + -z-(-B<9 x c)7- 2 (l + <9 x /3)^|.D x | - 5 + 72 0 ,2 (C.30) 

+ ^i 72 0 , 3 W 1 - 1 + ^ii? 2 ^ 1 , 

72o.i is defined in (C.26) with r = 1 , 72o, 2 is defined in (C.27) with r = 1, r = —1/2, 
and 72 o, 3 is defined in (C.28). The remainder R 3 is of order zero and it is estimated 
in Lemma C.l. Moreover, as it is immediate to verify, 05 = — |< 9 x a 4 . We choose 
Pi, Pi such that the highest order coefficient 04 is independent of x. This means 

04 ( 7 , x) = c(t, x + Pi(t, x)) (1 + d x Pi(t, x))“2 = m(t) Vx € T, 

66 


(C.31) 


for some function m[t) independent of x. Applying the inverse diffeomorphism, this 
is equivalent to 

c(t, x) (1 + d x /3i(t, x)) 2 = m(t) Vx € T. 

2 2 

This implies l + d x (3i(t,x) = m{t) 3c(t, x)~s , which, after an integration in dx, gives 

-j^ /» _3 

m{t) = (— j c(t,x)~^dx^ 

Hence m in (C.31) is determined. We fix j3\ as 

Pi(t,x) = d~ l [m{t)^c{t,x)~^ - 1], 

and then we fix j3(r,t,x) := r/3i(t,x). As a consequence, /3(r,t,y ), j3i are also 
determined. Since <74 (i, x) = m(t) is independent of x, it follows that a§ = — \d x a± = 
0 (as it was natural to expect, because the vector field in P is anti-selfadjoint and 
the transformation T preserves this structure). We have conjugated P to 

Pi := ’J'iH’I '/ 1 = d t + im(t)L + a 3 d x + R 3 . (C.32) 

We underline that the coefficient m(t) is a function of time, independent of space. 

Lemma C.l. There exists a universal constant 5 q € (0,1) such that if 

|| c(t) - l|| L oo < So 

then \\d x /3i(t)\\ L °° + \\d x Pi(t)\\ L °° < 1/2 and 

||9 x /3i (t )||+ \\d x Pi(t)\\wiJ.,oo < Cfj,\\c(t) — l||w>.°° V// > 0 

for some positive constant C fl depending only on p,. As a consequence, Ti(f ) _1 

are bounded transformations of H^{ T), with 

ll^iCOII c(h») + ll^ r i( i ) _ 1 |k(J?M) < Cfj,( 1 + ||c(t) — 1|| m) V/i > 0. 

Moreover \m(t) — 1| < C\\c(t) — l||#i, 

||a 3 (i)||/^ < C^(\\c(t) - 1 ||+ \\d t c(t)\\ H »-i) Mp > 1. 

The remainder R^ft) maps L 2 ( T) into itself, with 

\\P3(t)\\c{L 2 ) < C{\\ c (t) ~ 1|| H r + \\V{t)\\L°° + ||^tc(t) ||z,°o + \\R 2 (t)\\c(L 2 )) 1 

and, for all p > 1/2, i? 3 (f) also maps H^{ T) into itself, with 

11^3(01 < Cfj,(\\c(t) — 1|| H»+r + ||V(i)||i?^ + ||^c(t)||^ + 11^2(0 ||/:(h/*)) i 

where r > 0 is a universal constant. 

Proof. The estimates follow from the explicit formulas above, the usual estimates for 
the composition of functions (see, e.g., Appendix B in [7]) and Sobolev estimates for 
pseudo-differential operators (see (C.14)). The estimate of the pseudo-differential 
remainder term is the reason for which r further space-derivatives are required on c. 
The term dtc appears only in a\ and r\. The term V appears only in a 3 and R 3 where 
it is explicitly written, and nowhere else. The operator R 2 only appears in R 3 in the 
term All the other terms depend only on c and its space-derivatives. □ 
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C.5. Reparametrization of time. Now we want to replace the coefficient m(t) in 
(C.32) with a constant coefficient. We consider a diffeomorphism of the time interval 

: [0,T] ^ [0,Ti], V’(O) = 0, ip(T) = T 1 , ip'(t)> 0, 

where T\ > 0 has to be determined. We consider the pull-back ?/>* defined as 
(y>*/i)(i, x) := h(ip(t),x), and similar for its inverse if . Then we calculate the 
conjugate 

^{dt + im(t)L )V>* = if'{if~ l {t))d t + zm(^ _1 (t))L. 

The two time-dependent coefficients are equal if m(t) = if'{t) for all t E [0,T]. 

We define 

if(t) := f m(s)ds, T\ := f m(t)dt, p(t) := miif -1 (t)). (C.33) 

Jo Jo 

Since \m — 1| is small, then the ratio T\/T is close to 1, and also iff ft) is close to 1 
for all t. We have conjugate 

iff l P\if* = p(t)P 2 , P '2 '■= dt + iL + a 6 d x + R 4 , (C.34) 

where 

a 6 {t,x) := a , P 4 ■■= -^-if^R^if* (C.35) 

p{t) p{t) 

(and, more explicitly, (iff 1 R- 3 if*)(t) = Roi’f ’ -1 (£)))• Now the coefficient of the high¬ 
est order term L is constant. 

C.6. Translation of the space variable. The goal of this section is to eliminate 
the space-average of the coefficient a&(t, x) in front of d x . Consider a time-dependent 
change of the space variable which is simply a translation, 

y = tp(t, x) = x + p(t) <s> x = ^ _1 (t,y) = y -p(t), 

and its pull-back (ip*h)(t,x) = h(t,(p(t,x )) = h(t,x + p(t)), and similarly for ip~ l . 
Thus = dt + p'(t)d x , and </?* commutes with every Fourier multiplier like 

d x , \D x \ r ,L. We calculate the conjugate 

Pi '■= p * 1 P 2 ( P* = dt + iL + ajd x + R5, 

where 

a 7 := p'(t) + (pf l a & ), R 5 := . (C.36) 

Since preserve the space average, we fix 

p(t) :=- [ [ ae(s,x)dxds. (C.37) 

27r Jo J T 

It follows that fjaj(t,x)dx = 0 for all t € [0, Tj]. Note that 93* commutes with the 
multiplication operator h i->- p(t)h , because p(t) is independent of x. Moreover, by 
the change of time variable s = if(t), ds = m(t)dt in the integral, we get 

p(T\) = - [ f 06 (s,x)dxds = - [ [ as(t, x) dxdt. (C.38) 

27T J o J T 27r J 0 J T 

Proof of Proposition 5.1 concluded. The composition <f> := of the previ¬ 
ous three transformations conjugates P = Also note that <l> _1 (/w) = 

for all u. The transformation ’Ll is estimated in Lemma C.l. The estimates 
for if*,<p* are straightforward. Finally, rename W := a 7 and R 3 := R§. □ 

Notation. In the following Proposition we use the shorter notation ||o||t,x to denote 
the C°([0, T];X) norm of any u, with X = L 2 (T), L°°(T), iP*(T), £(L 2 (T)), etc. 
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Proposition C.2. Assume the hypotheses of Proposition 5.1. 

(i) (Regularity). In addition, let p > 1/2, let ||c— 1|| r,Hr <K< oo, and let 

A// := ||c— 1|| T,Hr+r + \\V\\t,Hv + ||^c||t,^ + 11-^21| T,C(Hr) < OO. 

Then R 3 maps C°([0,Ti]; H^(T)) into itself, with 

II-R3IIT1 ,C(H?) < (C.39) 

for some constant depending on p, K. For p > 1, 

< Cn (||c 1|| T,Hr + \\d t c\\ TjH v-i + ||R||t,.H>) (C.40) 

and 

< C^\\c\\T,H»\\u\\T,Hr, || ^~ l u\\r,Hr < Cfj\\c\\T,Hr 

for all u = u(t,x), for some constant C / depending only on p. 

(ii) (Stability). Consider another triple {c!, V ', R 2 ) such that d also satisfies (5.1), 
and A/o < 00 also for (d, V', R 2 ). Let <$', 4'' 1 , </4 , tji 1 *, T[, W', R' 3 be the corresponding 
objects obtained for the triple (d,V', R 2 ). Then for all u G L 2 ( T), all t G [0, T], 

||4'i(t)ri-4'' 1 (t)u|| i 2 + ||4'i(t) _1 ri-'I' , 1 (t) _1 tt|| i 2 < C\\c(t) - c'(t)\\ L 2 \\u\\ H i . (C.41) 

Let X := T\ jT[, and let T be the time-rescaling operator ( Tv)(t,x ) := v{\t,x). Then 
for all p > 0, all v = v(t, x), 

\\if*(p*v — ifW*{Tv)\\ T,Hr < CT^WdtvW^^r + IMlTi,/fe+i) A 0 (C.42) 

II ip'c 1 'ip'~ 1 v - T(^“V“ 1 'r)|| T ' ) ^ < CT(\\d t v\\ T ,Hr + \\v\\ TtH „+i)A 0 (C.43) 

where 

Ao := Il c — c'Ht,// 1 + \\dtc — dfd, V — \ '\\t,l 2 - 

Also, 

|1 — A| + ||m — m'\\co([ 0 j T]) — R|l c — c '|| (C.44) 

and, if 

M{x) :={1 + dJ 1 (T,x-p{T 1 ))}h, M’{x) := {1 + dJ[(T, x - p'(T[))}*, 

then 

||M — M r ||l°o(t) < C(||c — c / ||y i ^2 + ||0iC — dfd, V — bplr,/, 2 )- (C.45) 

For p > 1, if 

||c — 1|| T,Hr+i + \\dtc\\T,He- + || < 9| c ||:r,/// i -i + P IIt,.H>+i + \\dtV\\T,Hr < 1, (C.46) 
and (C.46) also holds for d ,V', then 

\\w' -TW\\ t{jH , < C^Wc-d\\ T ,H^ + \\dtc-d t d\\ T , H r-i + \\v -V'\\t,h^. (C.47) 
Moreover, if 

ll c — 1 ||t,//h-i + ll^i c llr,^+ 1 + II^ c IIt,l 2 + P Hr,// 1 + IPP IIt,l 2 

+ 11-^2 \\T,c(m)nC{L 2 ) + ll^^2||r,£(L 2 ) < 1 (C.48) 

and (C.48) also holds for d,V',R 2 , then 

11-^3 — (T-^3)|| t[,c(l 2 ) < C(ll c — P||t.j/h-i + ll^c — ^c'llr^i 

+ l|R - v'W t ,// 1 + II-R2 - ll' 2 \\ TC ( L 2)'). (c.49) 

Proof. To prove statement (ii) we make repeatedly use of triangular inequality and 
explicit formulas. In particular, to estimate p( , ip(Xt)) — p'(if'(t )), we use explicit 
formulas similar to (C.38). To estimate R ' 5 — (TR 5 ) we note that the rescaled 
operator (TR 5 ) is the composition T-R5T , and then we also use (C.42)-(C.43). 
Remember that we have renamed W := 07 and R 3 := R§. □ 
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